Genomic insights into the adaptation and evolution of the nautilus,an ancient but evolving “living fossil”

The nautilus, commonly known as a “living fossil,” is endangered and may be at risk of extinction. The lack of genomic information hinders a thorough understanding of its biology and evolution, which can shed light on the conservation of this endangered species. Here, we report the first high-quality chromosome-level genome assembly of Nautilus pompilius. The assembled genome size comprised 785.15 Mb. Comparative genomic analyses indicated that transposable elements (TEs) and large-scale genome reorganizations may have driven lineage-specific evolution in the cephalopods. Remarkably, evolving conserved genes and recent TE insertion activities were identified in N. pompilius, and we speculate that these findings reflect the strong adaptability and long-term survival of the nautilus. We also identified gene families that are potentially responsible for specific adaptation and evolution events. Our study provides unprecedented insights into the specialized biology and evolution of N. pompilius, and the results serve as an important resource for future conservation genomics of the nautilus and closely related species.     

Chromosome-level genome assembly of Scapharca kagoshimensis reveals the expanded molecular basis of heme biosynthesis in ark shells

Ark shells are commercially important clam species that inhabit in muddy sediments of shallow coasts in East Asia. For a long time, the lack of genome resources has hindered scientific research of ark shells. Here, we report a high-quality chromosome-level genome assembly of Scapharca kagoshimensis, with an aim to unravel the molecular basis of heme biosynthesis, and develop genomic resources for genetic breeding and population genetics in ark shells. Nineteen scaffolds corresponding to 19 chromosomes were constructed from 938 contigs (contig N50 = 2.01 Mb) to produce a final high-quality assembly with a total length of 1.11 Gb and scaffold N50 around 60.64 Mb. The genome assembly represents 93.4% completeness via matching 303 eukaryota core conserved genes. A total of 24,908 protein-coding genes were predicted and 24,551 genes (98.56%) of which were functionally annotated. The enrichment analyses suggested that genes in heme biosynthesis pathways were expanded and positive selection of the haemoglobin genes was also found in the genome of S. kagoshimensis, which gives important insights into the molecular mechanisms and evolution of the heme biosynthesis in mollusca. The valuable genome assembly of S. kagoshimensis would provide a solid foundation for investigating the molecular mechanisms that underlie the diverse biological functions and evolutionary adaptations of S. kagoshimensis.     

Chromosome-level genome assembly of an agricultural pest,the rice leaffolder Cnaphalocrocis exigua (Crambidae,Lepidoptera)

The rice leaffolder Cnaphalocrocis exigua (Crambidae, Lepidoptera) is an important agricultural pest that damages rice crops and other members of related grass families. C. exigua exhibits a very similar morphological phenotype and feeding behaviour to C. medinalis, another species of rice leaffolder whose genome was recently reported. However, genomic information for C. exigua remains extremely limited. Here, we used a hybrid strategy combining different sequencing technologies, including Illumina, PacBio, 10× Genomics, and Hi – C scaffolding, to generate a high-quality chromosome-level genome assembly of C. exigua. We initially obtained a 798.8 Mb assembly with a contig N50 size of 2.9 Mb, and the N50 size was subsequently increased to 25.7 Mb using Hi – C technology to anchor 1413 scaffolds to 32 chromosomes. We detected a total of 97.7% Benchmarking Universal Single-Copy Orthologues (BUSCO) in the genome assembly, which was comprised of ~52% repetitive sequence and annotated 14,922 protein-coding genes. Of note, the Z and W sex chromosomes were assembled and identified. A comparative genomic analysis demonstrated that despite the high synteny observed between the two rice leaffolders, the species have distinct genomic features associated with expansion and contraction of gene families and selection pressure. In summary, our chromosome-level genome assembly and comparative genomic analysis of C. exigua provide novel insights into the evolution and ecology of this rice insect pests and offer useful information for pest control.     

Chemical Compositions and Antioxidant Activities of Essential Oils Extracted from the Petals of Three Wild Tree Peony Species and Eleven Cultivars

