The Sinocyclocheilus cavefish genome provides insights into cave adaptation

Decoding the fish genome opens a new era in important trait research and molecular breeding in China

Aquaculture represents the fastest-growing global food production sector, as it has become an essential component of the global food supply. China has the world's largest aquaculture industry in terms of production volume. However, the sustainable development of fish culture is hindered by several concerns, including germplasm degradation and disease outbreaks. The practice of genomic breeding, which relies heavily on genome information and genotypephenotype relationships, has significant potential for increasing the efficiency of aquaculture production. In 2014, the completion of the genome sequencing and annotation of the Chinese tongue sole signified the beginning of the fish genomics era in China. Since then, domestic researchers have made dramatic progress in functional genomic studies. To date, the genomes of more than 60 species of fish in China have been assembled and annotated. Based on these reference genomes, evolutionary, comparative, and functional genomic studies have revolutionized our understanding of a wide range of biologically and economically important traits of fishes, including growth and development, sex determination, disease resistance, metamorphosis, and pigmentation. Furthermore, genomic tools and breeding techniques such as SNP arrays, genomic selection, and genome editing have greatly accelerated genetic improvement through the incorporation of functional genomic information into breeding activities. This review aims to summarize the current status, advances, and perspectives of the genome resources, genomic study of important traits, and genomic breeding techniques of fish in China. The review will provide aquaculture researchers, fish breeders, and farmers with updated information concerning fish genomic research and breeding technology. The summary will help to promote the genetic improvement of production traits and thus will support the sustainable development of fish aquaculture.

Genetic and selective constraints on the optimization of gene product diversity

RNA and protein expressed from the same gene can have diverse isoforms due to various post-transcriptional and post-translational modifications. For the vast majority of alternative isoforms, It is unknown whether they are adaptive or simply biological noise. As we cannot experimentally probe the function of each isoform, we can ask whether the distribution of isoforms across genes and across species is consistent with expectations from different evolutionary processes. However, there is currently no theoretical framework that can generate such predictions. To address this, we developed a mathematical model where isoform abundances are determined collectively by cis-acting loci, trans-acting factors, gene expression levels, and isoform decay rates to predict isoform abundance distributions across species and genes in the face of mutation, genetic drift, and selection. We found that factors beyond selection, such as effective population size and the number of cis-acting loci, significantly influence evolutionary outcomes. Notably, suboptimal phenotypes are more likely to evolve when the population is small and/or when the number of cis-loci is large. We also explored scenarios where modification processes have both beneficial and detrimental effects, revealing a non-monotonic relationship between effective population size and optimization, demonstrating how opposing selection pressures on cis- and trans-acting loci can constrain the optimization of gene product diversity. As a demonstration of the power of our theory, we compared the expected distribution of A-to-I RNA editing levels in coleoids and found this to be largely consistent with non-adaptive explanations.

Exudates of Microcystis aeruginosa on oxidative stress and inflammatory responses in gills of Sinocyclocheilus grahami

Early cyanobacterial blooms studies observed that exposure to blue-green algae led to fish gills impairment. The objective of this work was to evaluate the toxic mechanisms of exudates of Microcystis aeruginosa (MaE) on fish gills. In this study, the toxic mechanism of MaE (2×106 cells/mL) and one of its main components phytosphingosine (PHS) with two concentrations 2.9 ng/mL and 145 ng/mL were conducted by integrating histopathology, biochemical biomarkers, and transcriptomics techniques in Sinocyclocheilus grahami (S. grahami) for 96 h exposure. Damaged gill tissue with epithelial hyperplasia and hypertrophy, remarkable Na+/K+-ATPase (NKA) enzyme activity, disrupted the redox homeostats including lipid peroxidation and inflammatory responses were observed in the fish of MaE exposure group. Compare to MaE exposure, two concentrations of PHS exposure appeared to be a trend of lower degree of tissue damage, NKA activity and oxidative stress, but induced obviously lipid metabolism disorder with higher triglycerides, total cholesterol and total bile acid, which might be responsible for inflammation responses in fish gill. By transcriptome analysis, MaE exposure were primarily enriched in pathways related to gill function and immune response. PHS exposure, with higher number of differentially expressed genes (DEGs), were enriched in Toll-like receptor (TLR), Mitogen-Activated Protein Kinase (MAPK) and NOD-like receptor protein 3 (NLRP3) pathways. We concluded that MaE and PHS were induced the inflammatory responses, with oxidative stress-induced inflammation for MaE exposure but lipid metabolism disorder-induced inflammation for PHS exposure. The present study provided two toxin-induced gill inflammation response pathways under cyanobacterial blooms, which could be a scientific basis for the ecological and health risk assessment in the aquatic environment.

Chromosome-scale genome assemblies of sexually dimorphic male and female Acrossocheilus fasciatus

Acrossocheilus fasciatus is a stream-dwelling fish species of the Barbinae subfamily. It is valued for its colorfully striped appearance and delicious meat. This species is also characterized by apparent sexual dimorphism and toxic ovum. Biology and aquaculture researches of A. fasciatus are hindered by the lack of a high-quality reference genome. Here, we report chromosome-level genome assemblies of the male and female A. fasciatus. The HiFi-only genome assemblies for both female and male individuals were 899.13 Mb (N50 length of 32.58 Mb) and 885.68 Mb (N50 length of 33.06 Mb), respectively. Notably, a substantial proportion of the assembled sequences, accounting for 96.15% and 98.35% for female and male genomes, respectively, were successfully anchored onto 25 chromosomes utilizing Hi-C data. We annotated the female assembly as a reference genome and identified a total of 400.62 Mb (44.56%) repetitive sequences, 27,392 protein-coding genes, and 35,869 ncRNAs. The high-quality male and female reference genomes will provide genomic resources for developing sex-specific molecular markers, inform single-sex breeding, and elucidate genetic mechanisms of sexual dimorphism.

Cave beetle lineages gained genes before going down under: An example of repeated genomic exaptation?

The adaptation of animals to subterranean habitats like caves and aquifers stereotypically leads to dramatic trait-loss consequences like the lack of eyes and body pigmentation. These body plan regression trends are expected to be tied to gene loss as well. Indeed, previous studies documented the degeneration of vision genes in obligate cave dwellers. Contradicting this picture, the first broad-scale comparative transcriptome-wide study of gene content evolution in separate subterranean Australian and Mediterranean beetle clades unearthed evidence of global gene gain and retention. This suggests that the transition to cave life may be more contingent on gene repertoire expansion than contraction. Future studies, however, will need to examine how much the observed patterns of gene content evolution reflect subfunctionalization and fitness-securing genetic redundancy outcomes following gene duplication as opposed to adaptive trajectories.

Prevotella enterotype associates with diets supporting acidic faecal pH and production of propionic acid by microbiota

Metabolism of dietary fibres by colon microbiota plays an important role for human health. Personal data from a nutrition study (57 subjects) were analysed to elucidate quantitative associations between the diet, faecal microbiome, organic acid concentrations and pH. Ratios of the predominant acids acetate, butyrate and propionate ranged from 1:0.67:0.27 to 1:0.17:0.36. Pectin-rich diets resulted in higher faecal acetate concentrations. Negative correlation between faecal pH and BSS was observed. Higher faecal pH and lower acid concentrations were related to the higher abundance of amino acid degrading Clostridium, Odoribacter and Eubacterium coprostanoligenes, which are weak carbohydrate fermenting taxa. Propionic acid correlated especially to high abundance of Prevotella and low abundance of proteobacteria. The acetate to propionate ratio of the Prevotella enterotype was about half of that of the Bacteroides enterotype. Based on the results we suggest the measurement of faecal pH and organic acid composition for research and diagnostic purposes.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

Highly dynamic evolution of the chemosensory system driven by gene gain and loss across subterranean beetles