The aim of this study was to investigate the essential oil (EO) compositions and antioxidant activities from petals of three wild tree peony species (Paeonia delavayi, P. lutea, and P. rockii) and eleven P. suffruticosa cultivars from different cultivar groups. The EOs yields varied from 0.63% to 1.25% (v/v) among samples when using supercritical CO₂ extraction. One hundred and sixty‐three components were detected by GC/MS; and among them, linalool oxide, (Z)‐5‐dodecen‐1‐yl acetate, nonadecane, (Z)‐5‐nonadecene, heneicosane, phytol, and linoleic acid ethyl ester were dominant. According to hierarchical cluster analysis, principal component analysis and correspondence analysis, P. lutea, P. delavayi, and ‘High Noon’ were clustered in a group described as having a refreshing herbal‐like note due to high rates of phytol and linalool oxide. Notably, P. lutea and P. delavayi also had strong DPPH and ABTS radical scavenging activities. These results suggest that P. lutea and P. delavayi are the most promising candidates as useful sources of fragrances and natural antioxidants.     

Growth characteristics of Uranoscopus scaber Linnaeus 1758, inhabiting the Iskenderun Bay (Northeastern Mediterranean)

This study aimed to determine the age and some growth characteristics of Atlantic stargazer (Uranoscopus scaber) from Iskenderun Bay (Northeastern Mediterranean). For this purpose, a total of 150 Atlantic stargazer ranging in size from 9.1 to 28.0 cm in total length (weight: 11.7–345.7 g) were collected as by-catch from a commercial trawl fishing boat at a depth of 80–100 m between May 2015 and January 2016. The bottom trawl gear used was equipped with a 44 mm stretched mesh size net at the cod-end. The percentage of females and males were 46.7% and 53.3% respectively. The total length–weight relationships equation with coefficient of determination (R²) were found as W = 0.011*TL^3.131, R² = .9728, for all individuals, W = 0.015*TL^3.021, R² = .9512 for females and W = 0.0102*TL^3.136, R² = .9553 for males. By using the von Bertalanffy equation, the growth parameters of Atlantic stargazer were estimated to be L_∞ = 42.35 cm, k = 0.098, t₀ = −1.8474 for all individuals; L_∞ = 36.92 cm, k = 0.138, t₀ = −1.2693 for females and L_∞ = 38.77 cm, k = 0.1, t₀ = −2.334 for males. In this study, age reading was done by two independent readers and index of average percentage error (IAPE) was calculated as 6.1%. The highest condition factor value calculated as 1.81 in the age group 6 and the lowest condition factor value was calculated as 1.48 in the age group 1.     

Melanocortin 1 receptor (MC1R) polymorphisms’ influence on size and dermoscopic features of nevi

The melanocortin 1 receptor (MC1R) is a highly polymorphic gene. The loss‐of‐function MC1R variants (“R”) have been strongly associated with red hair color phenotype and an increased melanoma risk. We sequenced the MC1R gene in 175 healthy individuals to assess the influence of MC1R on nevus phenotype. We identified that MC1R variant carriers had larger nevi both on the back [p‐value = .016, adjusted for multiple parameters (adj. p‐value)] and on the upper limbs (adj. p‐value = .007). Specifically, we identified a positive association between the “R” MC1R variants and visible vessels in nevi [p‐value = .033, corrected using the FDR method for multiple comparisons (corrected p‐value)], dots and globules in nevi (corrected p‐value = .033), nevi with eccentric hyperpigmentation (corrected p‐value = .033), a high degree of freckling (adj. p‐value = .019), and an associative trend with presence of blue nevi (corrected p‐value = .120). In conclusion, the MC1R gene appears to influence the nevus phenotype.     

Genetic variation in nuclear and mitochondrial markers supports a large sex difference in lifetime reproductive skew in a lekking species

Sex differences in skews of vertebrate lifetime reproductive success are difficult to measure directly. Evolutionary histories of differential skew should be detectable in the genome. For example, male‐biased skew should reduce variation in the biparentally inherited genome relative to the maternally inherited genome. We tested this approach in lek‐breeding ruff (Class Aves, Philomachus pugnax) by comparing genetic variation of nuclear microsatellites (θ_n; biparental) versus mitochondrial D‐loop sequences (θ_m; maternal), and conversion to comparable nuclear (N_e) and female (N_ef) effective population size using published ranges of mutation rates for each marker (μ). We provide a Bayesian method to calculate N_e (θ_n = 4N_eμ_n) and N_ef (θ_m = 2N_efμ_m) using 95% credible intervals (CI) of θ_n and θ_m as informative priors, and accounting for uncertainty in μ. In 96 male ruffs from one population, N_e was 97% (79–100%) lower than expected under random mating in an ideal population, where N_e:N_ef = 2. This substantially lower autosomal variation represents the first genomic support of strong male reproductive skew in a lekking species.     

Sequencing of Camelina neglecta,a diploid progenitor of the hexaploid oilseed Camelina sativa

Camelina neglecta is a diploid species from the genus Camelina, which includes the versatile oilseed Camelina sativa. These species are closely related to Arabidopsis thaliana and the economically important Brassica crop species, making this genus a useful platform to dissect traits of agronomic importance while providing a tool to study the evolution of polyploids. A highly contiguous chromosome-level genome sequence of C. neglecta with an N50 size of 29.1 Mb was generated utilizing Pacific Biosciences (PacBio, Menlo Park, CA) long-read sequencing followed by chromosome conformation phasing. Comparison of the genome with that of C. sativa shows remarkable coincidence with subgenome 1 of the hexaploid, with only one major chromosomal rearrangement separating the two. Synonymous substitution rate analysis of the predicted 34 061 genes suggested subgenome 1 of C. sativa directly descended from C. neglecta around 1.2 mya. Higher functional divergence of genes in the hexaploid as evidenced by the greater number of unique orthogroups, and differential composition of resistant gene analogs, might suggest an immediate adaptation strategy after genome merger. The absence of genome bias in gene fractionation among the subgenomes of C. sativa in comparison with C. neglecta, and the complete lack of fractionation of meiosis-specific genes attests to the neopolyploid status of C. sativa. The assembled genome will provide a tool to further study genome evolution processes in the Camelina genus and potentially allow for the identification and exploitation of novel variation for Camelina crop improvement.     

The effects of environmental variables and human disturbance on woody species richness and diversity in a bamboo–deciduous forest in northeastern Thailand

Variations in species richness and diversity at a local scale are affected by a number of complex and interacting variables, including both natural environmental factors and human-made changes to the local environment. Here we identified the most important determinants of woody species richness and diversity at different growth stages (i.e. adult, sapling and seedling) in a bamboo–deciduous forest in northeast Thailand. A total of 20 environmental and human disturbance variables were used to determine the variation in species richness and diversity. In total, we identified 125 adult, 111 sapling (within fifty 20 × 20-m plots) and 89 seedling species (within one hundred and twenty 1 × 1-m subplots). Overall results from stepwise multiple regression analyses showed that environmental variables were by far the most important in explaining the variation in species richness and diversity. Forest structure (i.e. number of bamboo clumps and canopy cover) was important in determining the adult species richness and diversity (R ² = 0.48, 0.30, respectively), while topography (i.e. elevation) and human disturbance (i.e. number of tree stumps) were important in determining the sapling species richness and diversity (R ² = 0.55, 0.39, respectively). Seedling species richness and diversity were negatively related to soil phosphorus. Based on our results, we suggest that the presence of bamboos should be incorporated in management strategies for maintaining woody species richness and diversity in these forest ecosystems. Specifically, if bamboos cover the forest floor at high densities, it may be necessary to actively control these species for successful tree establishment.     

SPIKEPIPE: A metagenomic pipeline for the accurate quantification of eukaryotic species occurrences and intraspecific abundance change using DNA barcodes or mitogenomes