Chemical cues in subterranean habitats differ highly from those on the surface due to the contrasting environmental conditions, such as absolute darkness, high humidity or food scarcity. Subterranean animals underwent changes to their sensory systems to facilitate the perception of essential stimuli for underground lifestyles. Despite representing unique systems to understand biological adaptation, the genomic basis of chemosensation across cave-dwelling species remains unexplored from a macroevolutionary perspective. Here, we explore the evolution of chemoreception in three beetle tribes that underwent at least six independent transitions to the underground, through a phylogenomics spyglass. Our findings suggest that the chemosensory gene repertoire varies dramatically between species. Overall, no parallel changes in the net rate of evolution of chemosensory gene families were detected prior, during, or after the habitat shift among subterranean lineages. Contrarily, we found evidence of lineage-specific changes within surface and subterranean lineages. However, our results reveal key duplications and losses shared between some of the lineages transitioning to the underground, including the loss of sugar receptors and gene duplications of the highly conserved ionotropic receptors IR25a and IR8a, involved in thermal and humidity sensing among other olfactory roles in insects. These duplications were detected both in independent subterranean lineages and their surface relatives, suggesting parallel evolution of these genes across lineages giving rise to cave-dwelling species. Overall, our results shed light on the genomic basis of chemoreception in subterranean beetles and contribute to our understanding of the genomic underpinnings of adaptation to the subterranean lifestyle at a macroevolutionary scale.

Convergent gene losses and pseudogenizations in multiple lineages of stomachless fishes

The regressive evolution of independent lineages often results in convergent phenotypes. Several teleost groups display secondary loss of the stomach, and four gastric genes, atp4a, atp4b, pgc, and pga2 have been co-deleted in agastric (stomachless) fish. Analyses of genotypic convergence among agastric fishes showed that four genes, slc26a9, kcne2, cldn18a, and vsig1, were co-deleted or pseudogenized in most agastric fishes of the four major groups. kcne2 and vsig1 were also deleted or pseudogenized in the agastric monotreme echidna and platypus, respectively. In the stomachs of sticklebacks, these genes are expressed in gastric gland cells or surface epithelial cells. An ohnolog of cldn18 was retained in some agastric teleosts but exhibited an increased non-synonymous substitution when compared with gastric species. These results revealed novel convergent gene losses at multiple loci among the four major groups of agastric fish, as well as a single gene loss in the echidna and platypus.

A New Species of the Blind Cave Loach Genus Protocobitis (Cypriniformes: Cobitidae), Protocobitis longicostatus sp. nov., from Guangxi, China

Protocobitis species are typical cave-dwelling fish, exhibiting distinctive morphological adaptations such as colorless body, lack of eyes, and reduced scales and ribs in response to their extreme cave habitats. Distinct from the recorded species, P. anteroventris, P. polylepis, and P. typhlops, a new species, Protocobitis longicostatus sp. nov., is described from Guangxi Zhuang Autonomous Region, China. Protocobitis longicostatus sp. nov. can easily be distinguished from all known congeners by the following characteristics: whole body covered by scales except head, 12 branched caudal fin rays, and long ribs. These species face threats from habitat degradation, hydrological changes, and environmental pollution. Thus, the conservation of cavefish in China has become an urgent issue.

Full-length transcriptome reveals the circularly polarized light response-related molecular genetic characteristics of Oratosquilla oratoria

The mantis shrimp is the only animal that can recognize circularly polarized light (CPL), but its molecular genetic characteristics are unclear. Multi-tissue level full-length (FL) transcriptome sequencing of Oratosquilla oratoria, a representative widely distributed mantis shrimp, was performed in the present study. We used comparative transcriptomics to explore the critical genes of O. oratoria selected by CPL and the GNβ gene associated with CPL signal transduction was hypothesized to be positively selected. Furthermore, the FL transcriptomes of O. oratoria compound eyes under five light conditions were sequenced and used to detect alternative splicing (AS). The ASs associated with CPL recognition mainly occurred in the LWS, ARR and TRPC regions. The number of FL transcripts with AS events and annotation information also provided evidence that O. oratoria could recognize LCPL. Additionally, 51 sequences belonging to the LWS, UV and Peropsin gene families were identified based on conserved 7tm domains. The LWS, UV and Peropsin opsins have similar 3D structures with seven domains across the cell membrane and conserved KSLRTPSN, DRY, and QAKK motifs. In conclusion, these results are undoubtedly valuable for perfecting the vision theory of O. oratoria and other mantis shrimp.

Gut microbiota reflect adaptation of cave-dwelling tadpoles to resource scarcity

Gut microbiota are significant to the host's nutrition and provide a flexible way for the host to adapt to extreme environments. However, whether gut microbiota help the host to colonize caves, a resource-limited environment, remains unknown. The nonobligate cave frog Oreolalax rhodostigmatus completes its metamorphosis within caves for 3-5 years before foraging outside. Their tadpoles are occasionally removed from the caves by floods and utilize outside resources, providing a contrast to the cave-dwelling population. For both cave and outside tadpoles, the development-related reduction in their growth rate and gut length during prometamorphosis coincided with a shift in their gut microbiota, which was characterized by decreased Lactobacillus and Cellulosilyticum and Proteocatella in the cave and outside individuals, respectively. The proportion of these three genera was significantly higher in the gut microbiota of cave-dwelling individuals compared with those outside. The cave-dwellers' gut microbiota harbored more abundant fibrolytic, glycolytic, and fermentative enzymes and yielded more short-chain fatty acids, potentially benefitting the host's nutrition. Experimentally depriving the animals of food resulted in gut atrophy for the individuals collected outside the cave, but not for those from inside the cave. Imitating food scarcity reproduced some major microbial features (e.g. abundant Proteocatella and fermentative genes) of the field-collected cave individuals, indicating an association between the cave-associated gut microbiota and resource scarcity. Overall, the gut microbiota may reflect the adaptation of O. rhodostigmatus tadpoles to resource-limited environments. This extends our understanding of the role of gut microbiota in the adaptation of animals to extreme environments.

The two-way immunotoxicity in native fish induced by exudates of Microcystis aeruginosa: Immunostimulation and immunosuppression

Secondary metabolites of cyanobacterial blooms have caused serious risks to aquatic animals. The immune system is an important barrier for fish against pollutants in aquatic systems. The immunetoxic mechanism of the exudates of Microcystis aeruginosa (MaE) on fish was lacking due to the complex components of MaE. In this project, Sinocyclocheilus grahami was used as the model to study the immunotoxic effects of MaE and PHS (one of the main components of the MaE) in fish. The immunosuppression effects of MaE are mainly in, decreased head-kindey index, damaged tissue structure of head-kidney and downregulated NF-κB, IL-1β. PHS induce immunostimulation via, increasing spleen index, apparently increasing leucocytes, increasing the IgM and lysozyme levels in serum and skin mucus, upregulating protease in skin mucus, increasing pro-immunologic factors (IL-1β, IL-6, IL-8, IL-10, TNF-α and NF-κB), probably activating the TLRs/NF-κB, MAPK, FoxO1 and PPARγ signaling pathways. Therefore, our research identified potential data gaps that how the exudates of cyanobacteria induces immunostimulation and immunosuppression from immune organs level to skin mucus to blood cells to inflammatory factors to potential molecular initiating event of MaE and PHS. Further research is needed to obtain a deeper view of the molecular mechanisms involved in MaE and PHS immunotoxicity and its consequences in long-time exposures.

Expansion and contraction of lake basin shape the genetic structure of Sinocyclocheilus (Osteichthyes: Cypriniformes: Cyprinidae) populations in Central Yunnan, China

Geological events can strongly affect the genetic structures and differentiation of fish populations. Especially, as an endemic fish of the genus Sinocyclocheilus in the Yunnan-Guizhou Plateau, the effects of key geological events on the distributions and genetic structures remain poorly understood. Examining the phylogeographic patterns of Sinocyclocheilus fishes can be useful for elucidating the spatio-temporal dynamics of their population size, dispersal history and extent of geographical isolation, thereby providing a theoretical basis for their protection. Here, we used single nucleotide polymorphisms (SNP) method to investigate the phylogeographic patterns of Sinocyclocheilus fishes. Our analysis supports the endemicity of Sinocyclocheilus, but the samples of different regions of Sinocyclocheilus contain multiple ancestral components, which displayed more admixed and diversified genetic components, this may be due to the polymorphism of the ancestors themselves, or gene infiltration caused by hybridization between adjacent species of Sinocyclocheilus. We estimate that the most recent common ancestor (MRCA) of Sinocyclocheilus fish in the Central Yunnan Basin at approximately 3.75~3.11 Ma, and infer that the evolution of Sinocyclocheilus in the central Yunnan Basin is closely related to the formation of plateau lakes (around 4.0~0.02 Ma), and identifies the formation of Dianchi Lake and Fuxian Lake as key geological events shaping Sinocyclocheilus population structure. It is also the first time to prove that the altitude change has a great influence on the genetic variation among the populations of Sinocyclocheilus.