The accurate quantification of eukaryotic species abundances from bulk samples remains a key challenge for community ecology and environmental biomonitoring. We resolve this challenge by combining shotgun sequencing, mapping to reference DNA barcodes or to mitogenomes, and three correction factors: (a) a percent‐coverage threshold to filter out false positives, (b) an internal‐standard DNA spike‐in to correct for stochasticity during sequencing, and (c) technical replicates to correct for stochasticity across sequencing runs. The SPIKEPIPE pipeline achieves a strikingly high accuracy of intraspecific abundance estimates (in terms of DNA mass) from samples of known composition (mapping to barcodes R² = .93, mitogenomes R² = .95) and a high repeatability across environmental‐sample replicates (barcodes R² = .94, mitogenomes R² = .93). As proof of concept, we sequence arthropod samples from the High Arctic, systematically collected over 17 years, detecting changes in species richness, species‐specific abundances, and phenology. SPIKEPIPE provides cost‐efficient and reliable quantification of eukaryotic communities.     

Population genetic diversity and structure of two rare vernal pool grasses in central California

Vernal pool ecosystems are declining throughout California, with only 10% of historic habitat remaining. This has endangered many specialist endemic plant species, leaving extant populations fragmented, isolated, and threatened or endangered. Recovery plans for the increasing number of endangered vernal pool species require information on their genetic and ecological status to guide conservation and restoration efforts. Federally threatened Neostapfia colusana (Colusa grass) and federally endangered Tuctoria greenei (Greene’s tuctoria) are two endemic vernal pool grasses of high conservation concern in central California. Remaining populations are highly fragmented due to range-wide habitat destruction. Using five polymorphic microsatellite markers for each species, we performed genetic surveys of 240 individuals from eight vernal pools for N. colusana, and 317 individuals from 13 vernal pools for T. greenei. We detected high within-population genetic diversity for both species, with average allelic diversities of 24 alleles/locus (mean Hobs = 0.68, mean Hexp = 0.71) for N. colusana, and 19 alleles/locus (mean Hobs = 0.77, and mean Hexp = 0.79) for T. greenei. Bayesian clustering and AMOVA indicated two genetically distinct population groups for N. colusana (Fst = 0.268, P < 0.0001), and three for T. greenei (Fst = 0.11, P < 0.0001). We found very slight temporal genetic structure at one N. colusana (Fst = 0.013, P < 0.05) and two T. greenei (Fst = 0.015, Fst = 0.018, P < 0.05) pools. These estimates of population genetic diversity and structure are critical measures for both species that will help inform recovery management actions.     

Comparative Genome Analysis Reveals Cyanidiococcus gen. nov., A New Extremophilic Red Algal Genus Sister to Cyanidioschyzon (Cyanidioschyzonaceae,Rhodophyta)

The taxonomic placement of strains belonging to the extremophilic red alga Galdieria maxima has been controversial due to the inconsistent phylogenetic position inferred from molecular phylogenetic analyses. Galdieria maxima nom. inval. was classified in this genus based on morphology and molecular data in the early work, but some subsequent molecular phylogenetic analyses have inferred strains of G. maxima to be closely related to the genus Cyanidioschyzon. To address this controversy, an isolated strain identified as G. maxima using the rbcL gene sequence as the genetic barcode was examined using a comprehensive analysis across morphological, physiological, and genomic traits. Herein are reported the chloroplast-, mitochondrion-, and chromosome-level nuclear genome assemblies. Comparative analysis of orthologous gene clusters and genome arrangements suggested that the genome structure of this strain was more similar to that of the generitype of Cyanidioschyzon, C. merolae than to the generitype of Galdieria, G. sulphuraria. While the ability to uptake various forms of organic carbon for growth is an important physiological trait of Galdieria, this strain was identified as an ecologically obligate photoautotroph (i.e., the inability to utilize the natural concentrations of organic carbons) and lacked various gene models predicted as sugar transporters. Based on the genomic, morphological, and physiological traits, we propose this strain to be a new genus and species, Cyanidiococcus yangmingshanensis. Re-evaluation of the 18S rRNA and rbcL gene sequences of the authentic strain of G. maxima, IPPAS-P507, with those of C. yangmingshanensis suggests that the rbcL sequences of “G. maxima” deposited in GenBank correspond to misidentified isolates.     

Genomic analyses of a “living fossil”: The endangered dove‐tree

Davidia involucrata Baill, also known as the dove‐tree, is a living fossil and an endangered species currently restricted to the mountains of southwestern and central China. It has a beautiful and innovative trait of high horticultural value: two white bracts covering the flower caputila. Here, we report on the chromosome‐scale genome of this species using single‐molecule real‐time long reads and chromosome conformation capture (Hi‐C) techniques. This species has a larger genome size of 1,169 Mb and contains relatively more genes (42,554) than the closely related species Camptotheca acuminata (397 Mb and 31,825 genes). Both species shared one recent whole genome duplication before their divergence. The expansion of the repetitive elements after their divergence contributed greatly to the increase in the genome size of the dove‐tree. Photosynthesis‐related genes were almost absent or showed reduced expression in the bracts of the dove‐tree, while defence‐ and chemical‐related genes increased greatly, highlighting the important roles of the bracts in protecting flowers and attracting pollinators. The effective population size of the dove‐tree continuously decreased during the climate changes of the Quaternary. Such climate sensitivity should be fully considered in conservation efforts for this relict endangered species in the context of continuous climate warming in the future.     

Isolation and characterization of polymorphic microsatellite DNA markers in topmouth gudgeon,Pseudorasbora (Teleostei: Cyprinidae)

Due to accidental introductions, the distribution of the east Asian topmouth gudgeon, Pseudorasbora parva, has rapidly expanded over recent years, thereby threatening the endangered Japanese species, P. pumila pumila. Sixteen microsatellite loci were isolated and characterized for P. parva, 14 being highly polymorphic (mean H_E = 0.58, mean number of alleles = 4.4). Successful characterization of 10 of these loci was also achieved for P. pumila pumila. These markers should be useful in future investigation of colonization of P. parva and the development of conservation strategies for P. pumila pumila.     

Emperipolesis mediated by CD8+ T cells correlates with biliary epithelia cell injury in primary biliary cholangitis

Primary biliary cholangitis (PBC) is an autoimmune disease characterized by chronic destruction of the bile ducts. A major unanswered question regarding the pathogenesis of PBC is the precise mechanisms of small bile duct injury. Emperipolesis is one of cell‐in‐cell structures that is a potential histological hallmark associated with chronic hepatitis B. This study aimed to clarify the pathogenesis and characteristics of emperipolesis in PBC liver injury. Sixty‐six PBC patients, diagnosed by liver biopsy combined with laboratory test, were divided into early‐stage PBC (stages I and II, n = 39) and late‐stage PBC (stages III and IV, n = 27). Emperipolesis was measured in liver sections stained with haematoxylin‐eosin. The expressions of CK19, CD3, CD4, CD8, CD20, Ki67 and apoptosis of BECs were evaluated by immunohistochemistry or immunofluorescence double labelling. Emperipolesis was observed in 62.1% of patients with PBC, and BECs were predominantly host cells. The number of infiltrating CD3⁺ and CD8⁺ T cells correlated with the advancement of emperipolesis (R² = 0.318, P < .001; R² = 0.060, P < .05). The cell numbers of TUNEL‐positive BECs and double staining for CK19 and Ki67 showed a significant positive correlation with emperipolesis degree (R² = 0.236, P < .001; R² = 0.267, P < .001). We conclude that emperipolesis mediated by CD8⁺ T cells appears to be relevant to apoptosis of BEC and thus may aggravate the further injury of interlobular bile ducts.     

The relationships between fish size and otolith dimensions in the common sole (Solea solea (Linnaeus, 1758)) captured in the Northeastern Mediterranean