Gene Characterization of Nocturnin Paralogues in Goldfish: Full Coding Sequences, Structure, Phylogeny and Tissue Expression

The aim of this work is the full characterization of all the nocturnin (noc) paralogues expressed in a teleost, the goldfish. An in silico analysis of the evolutive origin of noc in Osteichthyes is performed, including the splicing variants and new paralogues appearing after teleostean 3R genomic duplication and the cyprinine 4Rc. After sequencing the full-length mRNA of goldfish, we obtained two isoforms for noc-a (noc-aa and noc-ab) with two splice variants (I and II), and only one for noc-b (noc-bb) with two transcripts (II and III). Using the splicing variant II, the prediction of the secondary and tertiary structures renders a well-conserved 3D distribution of four α-helices and nine β-sheets in the three noc isoforms. A synteny analysis based on the localization of noc genes in the patrilineal or matrilineal subgenomes and a phylogenetic tree of protein sequences were accomplished to stablish a classification and a long-lasting nomenclature of noc in goldfish, and valid to be extrapolated to allotetraploid Cyprininae. Finally, both goldfish and zebrafish showed a broad tissue expression of all the noc paralogues. Moreover, the enriched expression of specific paralogues in some tissues argues in favour of neo- or subfunctionalization.

Re-identification and characterization of grass carp Ctenopharyngodon idella TLR20

Toll-like receptors (TLRs) play a crucial role in the recognition of microbial-associated molecular patterns in the innate immune system. Fish TLRs have undergone significant gene expansion to adapt to complex aquatic environments. Among them, TLR20 from the TLR11 family actively responds to viral and bacterial invasions. Previous studies have reported two TLR20s in grass carp (Ctenopharyngodon idella), and in this study, we revised this conclusion. Based on the latest grass carp genome, we identified a new TLR20 member. These three TLR20s are arranged in tandem on chromosome 9, indicating that they are generated by gene duplication events. They were renamed CiTLR20.1 to CiTLR20.3 based on their chromosomal positions. The CiTLR20s in C. idella exhibit higher similarities with those in Danio rerio, Cyprinus carpio, and Megalobrama amblycephala, and lower similarities with those in other distantly related fish species. Selective pressure analysis revealed low conservation and negative evolution of TLR20s during evolution. The 3D structures of the three TLR20s showed significant differences, reflecting functional variations and different downstream adaptor molecule recruitment. Transcriptome data revealed tissue distribution differences of TLR20s, with TLR20.1 showing relatively low expression levels in all the tissues, while TLR20.2 and TLR20.3 showed higher expression in the head kidney, spleen, and gill. Additionally, TLR20.2 and TLR20.3 actively responded to GCRV-II infection, with higher upregulation of TLR20.2 in response to Aeromonas hydrophila challenge. In conclusion, this study corrected the number of grass carp TLR20 members and analyzed TLR20 from an evolutionary and structural perspective, exploring its role in antiviral and antibacterial defense. This study provides reference for future research on fish TLR20.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

A long-ignored unique ecosystem of cavefishes in the southern karst: achievements, challenges, and prospects

Cavefishes represent a taxon that has experienced natural selection pressures. This paper summarizes the results with respect to the taxonomy, diversity, phylogeny, and adaptation aspects of cavefishes research. It showed that: 1) These studies suggest that cavefishes play important roles in the study of geologic history and adaptation to extreme environments, but the mechanisms involved 168 species of cavefishes belonging to 17 genera, four families, and two orders have been recorded in China. Meanwhile, more new species are being discovered recently, and the species diversity of cavefishes are still underestimated, indicating the need to strengthen the survey in field. 2) The biogeography of cavefishes have focused on Sinocyclocheilus and Triplophysa, that have helped understand the geomorphology of karst areas in southern China and the spatial pattern of species diversity. These studies revealed the influences of evolution and geological history in Sinocyclocheilus, but there are still many species that have not been studied accordingly. 3) Some adaptive mechanistic studies have been conducted on cavefishes, primarily focusing on eye and body color degradation and energy metabolism in the genus Sinocyclocheilus to reveal adaptive mechanisms in the dark environment. 4) The IUCN list of protected cavefishes species in China only includes 21 species. The List of Key Protected Wild Animals for 2021 includes all species of Sinocyclocheilus as National Class II.It is necessary to strengthen the research on the biodiversity and adaptation and need consider the conservation actions for cavefishes.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

Metagenomic analysis reveals hidden links between gut microbes and habitat adaptation among cave and surface dwelling Sinocyclocheilus species

Intestinal microbes are closely related to vital host functions such as digestion and nutrient absorption, which play important roles in enhancing host adaptability. As a natural "laboratory", caves provide an outstanding model for understanding the significance of gut microbes and feeding habits in the habitat adaptability of hosts. However, research on the relationship between gut microbes, feeding habits, and the adaptability of troglobites remains insufficient. In this study, we compared the characteristics of the intestinal microbes of Sinocyclocheilus cavefish and surface fish and further established the relationship between intestinal and habitat microbes. Furthermore, we conducted environmental DNA (eDNA) (metabarcoding) analysis of environmental samples to clarify the composition of potential food resources in the habitats of the Sinocyclocheilus cavefish and surface fish. Results showed that the structure of the Sinocyclocheilus gut microbes was more related to ecological type (habitat type) than phylogenetic relationships. While horizontal transfer of habitat microbes was a source of gut microbes, hosts also showed strong selection for inherent microbes as dominant microorganisms. Differences in the composition and structure of gut microbes, especially dominant microbes, may enhance the adaptability of the two Sinocyclocheilus fish types from the perspectives of food intake, nutrient utilization, and harmful substance metabolism, suggesting that food resources, predation patterns, intestinal flora, digestive and absorptive capacity, and feeding habits and preferences are linked to habitat adaptability. These results should facilitate our understanding of the significance of fish gut microbes to habitat adaptation and provide a new perspective for studying the adaptive mechanisms of cavefish.

Mitochondrial phylogeography and molecular evolution of the rhodopsin visual pigment in troglobitic populations of Astyanax mexicanus (De Filippi, 1853)

Cave-adapted animals provide a unique opportunity to study the evolutionary mechanisms underlying phenotypic, metabolic, behavioral, and genetic evolution in response to cave environments. The Mexican tetra ( Astyanax mexicanus) is considered a unique model system as it shows both surface and cave-dwelling morphs. To date, at least 33 different cave populations have been identified, with phylogenetic studies suggesting an origin from at least two independent surface lineages, thereby providing a unique opportunity to study parallel evolution. In the present study, we carried out the most exhaustive phylogeographic study of A. mexicanus to date, including cave and surface localities, using two mitochondrial markers (cytochrome b (cyt b) and cytochrome c oxidase subunit I ( COI)) and nuclear rhodopsin visual pigment ( rho). Additionally, we inferred the molecular evolution of rho within the two contrasting environments (cave and surface) and across three geographic regions (Sierra de El Abra, Sierra de Guatemala, and Micos). In total, 267 individuals were sequenced for the two mitochondrial fragments and 268 individuals were sequenced for the rho visual pigment from 22 cave and 46 surface populations. Phylogeographic results based on the mitochondrial data supported the two-lineage hypothesis, except for the Pachón and Chica caves, whose introgression has been largely documented. The Sierra de El Abra region depicted the largest genetic diversity, followed by the Sierra de Guatemala region. Regarding the phylogeographic patterns of rho, we recovered exclusive haplogroups for the Sierra de El Abra (Haplogroup I) and Sierra de Guatemala regions (Haplogroup IV). Moreover, a 544 bp deletion in the rho gene was observed in the Escondido cave population from Sierra de Guatemala, reducing the protein from seven to three intramembrane domains. This change may produce a loss-of-function (LOF) but requires further investigation. Regarding nonsynonymous ( dN) and synonymous ( dS) substitution rates (omega values ω), our results revealed the prevailing influence of purifying selection upon the rho pigment for both cave and surface populations (ω<1), but relaxation at the El Abra region. Notably, in contrast to the other two regions, we observed an increase in the number of dN mutations for Sierra de El Abra. However, given that a LOF was exclusively identified in the Sierra de Guatemala region, we cannot dismiss the possibility of a pleiotropic effect on the Rho protein.