Fish specimens were captured by a commercial bottom trawler at a depth of 50–80 m from Iskenderun Bay (Hatay, Turkey) between December 2017 and May 2018. The bottom trawl gear used was equipped with a 44 mm stretched mesh size net at the cod-end. Blind side and eyed side otolith lengths (OL), otolith breadths (OB) and otolith weights (OW) were measured from each specimen to the nearest 0.001 mm and 0.0001g respectively. A total of 110 fish (43 females and 67 males) were collected. Total length ranged from 20.8 to 28.2 cm and 68.0 to 166.1 g (males) and 21.1 to −28.5 cm and 74.5 to 201.4 g (females). The coefficients of determination between fish weight and otolith weight, and total length and otolith weight (sexes combined) were found as R² = .7694 (0.77) and R² = .6274 (0.63), respectively. A moderate positive relationship between the total length-otolith dimensions, and fish weight-otolith dimensions, was also demonstrated.     

Relationship between semen characteristics and body size in Barbus barbus L. (Teleostei: Cyprinidae) and effects of ions and osmolality on sperm motility

The objectives of the present study were to determine the relationships among length and weight of males, sperm volume, spermatozoa concentration, total number of spermatozoa, ionic contents and osmolality of seminal plasma in Barbus barbus. The effect of osmolality on sperm motility parameters after activation in NaCl, KCl, or sucrose solutions was also examined. There were significant correlations between spermatozoa concentration – length (R = + 0.7) and – weight (R = + 0.8) of males. No significant correlations were observed between the total number of spermatozoa, sperm volume, and length and weight of males. Seminal plasma osmolality was higher when the total number of spermatozoa (R = + 0.6) and sperm volume (R = + 0.6) were higher. Sperm motility and velocity was positively correlated with osmolality (R = + 0.5). The correlation between sperm motility and K⁺ was negative (R = 0.5), but positively correlated with Ca²⁺ (R = 0.8), Na⁺ (R = 0.8), and Cl⁻ (R = 0.8). There was a rapid decrease (P < 0.05) in sperm motility parameters after sperm activation. Just after sperm activation, beating waves propagated along the full length of flagella. At latter stages post sperm activation, the waves appeared only in proximal part of the flagellum. The highest spermatozoa velocity and percentage of motility were observed at 215–235 mOsmol kg^− 1 in NaCl, KCl or sucrose. The tip of the flagellum became curled into a loop shape which shortened the flagellum after activation of sperm in distilled water. B. barbus sperm is very similar to that of other cyprinids in terms of ionic contents and osmolality of the seminal plasma, mechanism of sperm activation and behavior and motility of sperm during swimming period.     

Population genetic structure of island and mainland populations of the quokka, Setonix brachyurus (Macropodidae): a comparison of AFLP and microsatellite markers

Translocation and reintroduction are important tools for the conservation or recovery of species threatened with extinction in the wild. However, an understanding of the potential genetic consequences of mixing populations requires an understanding of the genetic variation within, and similarities among, donor and recipient populations. Genetic diversity was measured using two independent marker systems (microsatellites and AFLPs) for one island and four small remnant mainland populations of Setonix brachyurus, a threatened medium sized macropod restricted to fragmented habitat remnants and two off-shore islands in southwest Australia. Microsatellite diversity in the island population (R _s = 3.2, H _e = 71%) was similar to, or greater than, all mainland populations (R _s = 2.1–3.9, H _e = 34-71%). In contrast, AFLP diversity was significantly lower in the island population (PPL = 20.5; H _j = 0.118) compared to all mainland populations (mean PPL = 79.5–89.7; mean H _j = 0.23–0.29). Microsatellites differentiated all (mainland and island) populations from each other. However, AFLP only differentiated the island population from the mainland populations—all mainland populations were not significantly differentiated from each other for this marker. Given a known time since isolation of the island population from the mainland (6,000 years ago), and an overall more conservative rate of evolution of AFLP markers, our results are consistent with mainland populations fragmenting thousands of years ago (but <6,000 years), probably as a consequence of reduced rainfall and the constriction of the preferred mesic habitat of quokkas. Our results also support a recent history of severe population bottlenecks in mainland populations, and a long history of bottlenecks of the island population, but reflect a recent explosion in numbers since European occupation of the island. Our results indicate that translocation of island populations to supplement mainland populations would introduce genetically markedly differentiated, and possibly maladapted, individuals.