RNA polymerase II pausing is essential during spermatogenesis for appropriate gene expression and completion of meiosis

Male germ cell development requires precise regulation of gene activity in a cell-type and stage-specific manner, with perturbations in gene expression during spermatogenesis associated with infertility. Here, we use steady-state, nascent and single-cell RNA sequencing strategies to comprehensively characterize gene expression across male germ cell populations, to dissect the mechanisms of gene control and provide new insights towards therapy. We discover a requirement for pausing of RNA Polymerase II (Pol II) at the earliest stages of sperm differentiation to establish the landscape of gene activity across development. Accordingly, genetic knockout of the Pol II pause-inducing factor NELF in immature germ cells blocks differentiation to mature spermatids. Further, we uncover unanticipated roles for Pol II pausing in the regulation of meiosis during spermatogenesis, with the presence of paused Pol II associated with double strand break formation by SPO11, and disruption of SPO11 expression in germ cells lacking NELF.

Molecular evolution of the hemoglobin gene family across vertebrates

Adaptation to various altitudes and oxygen levels is a major aspect of vertebrate evolution. Hemoglobin is an erythrocyte protein belonging to the globin superfamily, and the α-, β-globin genes of jawed vertebrates encode tetrameric ((α₂β₂) hemoglobin, which contributes to aerobic metabolism by delivering oxygen from the respiratory exchange surfaces into cells. However, there are various gaps in knowledge regarding hemoglobin gene evolution, including patterns in cartilaginous fish and the roles of gene conversion in various taxa. Hence, we evaluated the evolutionary history of the vertebrate hemoglobin gene family by analyses of 97 species representing all classes of vertebrates. By genome-wide analyses, we extracted 879 hemoglobin sequences. Members of the hemoglobin gene family were conserved in birds and reptiles but variable in mammals, amphibians, and teleosts. Gene motifs, structures, and synteny were relatively well-conserved among vertebrates. Our results revealed that purifying selection contributed substantially to the evolution of all vertebrate hemoglobin genes, with mean d_N/d_S (ω) values ranging from 0.057 in teleosts to 0.359 in reptiles. In general, after the fish-specific genome duplication, the teleost hemoglobin genes showed variation in rates of evolution, and the β-globin genes showed relatively high ω values after a gene transposition event in amniotes. We also observed that the frequency of gene conversion was high in amniotes, with fewer hemoglobin genes and higher rates of evolution. Collectively, our findings provide detail insight into complex evolutionary processes shaping the vertebrate hemoglobin gene family, involving gene duplication, gene loss, purifying selection, and gene conversion.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Sinocyclocheilus longicornus (Cypriniformes, Cyprinidae), a new species of microphthalmic hypogean fish from Guizhou, Southwest China

Sinocyclocheilus longicornussp. nov. is described from the Pearl River basin in Hongguo Town, Panzhou City, Guizhou Province, Southwest China. Based on the presence of the long horn-like structure on the back of the head, Sinocyclocheilus longicornussp. nov. is assigned to the Sinocyclocheilus angularis species group. Sinocyclocheilus longicornussp. nov. is distinguished from its congeners by a combination of morphological characters: (1) presence of a single, relatively long horn-like structure on the back of the head; (2) pigmentation absent; (3) reduced eyes; (4) dorsal-fin rays, ii, 7; (5) pectoral-fin rays, i, 13; (6) anal-fin rays, iii, 5; (7) pelvic-fin rays, i, 7; (8) lateral line pores 38–49; (9) gill rakers well developed, nine on first gill arch; and (10) tip of adpressed pelvic fin not reaching anus.

Novel genome sequence of Chinese cavefish (Triplophysa rosa) reveals pervasive relaxation of natural selection in cavefish genomes

All cavefishes, living exclusively in caves across the globe, exhibit similar phenotypic traits, including the characteristic loss of eyes. To understand whether such phenotypic convergence shares similar genomic bases, here we investigated genome-wide evolutionary signatures of cavefish phenotypes by comparing whole-genome sequences of three pairs of cavefishes and their surface fish relatives. Notably, we newly sequenced and generated a whole-genome assembly of the Chinese cavefish Triplophysa rosa. Our comparative analyses revealed several shared features of cavefish genome evolution. Cavefishes had lower mutation rates than their surface fish relatives. In contrast, the ratio of nonsynonymous to synonymous substitutions (ω) was significantly elevated in cavefishes compared to in surface fishes, consistent with the relaxation of purifying selection. In addition, cavefish genomes had an increased mutational load, including mutations that alter protein hydrophobicity profiles, which were considered harmful. Interestingly, however, we found no overlap in positively selected genes among different cavefish lineages, indicating that the phenotypic convergence in cavefishes was not caused by positive selection of the same sets of genes. Analyses of previously identified candidate genes associated with cave phenotypes supported this conclusion. Genes belonging to the lipid metabolism functional ontology were under relaxed purifying selection in all cavefish genomes, which may be associated with the nutrient-poor habitat of cavefishes. Our work reveals previously uncharacterized patterns of cavefish genome evolution and provides comparative insights into the evolution of cave-associated phenotypic traits.

Sensory evolution in a cavefish radiation: patterns of neuromast distribution and associated behaviour in Sinocyclocheilus (Cypriniformes: Cyprinidae)

The genus Sinocyclocheilus, comprising a large radiation of freshwater cavefishes, are well known for their presence of regressive features (e.g. variable eye reduction). Fewer constructive features are known, such as the expansion of the lateral line system (LLS), which is involved in detecting water movements. The precise relationship between LLS expansion and cave adaptation is not well understood. Here, we examine morphology and LLS-mediated behaviour in Sinocyclocheilus species characterized by broad variation in eye size, habitat and geographical distribution. Using live-staining techniques and automated behavioural analyses, we examined 26 Sinocyclocheilus species and quantified neuromast organ number, density and asymmetry within a phylogenetic context. We then examined how these morphological features may relate to wall-following, an established cave-associated behaviour mediated by the lateral line. We show that most species demonstrated laterality (i.e. asymmetry) in neuromast organs on the head, often biased to the right. We also found that wall-following behaviour was distinctive, particularly among eyeless species. Patterns of variation in LLS appear to correlate with the degree of eye loss, as well as geographical distribution. This work reveals that constructive LLS evolution is convergent across distant cavefish taxa and may mediate asymmetric behavioural features that enable survival in stark subterranean microenvironments.

Differential transcriptomic responses to heat stress in surface and subterranean diving beetles

Subterranean habitats are generally very stable environments, and as such evolutionary transitions of organisms from surface to subterranean lifestyles may cause considerable shifts in physiology, particularly with respect to thermal tolerance. In this study we compared responses to heat shock at the molecular level in a geographically widespread, surface-dwelling water beetle to a congeneric subterranean species restricted to a single aquifer (Dytiscidae: Hydroporinae). The obligate subterranean beetle Paroster macrosturtensis is known to have a lower thermal tolerance compared to surface lineages (CT_max 38 °C cf. 42–46 °C), but the genetic basis of this physiological difference has not been characterized. We experimentally manipulated the thermal environment of 24 individuals to demonstrate that both species can mount a heat shock response at high temperatures (35 °C), as determined by comparative transcriptomics. However, genes involved in these responses differ between species and a far greater number were differentially expressed in the surface taxon, suggesting it can mount a more robust heat shock response; these data may underpin its higher thermal tolerance compared to subterranean relatives. In contrast, the subterranean species examined not only differentially expressed fewer genes in response to increasing temperatures, but also in the presence of the experimental setup employed here alone. Our results suggest P. macrosturtensis may be comparatively poorly equipped to respond to both thermally induced stress and environmental disturbances more broadly. The molecular findings presented here have conservation implications for P. macrosturtensis and contribute to a growing narrative concerning weakened thermal tolerances in obligate subterranean organisms at the molecular level.

Combined nanometric and phylogenetic analysis of unique endocytic compartments in Giardia lamblia sheds light on the evolution of endocytosis in Metamonada

Background

Giardia lamblia, a parasitic protist of the Metamonada supergroup, has evolved one of the most diverged endocytic compartment systems investigated so far. Peripheral endocytic compartments, currently known as peripheral vesicles or vacuoles (PVs), perform bulk uptake of fluid phase material which is then digested and sorted either to the cell cytosol or back to the extracellular space.

Results

Here, we present a quantitative morphological characterization of these organelles using volumetric electron microscopy and super-resolution microscopy (SRM). We defined a morphological classification for the heterogenous population of PVs and performed a comparative analysis of PVs and endosome-like organelles in representatives of phylogenetically related taxa, Spironucleus spp. and Tritrichomonas foetus. To investigate the as-yet insufficiently understood connection between PVs and clathrin assemblies in G. lamblia, we further performed an in-depth search for two key elements of the endocytic machinery, clathrin heavy chain (CHC) and clathrin light chain (CLC), across different lineages in Metamonada. Our data point to the loss of a bona fide CLC in the last Fornicata common ancestor (LFCA) with the emergence of a protein analogous to CLC (GlACLC) in the Giardia genus. Finally, the location of clathrin in the various compartments was quantified.

Conclusions

Taken together, this provides the first comprehensive nanometric view of Giardia’s endocytic system architecture and sheds light on the evolution of GlACLC analogues in the Fornicata supergroup and, specific to Giardia, as a possible adaptation to the formation and maintenance of stable clathrin assemblies at PVs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01402-3.

Genome resequencing clarifies phylogeny and reveals patterns of selection in the toxicogenomics model Pimephales promelas

Background

The fathead minnow (Pimephales promelas) is a model species for toxicological research. A high-quality genome reference sequence is available, and genomic methods are increasingly used in toxicological studies of the species. However, phylogenetic relationships within the genus remain incompletely known and little population-genomic data are available for fathead minnow despite the potential effects of genetic background on toxicological responses. On the other hand, a wealth of extant samples is stored in museum collections that in principle allow fine-scale analysis of contemporary and historical genetic variation.

Methods

Here we use short-read shotgun resequencing to investigate sequence variation among and within Pimephales species. At the genus level, our objectives were to resolve phylogenetic relationships and identify genes with signatures of positive diversifying selection. At the species level, our objective was to evaluate the utility of archived-sample resequencing for detecting selective sweeps within fathead minnow, applied to a population introduced to the San Juan River of the southwestern United States sometime prior to 1950.

Results

We recovered well-supported but discordant phylogenetic topologies for nuclear and mitochondrial sequences that we hypothesize arose from mitochondrial transfer among species. The nuclear tree supported bluntnose minnow (P. notatus) as sister to fathead minnow, with the slim minnow (P. tenellus) and bullhead minnow (P. vigilax) more closely related to each other. Using multiple methods, we identified 11 genes that have diversified under positive selection within the genus. Within the San Juan River population, we identified selective-sweep regions overlapping several sets of related genes, including both genes that encode the giant sarcomere protein titin and the two genes encoding the MTORC1 complex, a key metabolic regulator. We also observed elevated polymorphism and reduced differentation among populations (F_ST) in genomic regions containing certain immune-gene clusters, similar to what has been reported in other taxa. Collectively, our data clarify evolutionary relationships and selective pressures within the genus and establish museum archives as a fruitful resource for characterizing genomic variation. We anticipate that large-scale resequencing will enable the detection of genetic variants associated with environmental toxicants such as heavy metals, high salinity, estrogens, and agrichemicals, which could be exploited as efficient biomarkers of exposure in natural populations.

Swimming behavior and hydrodynamics of the Chinese cavefish Sinocyclocheilus rhinocerous and a possible role of its head horn structure

The blind troglobite cavefish Sinocyclocheilus rhinocerous lives in oligotrophic, phreatic subterranean waters and possesses a unique cranial morphology including a pronounced supra-occipital horn. We used a combined approach of laboratory observations and Computational Fluid Dynamics modeling to characterize the swimming behavior and other hydrodynamic aspects, i.e., drag coefficients and lateral line sensing distance of S. rhinocerous. Motion capture and tracking based on an Artificial Neural Network, complemented by a Particle Image Velocimetry system to map out water velocity fields, were utilized to analyze the motion of a live specimen in a laboratory aquarium. Computational Fluid Dynamics simulations on flow fields and pressure fields, based on digital models of S. rhinocerous, were also performed. These simulations were compared to analogous simulations employing models of the sympatric, large-eyed troglophile cavefish S. angustiporus. Features of the cavefish swimming behavior deduced from the both live-specimen experiments and simulations included average swimming velocities and three dimensional trajectories, estimates for drag coefficients and potential lateral line sensing distances, and mapping of the flow field around the fish. As expected, typical S. rhinocerous swimming speeds were relatively slow. The lateral line sensing distance was approximately 0.25 body lengths, which may explain the observation that specimen introduced to a new environment tend to swim parallel and near to the walls. Three-dimensional simulations demonstrate that just upstream from the region under the supra-occipital horn the equipotential of the water pressure and velocity fields are nearly vertical. Results support the hypothesis that the conspicuous cranial horn of S. rhinocerous may lead to greater stimulus of the lateral line compared to fish that do not possess such morphology.

Chromosome-level genome of Tibetan naked carp (Gymnocypris przewalskii) provides insights into Tibetan highland adaptation

Gymnocypris przewalskii, a cyprinid fish endemic to the Qinghai-Tibetan Plateau, has evolved unique morphological, physiological and genetic characteristics to adapt to the highland environment. Herein, we assembled a high-quality G. przewalskii tetraploid genome with a size of 2.03 Gb and scaffold N50 of 44.93 Mb, which was anchored onto 46 chromosomes. The comparative analysis suggested that gene families related to highland adaptation were significantly expanded in G. przewalskii. According to the G. przewalskii genome, we evaluated the phylogenetic relationship of 13 schizothoracine fishes, and inferred that the demographic history of G. przewalskii was strongly associated with geographic and eco-environmental alterations. We noticed that G. przewalskii experienced whole-genome duplication, and genes preserved post duplication were functionally associated with adaptation to high salinity and alkalinity. In conclusion, a chromosome-scale G. przewalskii genome provides an important genomic resource for teleost fish, and will particularly promote our understanding of the molecular evolution and speciation of fish in the highland environment.

The insulin gene as an energy homeostasis biomarker in Yangtze sturgeon (Acipenser dabryanus)

Insulin plays an important role in maintaining energy homeostasis and has the potential to be an indicator of energy homeostasis in the Yangtze sturgeon, Acipenser dabryanus. In this study, the Yangtze sturgeon insulin (Adinsulin) was cloned and characterized. To evaluate the possibility of insulin as an energy state assessment indicator, quantification real-time PCR (qRT-PCR) was used to evaluate expression changes in different tissues (the whole brain, esophagus, cardiac stomach, pyloric stomach, pyloric caeca, duodenum, valvula intestine, rectum, liver, pancreas, spleen, kidney, heart, muscle, gill and eye) from 6 fish (average weight 325.7 ± 22.3 g) and in three experiments including postprandial, fasting and re-feeding, and glucose tolerance treatment in which fish were divided into two groups including a group that administered a glucose solution (1 ul/g body weight) and another group that administered sterile water as control. In these three experiments, 6 fish were sampled, respectively, then been used to evaluate expression changes of insulin. All fish in feeding groups were fed in tanks (60.0 cm × 50.0 cm × 40.0 cm) with a commercial diet (crude protein ≥ 40%, crude fat ≥ 12%, coarse fiber ≤ 6%, crude ash ≤ 18%; TONGWEI CO., LTD, China) once a day at 16:00. The result showed that Adinsulin was highly expressed in the pancreas, which was the basis for the next experiment to use the pancreas as the test target. Adinsulin expression significantly increased 1 h after feeding and decreased rapidly after 3 h of feeding, but it was still significantly higher than that of the group without feeding (P < 0.01). Compared to the feeding group, the expression of Adinsulin was significantly reduced in the fasting group of 3 days (P < 0.01), 6 days (P < 0.01), 10 days (P < 0.05), 11 days (P < 0.05) and 13 days (P < 0.01) and was no significant difference in re-feeding for 1st day, 2nd day and 4th day, but there was difference between re-feeding group and fasting group. After glucose tolerance treatment, serum glucose levels increased significantly (P < 0.05), accompanied by a significant increase (P < 0.001) in insulin expression. This study result shows that insulin has the capacity to measure the energy homeostasis of Yangtze sturgeon. Further development of detection methods for sturgeon plasma or serum insulin will avoid slaughtering animals and is more practical in energy homeostasis assessment.

Evolutionary Genomics Reveals Multiple Functions of Arylalkylamine N-Acetyltransferase in Fish

As an important hormone, melatonin participates in endocrine regulation of diverse functions in vertebrates. Its biosynthesis is catalyzed by four cascaded enzymes, among them, arylalkylamine N-acetyltransferase (AANAT) is the most critical one. Although only single aanat gene has been identified in most groups of vertebrates, researchers including us have determined that fish have the most diverse of aanat genes (aanat1a, aanat1b, and aanat2), playing various potential roles such as seasonal migration, amphibious aerial vision, and cave or deep-sea adaptation. With the rapid development of genome and transcriptome sequencing, more and more putative sequences of fish aanat genes are going to be available. Related phylogeny and functional investigations will enrich our understanding of AANAT functions in various fish species.

Heterochrony and parallel evolution of echinoderm, hemichordate and cephalochordate internal bars

Deuterostomes comprise three phyla with radically different body plans. Phylogenetic bracketing of the living deuterostome clades suggests the latest common ancestor of echinoderms, hemichordates and chordates was a bilaterally symmetrical worm with pharyngeal openings, with these characters lost in echinoderms. Early fossil echinoderms with pharyngeal openings have been described, but their interpretation is highly controversial. Here, we critically evaluate the evidence for pharyngeal structures (gill bars) in the extinct stylophoran echinoderms Lagynocystis pyramidalis and Jaekelocarpus oklahomensis using virtual models based on high-resolution X-ray tomography scans of three-dimensionally preserved fossil specimens. Multivariate analyses of the size, spacing and arrangement of the internal bars in these fossils indicate they are substantially more similar to gill bars in modern enteropneust hemichordates and cephalochordates than to other internal bar-like structures in fossil blastozoan echinoderms. The close similarity between the internal bars of the stylophorans L. pyramidalis and J. oklahomensis and the gill bars of extant chordates and hemichordates is strong evidence for their homology. Differences between these internal bars and bar-like elements of the respiratory systems in blastozoans suggest these structures might have arisen through parallel evolution across deuterostomes, perhaps underpinned by a common developmental genetic mechanism.

Odor-Induced Vomiting Is Combinatorially Triggered by Palp Olfactory Receptor Neurons That Project to the Lobus Glomerulatus in Locust Brain

Although vomiting is commonly recognized as a protective reaction in response to toxic stimuli, the elaborate sensory processes and necessary molecular components are not fully clear, which is due to a lack of appropriate experimental animal models. Vomiting reflex to volatile chemicals renders locust one candidate for vomiting model. Here, we identified a panel of chemical cues that evoked evident vomiting in locust nymphs and demonstrated the selected combinatorial coding strategy that palps but not antennae olfactory receptor neurons (ORNs) employed. Specifically, knocking down individual palp odorant receptors (ORs) such as OR17, OR21, and OR22 attenuated the vomiting intensity evoked by E-2-hexenal and hexanal, while suppression of OR12 and OR22 augmented vomiting to E-2-hexenal and 2-hexanone, respectively. Furthermore, dual-RNAi treatment against OR17 or OR21 together with OR22 resulted in a much lower response intensity than that of individual OR suppression. Furthermore, OR12 was revealed in palp sensilla basiconica (pb) subtype 3 to tune the neuronal decaying activity to E-2-hexenal. Finally, anterograde labeling indicated that palp ORNs primarily projected into the lobus glomerulatus (LG), and the projection neurons (PNs) in the LG further projected into the accessary calyx (ACA). Together, the establishment of an olfaction-inducible vomiting model in locusts deepens the understanding of olfactory coding logics and provides an opportunity to clarify the neural basis underlying animal vomiting.

Chromosome-level genome assembly of grass carp (Ctenopharyngodon idella) provides insights into its genome evolution

Background

The grass carp has great economic value and occupies an important evolutionary position. Genomic information regarding this species could help better understand its rapid growth rate as well as its unique body plan and environmental adaptation.

Results

We assembled the chromosome-level grass carp genome using the PacBio sequencing and chromosome structure capture technique. The final genome assembly has a total length of 893.2 Mb with a contig N50 of 19.3 Mb and a scaffold N50 of 35.7 Mb. About 99.85% of the assembled contigs were anchored into 24 chromosomes. Based on the prediction, this genome contained 30,342 protein-coding genes and 43.26% repetitive sequences. Furthermore, we determined that the large genome size can be attributed to the DNA-mediated transposable elements which accounted for 58.9% of the repetitive sequences in grass carp. We identified that the grass carp has only 24 pairs of chromosomes due to the fusion of two ancestral chromosomes. Enrichment analyses of significantly expanded and positively selected genes reflected evolutionary adaptation of grass carp to the feeding habits. We also detected the loss of conserved non-coding regulatory elements associated with the development of the immune system, nervous system, and digestive system, which may be critical for grass carp herbivorous traits.

Conclusions

The high-quality reference genome reported here provides a valuable resource for the genetic improvement and molecular-guided breeding of the grass carp.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08503-x.

Adapting to Novel Environments Together: Evolutionary and Ecological Correlates of the Bacterial Microbiome of the World's Largest Cavefish Diversification (Cyprinidae, Sinocyclocheilus)

The symbiosis between a host and its microbiome is essential for host fitness, and this association is a consequence of the host’s physiology and habitat. Sinocyclocheilus, the largest cavefish diversification of the world, an emerging multi-species model system for evolutionary novelty, provides an excellent opportunity for examining correlates of host evolutionary history, habitat, and gut-microbial community diversity. From the diversification-scale patterns of habitat occupation, major phylogenetic clades (A–D), geographic distribution, and knowledge from captive-maintained Sinocyclocheilus populations, we hypothesize habitat to be the major determinant of microbiome diversity, with phylogeny playing a lesser role. For this, we subject environmental water samples and fecal samples (representative of gut-microbiome) from 24 Sinocyclocheilus species, both from the wild and after being in captivity for 6 months, to bacterial 16S rRNA gene profiling using Illumina sequencing. We see significant differences in the gut microbiota structure of Sinocyclocheilus, reflective of the three habitat types; gut microbiomes too, were influenced by host-related factors. There is no significant association between the gut microbiomes and host phylogeny. However, there is some microbiome related structure at the clade level, with the most geographically distant clades (A and D) being the most distinct, and the two overlapping clades (B and C) showing similarities. Microbes inhabiting water were not a cause for significant differences in fish-gut microbiota, but water quality parameters were. Transferring from wild to captivity, the fish microbiomes changed significantly and became homogenized, signifying plastic changes and highlighting the importance of environmental factors (habitat) in microbiome community assembly. The core microbiome of this group, at higher taxonomic scale, resembled that of other teleost fishes. Our results suggest that divergent natural environments giving rise to evolutionary novelties underlying host adaptations, also includes the microbiome of these fishes.

A near-chromosome level genome assembly of the European hoverfly, Sphaerophoria rueppellii (Diptera: Syrphidae), provides comparative insights into insecticide resistance-related gene family evolution

BACKGROUND

Sphaerophoria rueppellii, a European species of hoverfly, is a highly effective beneficial predator of hemipteran crop pests including aphids, thrips and coleopteran/lepidopteran larvae in integrated pest management (IPM) programmes. It is also a key pollinator of a wide variety of important agricultural crops. No genomic information is currently available for S. rueppellii. Without genomic information for such beneficial predator species, we are unable to perform comparative analyses of insecticide target-sites and genes encoding metabolic enzymes potentially responsible for insecticide resistance, between crop pests and their predators. These metabolic mechanisms include several gene families - cytochrome P450 monooxygenases (P450s), ATP binding cassette transporters (ABCs), glutathione-S-transferases (GSTs), UDP-glycosyltransferases (UGTs) and carboxyl/choline esterases (CCEs).

METHODS AND FINDINGS

In this study, a high-quality near-chromosome level de novo genome assembly (as well as a mitochondrial genome assembly) for S. rueppellii has been generated using a hybrid approach with PacBio long-read and Illumina short-read data, followed by super scaffolding using Hi-C data. The final assembly achieved a scaffold N50 of 87Mb, a total genome size of 537.6Mb and a level of completeness of 96% using a set of 1,658 core insect genes present as full-length genes. The assembly was annotated with 14,249 protein-coding genes. Comparative analysis revealed gene expansions of CYP6Zx P450s, epsilon-class GSTs, dietary CCEs and multiple UGT families (UGT37/302/308/430/431). Conversely, ABCs, delta-class GSTs and non-CYP6Zx P450s showed limited expansion. Differences were seen in the distributions of resistance-associated gene families across subfamilies between S. rueppellii and some hemipteran crop pests. Additionally, S. rueppellii had larger numbers of detoxification genes than other pollinator species.

CONCLUSION AND SIGNIFICANCE

This assembly is the first published genome for a predatory member of the Syrphidae family and will serve as a useful resource for further research into selectivity and potential tolerance of insecticides by beneficial predators. Furthermore, the expansion of some gene families often linked to insecticide resistance and selectivity may be an indicator of the capacity of this predator to detoxify IPM selective insecticides. These findings could be exploited by targeted insecticide screens and functional studies to increase effectiveness of IPM strategies, which aim to increase crop yields by sustainably and effectively controlling pests without impacting beneficial predator populations.

The first de novo genome assembly and sex marker identification of Pluang Chomphu fish (Tor tambra) from Southern Thailand

The Tor genus belongs to the group of cyprinid fish commonly known as mahseer. Although Tor species are rapidly declining in the wild, and some face extinction, ambiguities in species identification hinder their collection and conservation. We conducted a genome survey of male and female Tor tambra collected in Thailand. The genome sizes of the male and female fish were approximately 1,671 and 1,645 Mb, respectively, with repeat contents of approximately 33%. The heterozygosity ratios of the male and female fish, which were 0.34% and 0.39%, respectively, suggested that the sex of T. tambra is determined by the ZW system. A sex marker was identified in silico and confirmed by PCR amplification. The result indicated that T. tambra has a ZZ/ZW sex determination system. Subsequently, comparative genomic and phylogenetic analyses of T. tambra and other fish in the Cyprinidae family were performed to explore the genetic diversity and evolution of the species. We also assembled the complete mitochondrial genome sequences of the T. tambra collected in Thailand. A phylogenetic tree of different Tor species, constructed based on mitochondrial genome sequences, indicated that T. tambra was closely related to T. tambroides. We believe this is the first genome survey of a species from the Tor genus or Mahseer group. Our results may help identify Tor species, providing a reference for genetic studies of the Tor genus and other mahseer fish.

Hybridization underlies localized trait evolution in cavefish

Introgressive hybridization may play an integral role in local adaptation and speciation (Taylor and Larson, 2019). In the Mexican tetra Astyanax mexicanus, cave populations have repeatedly evolved traits including eye loss, sleep loss, and albinism. Of the 30 caves inhabited by A. mexicanus, Chica cave is unique because it contains multiple pools inhabited by putative hybrids between surface and cave populations (Mitchell et al., 1977), providing an opportunity to investigate the impact of hybridization on complex trait evolution. We show that hybridization between cave and surface populations may contribute to localized variation in traits associated with cave evolution, including pigmentation, eye development, and sleep. We also uncover an example of convergent evolution in a circadian clock gene in multiple cavefish lineages and burrowing mammals, suggesting a shared genetic mechanism underlying circadian disruption in subterranean vertebrates. Our results provide insight into the role of hybridization in facilitating phenotypic evolution.

•

Hybridization leads to a localized difference in sleep duration within a single cave

•

Genomic analysis identifies coding differences in Cry1A across cave pools

•

Changes in Cry1A appear to be conserved in cavefish and burrowing mammals

Fecal and soil microbiota composition of gardening and non-gardening families

Historically, humans have interacted with soils, which contain a rich source of microorganisms. Fruit and vegetable gardening is the primary interaction humans have with soil today. Animal research reveals that soil microorganisms can be transferred to the rodent intestine. However, studies on fecal and soil microbial changes associated with gardening in humans are lacking. The current case-controlled cohort study aimed to characterize the fecal and soil microbiota of gardening families (n = 10) and non-gardening (control) families (n = 9). Families included two adults and one child (5-18 years) for a total of 56 participants. All participants provided a fecal sample, soil sample, and diet history questionnaires before the gardening season (April) and during the peak of the gardening season (August). Healthy Eating Index (HEI-2015) scores and nutrient analysis were performed. Fecal and soil DNA were extracted and amplified. Sequence data were then processed and analyzed. Peak season gardening families tended to have greater fecal operational features, a greater Faith's Phylogenetic Diversity score, greater fiber intake, and higher abundances of fiber fermenting bacteria than peak control families. Soil endemic microbes were also shared with gardening participant's fecal samples. This study revealed that the fecal microbiota of gardening families differs from non-gardening families, and that there are detectable changes in the fecal microbial community of gardeners and their family members over the course of the gardening season. Additional research is necessary to determine if changes induced by gardening on the gut microbiota contribute to human health.

Dominance of Asia II 1 species of Bemisia tabaci in Pakistan and beyond

Globally, Whitefly (Bemisia tabaci) is one of the most important insect pests of crops that causes huge economical losses. The current study was designed to exclusively screen the B. tabaci species in the cotton field of Pakistan during 2017-2020 and have to conduct comparative analysis of B. tabaci species in Asia where Asia II 1 has been reported. A total of 5142 B. tabaci sequences of mitochondrial cytochrome oxidase 1 (mtCO1) from Asian countries were analyzed to determine the species and their distribution in the region. Our analysis over time and space showed that Asia II 1 has gradually dominated over Asia 1 in Punjab Province and over both Asia 1 and MEAM1 in Sindh Province. Asia has been divided into three regions i.e., South Asia (2524 sequences), Southeast Asia (757 sequences) and East Asia (1569 sequences) and dominance of different species of B. tabaci has been determined by calculating the relative percentage of each species. Interestingly, Asia II 1 has been found dominant in the neighboring region (northern zone) of India and also being dominant in its central zone. The dominance of Asia II 1 in Pakistan and northern India explains whitefly epidemic being reported in recent years.

Evolving in the darkness: Phylogenomics of Sinocyclocheilus cavefishes highlights recent diversification and cryptic diversity

Troglomorphism-any morphological adaptation enabling life to the constant darkness of caves, such as loss of pigment, reduced eyesight or blindness, over-developed tactile and olfactory organs-has long intrigued biologists. However, inferring the proximate and ultimate mechanisms driving the evolution of troglomorphism (stygomorphism) in freshwater fish requires a sound understanding of the evolutionary relationships between surface and stygomorphic lineages. We use Restriction Site Associated DNA Sequencing (RADseq) to better understand the evolution of the Sinocyclocheilus fishes of China. With a remarkable array of derived stygomorphic traits, they comprise the largest cavefish diversification in the world, emerging as a multi-species model system to study evolutionary novelty. We sequenced a total of 120 individuals throughout the Sinocyclocheilus distribution. The data comprised a total of 646,497 bp per individual, including 4378 loci and 67,983 SNPs shared across a minimum of 114 individuals at a given locus. Phylogenetic analyses using either the concatenated RAD loci (RAxML) or the SNPs under a coalescent model (SVDquartets, SNAPP) showed a high degree of congruence with similar topologies and high node support (>95 for most nodes in the phylogeny). The major clades recovered conform to a pattern previously established using Sanger-based mt-DNA sequences, with a few notable exceptions. We now recognize six major clades in this group, elevating the blind cavefish S. tianlinensis and the micro-eyed S. microphthalmus as two new distinct clades due to their deep divergence from other clades. PCA plots of the SNP data also support the recognition of six major clusters of species congruent with the identified clades in ordination space. A Bayes factor delimitation (BFD) analysis showed support for 21 species, recognizing 19 previously described species and two putative new cryptic ones. Two species whose identities were previously disputed, S. furcodorsalis and S. tianeensis, are supported here as distinct species. In addition, our multi-species calibrated tree in SNAPP suggests that the genus Sinocyclocheilus originated around 10.16 Mya, with most speciation events occurring in the last 2 Mya, likely favored by the uplift of the Qinghai-Tibetan Plateau and cave occupation induced by climate-driven aridification during this period. These results provide a firm basis for future comparative studies on the evolution of Sinocyclocheilus and its adaptations to cave life.

A scaffold-level genome assembly of a minute pirate bug, Orius laevigatus (Hemiptera: Anthocoridae), and a comparative analysis of insecticide resistance-related gene families with hemipteran crop pests

BACKGROUND

Orius laevigatus, a minute pirate bug, is a highly effective beneficial predator of crop pests including aphids, spider mites and thrips in integrated pest management (IPM) programmes. No genomic information is currently available for O. laevigatus, as is the case for the majority of beneficial predators which feed on crop pests. In contrast, genomic information for crop pests is far more readily available. The lack of publicly available genomes for beneficial predators to date has limited our ability to perform comparative analyses of genes encoding potential insecticide resistance mechanisms between crop pests and their predators. These mechanisms include several gene/protein families including cytochrome P450s (P450s), ATP binding cassette transporters (ABCs), glutathione S-transferases (GSTs), UDP-glucosyltransferases (UGTs) and carboxyl/cholinesterases (CCEs).

METHODS AND FINDINGS

In this study, a high-quality scaffold level de novo genome assembly for O. laevigatus has been generated using a hybrid approach with PacBio long-read and Illumina short-read data. The final assembly achieved a scaffold N50 of 125,649 bp and a total genome size of 150.98 Mb. The genome assembly achieved a level of completeness of 93.6% using a set of 1658 core insect genes present as full-length genes. Genome annotation identified 15,102 protein-coding genes - 87% of which were assigned a putative function. Comparative analyses revealed gene expansions of sigma class GSTs and CYP3 P450s. Conversely the UGT gene family showed limited expansion. Differences were seen in the distributions of resistance-associated gene families at the subfamily level between O. laevigatus and some of its targeted crop pests. A target site mutation in ryanodine receptors (I4790M, PxRyR) which has strong links to diamide resistance in crop pests and had previously only been identified in lepidopteran species was found to also be present in hemipteran species, including O. laevigatus.

CONCLUSION AND SIGNIFICANCE

This assembly is the first published genome for the Anthocoridae family and will serve as a useful resource for further research into target-site selectivity issues and potential resistance mechanisms in beneficial predators. Furthermore, the expansion of gene families often linked to insecticide resistance may be an indicator of the capacity of this predator to detoxify selective insecticides. These findings could be exploited by targeted pesticide screens and functional studies to increase effectiveness of IPM strategies, which aim to increase crop yields by sustainably, environmentally-friendly and effectively control pests without impacting beneficial predator populations.

Albinism in the largest extant amphibian: A metabolic, endocrine, or immune problem?

BACKGROUND

Pigment regression is an intriguing phenomenon that can be caused by disorders in melanin metabolism or endocrine regulation, or by autoimmune disorders. Albino animals serve as excellent models for the study of the genetic determination of morphology, particularly the evolution of and molecular mechanisms underlying chromatophore-related diseases in animals and humans.

MATERIAL AND METHODS

The artificial culture of Andrias davidianus, the largest extant amphibian, is flourishing in China due to the great ecological and economic value of this animal. Approximately 0.1% of individuals express an albino phenotype accompanied by delayed somatic growth and mortality at early developmental stages. In this study, brain and skin transcriptomics were conducted to study the underlying molecular basis of the phenotype.

RESULTS

The results indicated decreased transcription of genes of melanin synthesis. Interestingly, MHC I isotypes and immune-related pathways accounted for the primary transcriptional differences between groups, suggesting that the albino phenotype represents a systematic immune problem to a far greater extent than a pigmentation defect. Albino individuals exhibited shifted transcription of MHC I isotypes, and the albino-specific isotype was characterized by increased charges and decreased space in the antigen- binding pocket, implying a drastic change in antigen specificity and a potential risk of autoimmune disorders.

CONCLUSION

These results suggest an association between the albino phenotype and MHC I variants in A. davidianus, which could serve as a convenient model for vitiligo or other autoimmune diseases.

Molecular phylogeny and historical biogeography of the cave fish genus Sinocyclocheilus (Cypriniformes: Cyprinidae) in southwest China

Modern accumulations of genetic data offer unprecedented opportunities for understanding the systematic classification and origins of specific groups of organisms. The genus Sinocyclocheilus is among the most cave abundant genera in Cyprinidae, with 76 recognized species, belonging to four species groups. Recent phylogenetic studies have shown that the classification of species groups within the genus Sinocyclocheilus remains controversial. In this study, we constructed a sequence supermatrix of 26 species from four species groups of the genus Sinocyclocheilus using the mitochondrial genome to reveal phylogenetic relationships, historical biogeography and patterns of species diversification in the genus Sinocyclocheilus. Phylogenetic analysis strongly supports the monophyletic groups of the three species groups (S. jii, S. cyphotergous, and S. tingi groups) except the S. angularis group. Phylogenetic analysis showed that S. anshuiensis and S. microphthalmus, which were recognized as numbers of S. angularis group, formed a strongly supported independent clade. Therefore, we propose a new species group, the S. microphthalmus group, which contains S. anshuiensis and S. microphthalmus. Biogeographic reconstruction suggests that the living Sinocyclocheilus may have originated in north-central Guangxi at the late Eocene and dispersed outward after a vicariance at 32.31 Million years ago (Ma). Early diversification is focused on the late Oligocene (ca. 25 Ma), which is related to the second uplift of the Qinghai-Tibetan Plateau and the lateral extrusion of the Indochina at the Oligocene/Miocene boundary. Our results suggest that two uplifts of the Qinghai-Tibetan Plateau and climate change in the Miocene may have influenced the diversification of the Sinocyclocheilus lineage. This article is protected by copyright. All rights reserved.

Integrative mRNA-miRNA interaction analysis reveals the molecular mechanism of skin color variation between wild-type and yellow mutant rainbow trout (Oncorhynchus mykiss)

Rainbow trout (Oncorhynchus mykiss) is an important economic fish in China. Skin color affects the economic value of trout. However, the molecular mechanism of the skin color variation between wild-type (WR) and yellow mutant rainbow trout (YR) is unclear. We sequenced mRNAs and miRNAs of dorsal skin to identify key color variation-associated mRNAs and miRNAs between WR and YR. Overall, 2060 out of 3625 differentially expressed genes were upregulated in YR, and 196 out of 275 differentially expressed miRNAs were downregulated in WR. We identified three key YR-upregulated genes related to the formation of xanthophores (GCH1, SLC2A11, and SOX10). Interestingly, several genes related to melanogenesis (TYR, TYRP1, TYRP2, MC1R, MITF, PMEL, SLC45A2, and OCA2) were downregulated in WR. Integrated analysis identified five miRNAs that target at least two skin color-related genes (miR-495-y, miR-543-y, miR-665-z, miR-433-y, and miR-382-x). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses of target genes identified noncoding RNA metabolic process as the most significantly enriched GO term, and several metabolic pathways associated with skin color were enriched significantly, such as tyrosine metabolism, histidine metabolism, and vitamin B6 metabolism. Quantitative real-time PCR of selected mRNAs and miRNAs validated the reliability of the integrated analysis. This study provides in-depth insights into the molecular mechanism of skin color variation between WR and YR, which will accelerate the genetic selection and breeding of rainbow trout with consumer-favored traits.