Genomic and transcriptomic insights into molecular basis of sexually dimorphic nuptial spines in Leptobrachium leishanense

Diversity and Molecular Evolution of Nonvisual Opsin Genes across Environmental, Developmental, and Morphological Adaptations in Frogs

Nonvisual opsins are transmembrane proteins expressed in the eyes and other tissues of many animals. When paired with a light-sensitive chromophore, nonvisual opsins form photopigments involved in various nonvisual, light-detection functions including circadian rhythm regulation, light-seeking behaviors, and seasonal responses. Here, we investigate the molecular evolution of nonvisual opsin genes in anuran amphibians (frogs and toads). We test several evolutionary hypotheses including the predicted loss of nonvisual opsins due to nocturnal ancestry and potential functional differences in nonvisual opsins resulting from environmental light variation across diverse anuran ecologies. Using whole-eye transcriptomes of 81 species, combined with genomes, multitissue transcriptomes, and independently annotated genes from an additional 21 species, we identify which nonvisual opsins are present in anuran genomes and those that are also expressed in the eyes, compare selective constraint among genes, and test for potential adaptive evolution by comparing selection between discrete ecological classes. At the genomic level, we recovered all 18 ancestral vertebrate nonvisual opsins, indicating that anurans demonstrate the lowest documented amount of opsin gene loss among ancestrally nocturnal tetrapods. We consistently found expression of 14 nonvisual opsins in anuran eyes and detected positive selection in a subset of these genes. We also found shifts in selective constraint acting on nonvisual opsins in frogs with differing activity periods, habitats, distributions, life histories, and pupil shapes, which may reflect functional adaptation. Although many nonvisual opsins remain poorly understood, these findings provide insight into the diversity and evolution of these genes across anurans, filling an important gap in our understanding of vertebrate opsins and setting the stage for future research on their functional evolution across taxa.

Highly Stretchable Polyurethane Porous Membranes with Adjustable Morphology for Advanced Lithium Metal Batteries

A highly flexible, tunable morphology membrane with excellent thermal stability and ionic conductivity can endow lithium metal batteries with high power density and reduced dendrite growth. Herein, a porous Polyurethane (PU) membrane with an adjustable morphology was prepared by a simple nonsolvent-induced phase separation technique. The precise control of the final morphology of PU membranes can be achieved through appropriate selection of a nonsolvent, resulting a range of pore structures that vary from finger-like voids to sponge-like pores. The implementation of combinatorial DFT and experimental analysis has revealed that spongy PU porous membranes, especially PU-EtOH, show superior electrolyte wettability (472%), high porosity (75%), good mechanical flexibility, robust thermal dimensional stability (above 170 °C), and elevated ionic conductivity (1.38 mS cm-1) in comparison to the polypropylene (PP) separator. The use of PU-EtOH in Li//Li symmetric cell results in a prolonged lifespan of 800 h, surpasing the longevity of PU or PP cells. Moreover, when subjected to a high rate of 5 C, the LiFePO4/Li half-cell with a PU-EtOH porous membrane displayed better cycling performance (115.4 mAh g-1) compared to the PP separator (104.4 mAh g-1). Finally, the prepared PU porous membrane exhibits significant potential for improving the efficiency and safety of LMBs.

Sex Differences in Cochlear Transcriptomes in Horseshoe Bats

Sexual dimorphism of calls is common in animals, whereas studies on the molecular basis underlying this phenotypic variation are still scarce. In this study, we used comparative transcriptomics of cochlea to investigate the sex-related difference in gene expression and alternative splicing in four Rhinolophus taxa. Based on 31 cochlear transcriptomes, we performed differential gene expression (DGE) and alternative splicing (AS) analyses between the sexes in each taxon. Consistent with the degree of difference in the echolocation pulse frequency between the sexes across the four taxa, we identified the largest number of differentially expressed genes (DEGs) and alternatively spliced genes (ASGs) in R. sinicus. However, we also detected multiple DEGs and ASGs in taxa without sexual differences in echolocation pulse frequency, suggesting that these genes might be related to other parameters of echolocation pulse rather than the frequency component. Some DEGs and ASGs are related to hearing loss or deafness genes in human or mice and they can be considered to be candidates associated with the sexual differences of echolocation pulse in bats. We also detected more than the expected overlap of DEGs and ASGs in two taxa. Overall, our current study supports the important roles of both DGE and AS in generating or maintaining sexual differences in animals.

The paddy frog genome provides insight into the molecular adaptations and regulation of hibernation in ectotherms

Amphibians, like the paddy frog (Fejervarya multistriata), have played a critical role in the transition from water to land. Hibernation is a vital survival adaptation in cold environments with limited food resources. We decoded the paddy frog genome to reveal the molecular adaptations linked to hibernation in ectotherms. The genome contained 13 chromosomes, with a significant proportion of repetitive sequences. We identified the key genes encoding the proteins of AANAT, TRPM8, EGLN1, and VEGFA essential for circadian rhythms, thermosensation, and hypoxia during hibernation by comparing the hibernator and non-hibernator genomes. Examining organ changes during hibernation revealed the central regulatory role of the brain. We identified 21 factors contributing to hibernation, involving hormone biosynthesis, protein digestion, DNA replication, and the cell cycle. These findings provide deeper insight into the complex mechanisms of ectothermic hibernation and contribute to our understanding of the broader significance of this evolutionary adaptation.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Tracking the Diversity and Chromosomal Distribution of the Olfactory Receptor Gene Repertoires of Three Anurans Species

Olfaction is a crucial capability for most vertebrates and is realized through olfactory receptors in the nasal cavity. The enormous diversity of olfactory receptors has been created by gene duplication, following a birth-and-death model of evolution. The olfactory receptor genes of the amphibians have received relatively little attention up to now, although recent studies have increased the number of species for which data are available. This study analyzed the diversity and chromosomal distribution of the OR genes of three anuran species (Engystomops pustulosus, Bufo bufo and Hymenochirus boettgeri). The OR genes were identified through searches for homologies, and sequence filtering and alignment using bioinformatic tools and scripts. A high diversity of OR genes was found in all three species, ranging from 917 in B. bufo to 1194 in H. boettgeri, and a total of 2076 OR genes in E. pustulosus. Six OR groups were recognized using an evolutionary gene tree analysis. While E. pustulosus has one of the highest numbers of genes of the gamma group (which detect airborne odorants) yet recorded in an anuran, B. bufo presented the smallest number of pseudogene sequences ever identified, with no pseudogenes in either the beta or epsilon groups. Although H. boettgeri shares many morphological adaptations for an aquatic lifestyle with Xenopus, and presented a similar number of genes related to the detection of water-soluble odorants, it had comparatively far fewer genes related to the detection of airborne odorants. This study is the first to describe the complete OR repertoire of the three study species and represents an important contribution to the understanding of the evolution and function of the sense of smell in vertebrates.

Association between disease-modifying antirheumatic drugs and bone turnover biomarkers

Rheumatoid arthritis (RA) has been linked to an increased risk of osteoporosis as well as fractures. Patients diagnosed with RA had a 25% increased risk of osteoporotic fracture, according to a recent population-based cohort study that compared them to people without RA. Several studies have found a correlation between osteoporosis and the presence of pro-inflammatory cytokines, such as tumor necrosis factor (TNF)-α, interleukin (IL)-1, and 6. These cytokines play a crucial part in the process of bone resorption by boosting osteoclast activation and encouraging osteoclast differentiation. Based on the correlation between RA, osteoporosis, and inflammation, it is possible that systemic immunosuppression with disease-modifying antirheumatic drugs (DMARDs) can help individuals with RA have a lower chance of developing osteoporosis and osteoporotic fractures. There is little information on how different DMARDs, biologic or non-biologic, affect RA patients' bone metabolism. In this study, we present an overview of the influence that targeted therapies, such as biologics, non-biologics, and small molecule inhibitors, have on bone homeostasis in RA patients.

Chromosome-level genome assembly of a high-altitude-adapted frog (Rana kukunoris) from the Tibetan plateau provides insight into amphibian genome evolution and adaptation

BACKGROUND

The high-altitude-adapted frog Rana kukunoris, occurring on the Tibetan plateau, is an excellent model to study life history evolution and adaptation to harsh high-altitude environments. However, genomic resources for this species are still underdeveloped constraining attempts to investigate the underpinnings of adaptation.

RESULTS

The R. kukunoris genome was assembled to a size of 4.83 Gb and the contig N50 was 1.80 Mb. The 6555 contigs were clustered and ordered into 12 pseudo-chromosomes covering ~ 93.07% of the assembled genome. In total, 32,304 genes were functionally annotated. Synteny analysis between the genomes of R. kukunoris and a low latitude species Rana temporaria showed a high degree of chromosome level synteny with one fusion event between chr11 and chr13 forming pseudo-chromosome 11 in R. kukunoris. Characterization of features of the R. kukunoris genome identified that 61.5% consisted of transposable elements and expansions of gene families related to cell nucleus structure and taste sense were identified. Ninety-five single-copy orthologous genes were identified as being under positive selection and had functions associated with the positive regulation of proteins in the catabolic process and negative regulation of developmental growth. These gene family expansions and positively selected genes indicate regions for further interrogation to understand adaptation to high altitude.

CONCLUSIONS

Here, we reported a high-quality chromosome-level genome assembly of a high-altitude amphibian species using a combination of Illumina, PacBio and Hi-C sequencing technologies. This genome assembly provides a valuable resource for subsequent research on R. kukunoris genomics and amphibian genome evolution in general.

Draft genome and multi-tissue transcriptome assemblies of the Neotropical leaf-frog Phyllomedusa bahiana

Amphibians are increasingly threatened worldwide, but the availability of genomic resources that could be crucial for implementing informed conservation practices lags well behind that for other vertebrate groups. Here, we describe draft de novo genome, mitogenome, and transcriptome assemblies for the Neotropical leaf-frog Phyllomedusa bahiana native to the Brazilian Atlantic Forest and Caatinga. We used a combination of PacBio long reads and Illumina sequencing to produce a 4.74-Gbp contig-level genome assembly, which has a contiguity comparable to other recent nonchromosome level assemblies. The assembled mitogenome comprises 16,239 bp and the gene content and arrangement are similar to other Neobratrachia. RNA-sequencing from 8 tissues resulted in a highly complete (86.3%) reference transcriptome. We further use whole-genome resequencing data from P. bahiana and from its sister species Phyllomedusa burmeisteri, to demonstrate how our assembly can be used as a backbone for population genomics studies within the P. burmeisteri species group. Our assemblies thus represent important additions to the catalog of genomic resources available from amphibians.

Whole-genome sequencing identifies new candidate genes for nonobstructive azoospermia

BACKGROUND

Genetic causes that lead to spermatogenetic failure in patients with nonobstructive azoospermia (NOA) have not been yet completely established.

OBJECTIVE

To identify low-frequency NOA-associated single nucleotide variants (SNVs) using whole-genome sequencing (WGS).

MATERIALS AND METHODS

Men with various types of NOA (n = 39), including samples that had been previously tested with whole-exome sequencing (WES; n = 6) and did not result in diagnostic conclusions. Variants were annotated using the Ensembl Variant Effect Predictor, utilizing frequencies from GnomAD and other databases to provide clinically relevant information (ClinVar), conservation scores (phyloP), and effect predictions (i.e., MutationTaster). Structural protein modeling was also performed.

RESULTS

Using WGS, we revealed potential NOA-associated SNVs, such as: TKTL1, IGSF1, ZFPM2, VCX3A (novel disease causing variants), ESX1, TEX13A, TEX14, DNAH1, FANCM, QRICH2, FSIP2, USP9Y, PMFBP1, MEI1, PIWIL1, WDR66, ZFX, KCND1, KIAA1210, DHRSX, ZMYM3, FAM47C, FANCB, FAM50B (genes previously known to be associated with infertility) and ALG13, BEND2, BRWD3, DDX53, TAF4, FAM47B, FAM9B, FAM9C, MAGEB6, MAP3K15, RBMXL3, SSX3 and FMR1NB genes, which may be involved in spermatogenesis.

DISCUSSION AND CONCLUSION

In this study, we identified novel potential candidate NOA-associated genes in 29 individuals out of 39 azoospermic males. Note that in 5 out of 6 patients subjected previously to WES analysis, which did not disclose potentially causative variants, the WGS analysis was successful with NOA-associated gene findings.

Nanomedicine and versatile therapies for cancer treatment

The higher prevalence of cancer is related to high rates of mortality and morbidity worldwide. By virtue of the properties of matter at the nanoscale, nanomedicine is proven to be a powerful tool to develop innovative drug carriers with greater efficacies and fewer side effects than conventional therapies. In this review, different nanocarriers for controlled drug release and their routes of administration have been discussed in detail, especially for cancer treatment. Special emphasis has been given on the design of drug delivery vehicles for sustained release and specific application methods for targeted delivery to the affected areas. Different polymeric vehicles designed for the delivery of chemotherapeutics have been discussed, including graft copolymers, liposomes, hydrogels, dendrimers, micelles, and nanoparticles. Furthermore, the effect of dimensional properties on chemotherapy is vividly described. Another integral section of the review focuses on the modes of administration of nanomedicines and emerging therapies, such as photothermal, photodynamic, immunotherapy, chemodynamic, and gas therapy, for cancer treatment. The properties, therapeutic value, advantages, and limitations of these nanomedicines are highlighted, with a focus on their increased performance versus conventional molecular anticancer therapies.

Sequencing of laser captured Z and W chromosomes of the tocantins paradoxical frog (Pseudis tocantins) provides insights on repeatome and chromosomal homology

Pseudis tocantins is the only frog species of the hylid genus Pseudis that possesses highly heteromorphic sex chromosomes. Z and W chromosomes of Ps. tocantins differ in size, morphology, position of the nucleolar organizer region (NOR) and the amount and distribution of heterochromatin. A chromosomal inversion and heterochromatin amplification on the W chromosome were previously inferred to be involved in the evolution of this sex chromosome pair. Despite these findings, knowledge related to the molecular composition of the large heterochromatic band of this W chromosome is restricted to the PcP190 satellite DNA, and no data are available regarding the gene content of either the W or the Z chromosome of Ps. tocantins. Here, we sequenced microdissected Z and W chromosomes of this species to further resolve their molecular composition. Comparative genomic analysis suggests that Ps. tocantins sex chromosomes are likely homologous to chromosomes 4 and 10 of Xenopus tropicalis. Analyses of the repetitive DNA landscape in the Z and W assemblies allowed for the identification of several transposable elements and putative satellite DNA sequences. Finally, some transposable elements from the W assembly were found to be highly diverse and divergent from elements found elsewhere in the genome, suggesting a rapid amplification of these elements on the W chromosome.

Effect of green synthesized ZnO nanoparticles using Paspalum scrobiculatum grains extract in biological applications

In this report, ZnO nanoparticles were biosynthesized using Paspalum scrobiculatum grains extract for the first time. GC-MS analysis explicated that diethyl phthalate was the major phytocompound with 94.09% in aqueous extract. ZnO nanoparticles formation was confirmed by various physicochemical analyses. HR-TEM images showed the hexagonal, rectangular shaped nanoparticles in 15-30 nm size. The antioxidant, anti-inflammatory, and anti-diabetic analyses showed the effective bioactivity of ZnO nanoparticles in 80 μg/ml concentration with 95.36%, 94.08%, and 91.96%, respectively. The morphological and tissue changes witnessed in larvicidal and insecticidal activities against Culex tritaeniorhynchus and Tribolium castaneum revealed the efficient nature of ZnO nanoparticles in 100 ppm at 48 h and 100 μg/kg at 72 h, respectively. The morphological changes in antibacterial activity demonstrated the bactericidal nature of ZnO nanoparticles against Salmonella typhi and Staphylococcus aureus in 150 μg/ml concentration. The morphological observations in anticancer activity against HepG2 liver cancer cells showed the potent drug features of ZnO nanoparticles in 100 μg/ml concentration with 97.18% of cytotoxicity. The ZnO nanoparticles showed no toxicity against HDF normal cells in lower concentrations and it explicated the biocompatible features of nanoparticles. The Vigna radiata plant growth was efficiently promoted by low (60 ppm) concentration of nanoparticles. The ZnO nanoparticles divulged effective degradation of IPA, EDTA, BQ, and DPBF in 75%, 45%, 55%, and 80% through ROS formation, respectively. Thus, the synthesized ZnO nanoparticles are biocompatible and inexpensive material compared to the traditional one and can be utilized as an efficient material in biological fields.

Genomic adaptations for arboreal locomotion in Asian flying treefrogs

SignificanceTo adapt to arboreal lifestyles, treefrogs have evolved a suite of complex traits that support vertical movement and gliding, thus presenting a unique case for studying the genetic basis for traits causally linked to vertical niche expansion. Here, based on two de novo-assembled Asian treefrog genomes, we determined that genes involved in limb development and keratin cytoskeleton likely played a role in the evolution of their climbing systems. Behavioral and morphological evaluation and time-ordered gene coexpression network analysis revealed the developmental patterns and regulatory pathways of the webbed feet used for gliding in Rhacophorus kio.

Evolutionary analyses of visual opsin genes in frogs and toads: Diversity, duplication, and positive selection

Among major vertebrate groups, anurans (frogs and toads) are understudied with regard to their visual systems, and little is known about variation among species that differ in ecology. We sampled North American anurans representing diverse evolutionary and life histories that likely possess visual systems adapted to meet different ecological needs. Using standard molecular techniques, visual opsin genes, which encode the protein component of visual pigments, were obtained from anuran retinas. Additionally, we extracted the visual opsins from publicly available genome and transcriptome assemblies, further increasing the phylogenetic and ecological diversity of our dataset to 33 species in total. We found that anurans consistently express four visual opsin genes (RH1, LWS, SWS1, and SWS2, but not RH2) even though reported photoreceptor complements vary widely among species. The proteins encoded by these genes showed considerable sequence variation among species, including at sites known to shift the spectral sensitivity of visual pigments in other vertebrates and had conserved substitutions that may be related to dim-light adaptation. Using molecular evolutionary analyses of selection (dN/dS) we found significant evidence for positive selection at a subset of sites in the dim-light rod opsin gene RH1 and the long wavelength sensitive cone opsin LWS. The function of sites inferred to be under positive selection are largely unknown, but a few are likely to affect spectral sensitivity and other visual pigment functions based on proximity to previously identified sites in other vertebrates. We also found the first evidence of visual opsin duplication in an amphibian with the duplication of the LWS gene in the African bullfrog, which had distinct LWS copies on the sex chromosomes suggesting the possibility of sex-specific visual adaptation. Taken together, our results indicate that ecological factors, such as habitat and life history, as well as behavior, may be driving changes to anuran visual systems.

Genetic Diversity and Population Dynamics of Leptobrachium leishanense (Anura: Megophryidae) as Determined by Tetranucleotide Microsatellite Markers Developed from Its Genome

Persisting declination of amphibians around the world has resulted in the public attaching importance to the conservation of their biodiversity. Genetic data can be greatly helpful in conservation planning and management, especially in species that are small in size and hard to observe. It is essential to perform genetic assessments for the conservation of Leptobrachium leishanense, an endangered toad and receiving secondary protection on the list of state-protected wildlife in China. However, current molecular markers with low reliability and efficiency hinder studies. Here, we sampled 120 adult toes from the population in the Leishan Mountain, 23 of which were used to develop tetranucleotide microsatellite markers based on one reference L. leishanense genome. After primer optimization, stability detection, and polymorphism detection, we obtained 12 satisfactory microsatellite loci. Then, we used these loci to evaluate the genetic diversity and population dynamics of the 120 individuals. Our results show that there is a low degree of inbreeding in the population, and it has a high genetic diversity. Recently, the population has not experienced population bottlenecks, and the estimated effective population size was 424.3. Accordingly, stabilizing genetic diversity will be key to population sustainability. Recovering its habitat and avoiding intentional human use will be useful for conservation of this species.

A dense linkage map for a large repetitive genome: discovery of the sex-determining region in hybridizing fire-bellied toads (Bombina bombina and Bombina variegata)

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7–10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

Pseudogenized Amelogenin Reveals Early Tooth Loss in True Toads (Anura: Bufonidae)

Extant anurans (frogs and toads) exhibit reduced dentition, ranging from a lack of mandibular teeth to complete edentulation, as observed in the true toads of the family Bufonidae. The evolutionary time line of these reductions remains vague due to a poor fossil record. Previous studies have demonstrated an association between the lack of teeth in edentulous vertebrates and the pseudogenization of the major tooth enamel gene amelogenin (AMEL) through accumulation of deleterious mutations and the disruption of its coding sequence. In this study, we have harnessed the pseudogenization of AMEL as a molecular dating tool to correlate loss of dentition with genomic mutation patterns during the rise of the family Bufonidae. Specifically, we have utilized AMEL pseudogenes in three members of the family as a tool to estimate the putative date of edentulation in true toads. Comparison of AMEL sequences from Rhinella marina, Bufo gargarizans and Bufo bufo, with nine extant, dentulous frogs, revealed mutations confirming AMEL inactivation in Bufonidae. AMEL pseudogenes in modern bufonids also exhibited remarkably high 86-93% sequence identity among each other, with only a slight increase in substitution rate and relaxation of selective pressure, in comparison with functional copies in other anurans. Moreover, using selection intensity estimates and synonymous substitution rates, analysis of functional and pseudogenized AMEL resulted in an estimated inactivation window of 46-60 million years ago in the lineage leading to modern true toads, a time line that coincides with the rise of the family Bufonidae.

Transcriptomic Profiling Reveals the Antiapoptosis and Antioxidant Stress Effects of Fos in Ischemic Stroke

Arterial hypertension is considered the most prevalent risk factor for stroke. Both pathophysiologic and clinical data previously acquired suggest a strong correlation between the hemodynamic nature of arterial hypertension and an increase in the risk of ischemic insult to tissues. However, the knowledge of specific molecular interactions between hypertension and ischemic stroke (IS) is limited. In this study, we performed systematic bioinformatics analysis of stroke-prone spontaneous hypertensive brain tissue samples of rats (GSE41452), middle cerebral artery occlusion of brain tissue samples of rats (GSE97537), and peripheral blood array data of IS patients (GSE22255). We identified that Fos, an immediate-early gene (IEG) that responds to alterations in arterial blood pressure, has a strong correlation with the occurrence and prognosis of IS. To further evaluate the potential function of Fos, the oxygen–glucose deprivation model and RNA sequencing of HT22 neuronal cells were performed. Consistent with the sequencing results, real-time quantitative PCR and Western blot indicate that Fos was elevated at 3 h and returned to normal levels at 6 h after oxygen–glucose deprivation. Knock-down of Fos by lentivirus significantly increased the oxidative stress level, neuronal apoptosis, and inhibited the mitochondrial function. In conclusion, Fos acts as an important link between hypertension and IS. Furthermore, Fos can be used as a potential biomarker for target therapy in the prevention of stroke among hypertensive patients and also potential treatment targeting apoptosis and oxidative stress after its onset.

Coding and Noncoding Variation in LRRK2 and Parkinson's Disease Risk

Background

The leucine‐rich repeat kinase 2 (LRRK2) gene harbors both rare highly damaging missense variants (eg, p.G2019S) and common noncoding variants (eg, rs76904798) with lower effect sizes that are associated with Parkinson's disease (PD) risk.

Objectives

This study aimed to investigate in a large meta‐analysis whether the LRRK2 Genome‐Wide Association Study (GWAS) signal represented by rs76904798 is independently associated with PD risk from LRRK2 coding variation and whether complex linkage disequilibrium structures with p.G2019S and the 5′ noncoding haplotype account for the association of LRRK2 coding variants.

Methods

We performed a meta‐analysis using imputed genotypes from 17,838 patients, 13,404 proxy patients, and 173,639 healthy controls of European ancestry. We excluded carriers of p.G2019S and/or rs76904798 to clarify the role of LRRK2 coding variation in mediating disease risk and excluded carriers of relatively rare LRRK2 coding variants to assess the independence of rs76904798. We also investigated the co‐inheritance of LRRK2 coding variants with p.G2019S, rs76904798, and p.N2081D.

Results

LRRK2 rs76904798 remained significantly associated with PD after excluding the carriers of relatively rare LRRK2 coding variants. LRRK2 p.R1514Q and p.N2081D were frequently co‐inherited with rs76904798, and the allele distribution of p.S1647T significantly changed among patients after removing rs76904798 carriers.

Conclusions

These data suggest that the LRRK2 coding variants previously related to PD (p.N551K, p.R1398H, p.M1646T, and p.N2081D) do not drive the 5′ noncoding GWAS signal. These data, however, do not preclude the independent association of the haplotype p.N551K‐p.R1398H and p.M1646T with altered disease risk. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson Movement Disorder Society. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

A brief review of vertebrate sex evolution with a pledge for integrative research: towards 'sexomics'

Triggers and biological processes controlling male or female gonadal differentiation vary in vertebrates, with sex determination (SD) governed by environmental factors or simple to complex genetic mechanisms that evolved repeatedly and independently in various groups. Here, we review sex evolution across major clades of vertebrates with information on SD, sexual development and reproductive modes. We offer an up-to-date review of divergence times, species diversity, genomic resources, genome size, occurrence and nature of polyploids, SD systems, sex chromosomes, SD genes, dosage compensation and sex-biased gene expression. Advances in sequencing technologies now enable us to study the evolution of SD at broader evolutionary scales, and we now hope to pursue a sexomics integrative research initiative across vertebrates. The vertebrate sexome comprises interdisciplinary and integrated information on sexual differentiation, development and reproduction at all biological levels, from genomes, transcriptomes and proteomes, to the organs involved in sexual and sex-specific processes, including gonads, secondary sex organs and those with transcriptional sex-bias. The sexome also includes ontogenetic and behavioural aspects of sexual differentiation, including malfunction and impairment of SD, sexual differentiation and fertility. Starting from data generated by high-throughput approaches, we encourage others to contribute expertise to building understanding of the sexomes of many key vertebrate species. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Genome of the butterfly hillstream loach provides insights into adaptations to torrential mountain stream life

Butterfly hillstream loach (Beaufortia kweichowensis), a benthic fish in the torrential mountain streams, possesses a totally flat ventrum, flattened craniofacial and body skeletons, and enlarged paired fins covered by substantially small keratinous structures. However, little is known about the genetic basis of these specialized morphological adaptations. Here we present a 448.52-Mb genome assembly with contig N50 length of 5.53 Mb by integrating Illumina short-read sequencing, Nanopore long-read sequencing and HiC-based chromatin map. Demographic history reconstruction of the butterfly hillstream loach reveals that the population dynamics is correlated with the different stages of uplifting of the Tibetan Plateau. Comparative genomic analysis finds evidence of six keratin genes in butterfly hillstream loach evolving under positive selection. Within these genes, two keratin genes exhibit species-specific and divergent amino acid changes, suggesting a role in the formation of the unculi. Additionally, a series of positively selected genes, rapid evolving genes, specific variant genes and expanded gene families are found, including genes related to Hedgehog, Notch and BMP pathways, which may be involved in craniofacial development. These findings may have important implications for understanding the genetic basis of phenotypic adaptation to torrential mountain stream life.

Rapid genetic adaptation to recently colonized environments is driven by genes underlying life history traits

Background

Uncovering the mechanisms underlying rapid genetic adaptation can provide insight into adaptive evolution and shed light on conservation, invasive species control, and natural resource management. However, it can be difficult to experimentally explore rapid adaptation due to the challenges associated with propagating and maintaining species in captive environments for long periods of time. By contrast, many introduced species have experienced strong selection when colonizing environments that differ substantially from their native range and thus provide a “natural experiment” for studying rapid genetic adaptation. One such example occurred when sea lamprey (Petromyzon marinus), native to the northern Atlantic, naturally migrated into Lake Champlain and expanded their range into the Great Lakes via man-made shipping canals.

Results

Utilizing 368,886 genome-wide single nucleotide polymorphisms (SNPs), we calculated genome-wide levels of genetic diversity (i.e., heterozygosity and π) for sea lamprey collected from native (Connecticut River), native but recently colonized (Lake Champlain), and invasive (Lake Michigan) populations, assessed genetic differentiation between all populations, and identified candidate genes that responded to selection imposed by the novel environments. We observed a 14 and 24% reduction in genetic diversity in Lake Michigan and Lake Champlain populations, respectively, compared to individuals from the Connecticut River, suggesting that sea lamprey populations underwent a genetic bottleneck during colonization. Additionally, we identified 121 and 43 outlier genes in comparisons between Lake Michigan and Connecticut River and between Lake Champlain and Connecticut River, respectively. Six outlier genes that contained synonymous SNPs in their coding regions and two genes that contained nonsynonymous SNPs may underlie the rapid evolution of growth (i.e., GHR), reproduction (i.e., PGR, TTC25, STARD10), and bioenergetics (i.e., OXCT1, PYGL, DIN4, SLC25A15).

Conclusions

By identifying the genomic basis of rapid adaptation to novel environments, we demonstrate that populations of invasive species can be a useful study system for understanding adaptive evolution. Furthermore, the reduction in genome-wide levels of genetic diversity associated with colonization coupled with the identification of outlier genes underlying key life history traits known to have changed in invasive sea lamprey populations (e.g., growth, reproduction) illustrate the utility in applying genomic approaches for the successful management of introduced species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07553-x.

Comparative genomics of Glandirana rugosa using unsupervised AI reveals a high CG frequency

Genome sequencing of a unique frog (Glandirana rugosa) having XY/ZW systems within the species and comparative genomics with other six frogs were performed using a batch-learning self-organizing map, which is unsupervised AI for oligonucleotide compositions, to clarify its genome characteristics.

The Japanese wrinkled frog (Glandirana rugosa) is unique in having both XX-XY and ZZ-ZW types of sex chromosomes within the species. The genome sequencing and comparative genomics with other frogs should be important to understand mechanisms of turnover of sex chromosomes within one species or during a short period. In this study, we analyzed the newly sequenced genome of G. rugosa using a batch-learning self-organizing map which is unsupervised artificial intelligence for oligonucleotide compositions. To clarify genome characteristics of G. rugosa, we compared its short oligonucleotide compositions in all 1-Mb genomic fragments with those of other six frog species (Pyxicephalus adspersus, Rhinella marina, Spea multiplicata, Leptobrachium leishanense, Xenopus laevis, and Xenopus tropicalis). In G. rugosa, we found an Mb-level large size of repeat sequences having a high identity with the W chromosome of the African bullfrog (P. adspersus). Our study concluded that G. rugosa has unique genome characteristics with a high CG frequency, and its genome is assumed to heterochromatinize a large size of genome via methylataion of CG.

A large genome with chromosome-scale assembly sheds light on the evolutionary success of a true toad (Bufo gargarizans)

We present a high-quality genome assembly for the Asiatic toad (Bufo gargarizans) and explore the evolution of several large gene families in amphibians. With a large genome assembly size of 4.55 Gb, the chromosome-scale assembly includes 747 scaffolds with an N50 of 539.8 Mb and 1.79% gaps. Long terminal repeats (LTRs) constitute a high proportion of the genome and their expansion is a key contributor to the inflated genome size in this species. This is very different from other small amphibian genomes, but similar to that of the enormous axolotl genome. The genome retains a large number of duplicated genes, with tandem (TD) and proximal duplications (PD) the predominant mode of duplication. A total of 122 gene families have undergone significant expansion and were mainly enriched in sensory perception of smell and bitter taste. The CYP2C subfamily, which plays an important role in metabolic detoxification, specifically expanded via TD and PD in the Asiatic toad and the cane toad (true toads). Most of Na⁺ /K⁺ -ATPase genes experienced accelerated evolution along Bufonid lineages and two amino acid sites involving toad-toxin resistance were found to experience positive selection. We also revealed a dynamic evolution of olfactory and vomeronasal receptor gene families which was likely driven by the water-to-land transition. The high-quality genome of the Asiatic toad will provide a solid foundation to understand the genetic basis of its many biological processes.

Perspectives on studying molecular adaptations of amphibians in the genomic era

Understanding the genetic mechanisms underlying particular adaptations/phenotypes of organisms is one of the core issues of evolutionary biology. The use of genomic data has greatly advanced our understandings on this issue, as well as other aspects of evolutionary biology, including molecular adaptation, speciation, and even conservation of endangered species. Despite the well-recognized advantages, usages of genomic data are still limited to non-mammal vertebrate groups, partly due to the difficulties in assembling large or highly heterozygous genomes. Although this is particularly the case for amphibians, nonetheless, several comparative and population genomic analyses have shed lights into the speciation and adaptation processes of amphibians in a complex landscape, giving a promising hope for a wider application of genomics in the previously believed challenging groups of organisms. At the same time, these pioneer studies also allow us to realize numerous challenges in studying the molecular adaptations and/or phenotypic evolutionary mechanisms of amphibians. In this review, we first summarize the recent progresses in the study of adaptive evolution of amphibians based on genomic data, and then we give perspectives regarding how to effectively identify key pathways underlying the evolution of complex traits in the genomic era, as well as directions for future research.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

A high-quality genome of taro (Colocasia esculenta (L.) Schott), one of the world's oldest crops

Taro (Colocasia esculenta (L.), Schott), from the Araceae family, is one of the oldest crops with important edible, medicinal, nutritional and economic value. Taro is a highly polymorphic species including diverse genotypes adapted to a broad range of environments, but the taro genome has rarely been investigated. Here, a high-quality chromosome-level genome of C. esculenta was assembled using data sequenced by Illumina, PacBio and Nanopore platforms. The assembled genome size was 2,405 Mb with a contig N50 of 400.0 kb and a scaffold N50 of 159.4 Mb. In total, 2,311 Mb (96.09%) of the contig sequences was anchored onto 14 chromosomes to form pseudomolecules, and 2,126 Mb (88.43%) was annotated as repetitive sequences. Of the 28,695 predicted protein-coding genes, 26,215 genes (91.4%) could be functionally annotated. On the basis of phylogenetic analysis using 769 genes, C. esculenta and Spirodela polyrhiza were placed on one branch of the tree that diverged approximately 73.23 million years ago. The synteny analyses showed that there have been two whole-genome duplication events in C. esculenta separated by a relatively short gap. According to comparative genome analysis, a larger number (1,189) of distinct gene families and long terminal repeats were enriched in C. esculenta. Our high-quality taro genome will provide valuable resources for further genetic, ecological and evolutionary analyses of taro or other species in the Araceae.

Transcriptional Inhibition of Sp-IAG by Crustacean Female Sex Hormone in the Mud Crab, Scylla paramamosain

In crustaceans, the regulation of sex differentiation is mediated by insulin-like androgenic hormone (IAG) and crustacean female sex hormone (CFSH). CFSH is reported to inhibit IAG gene (Sp-IAG) expression in the mud crab Scylla paramamosain, but the regulatory mechanism is not well understood. A 2674 bp 5′ flanking Sp-IAG contains many potential transcription factor binding sites. In this study, analysis of serially deleted 5′ flanking Sp-IAG and site-directed mutation (SDM) of transcription factor binding sites of the same gene showed that the promoter activity of reporter vectors with Sox-5-binding site, signal transducers and activators of transcription (STAT)-binding site and activator protein 1 (AP-1)-binding site were significantly higher than that of vectors without these regions, suggesting that they were involved in transcriptional regulation of Sp-IAG expression. The expression analysis of these transcription factor showed that there was no difference in the level of mRNA in Sox-5 and AP-1 in androgenic gland treated with recombinant CFSH, but expression of Sp-STAT was significantly reduced, suggesting that CFSH regulates the expression of Sp-STAT, inhibiting its function to regulate Sp-IAG. Further experiment revealed that RNAi mediated Sp-STAT gene knockdown reduced the expression of Sp-IAG. These results suggested that Sp-CFSH regulates Sp-IAG by inhibiting STAT. This is a pioneering finding on the transcriptional mechanism of IAG gene in crustaceans.

A draft genome assembly of the eastern banjo frog Limnodynastes dumerilii dumerilii (Anura: Limnodynastidae)

Amphibian genomes are usually challenging to assemble due to their large genome size and high repeat content. The Limnodynastidae is a family of frogs native to Australia, Tasmania and New Guinea. As an anuran lineage that successfully diversified on the Australian continent, it represents an important lineage in the amphibian tree of life but lacks reference genomes. Here we sequenced and annotated the genome of the eastern banjo frog Limnodynastes dumerilii dumerilii to fill this gap. The total length of the genome assembly is 2.38 Gb with a scaffold N50 of 285.9 kb. We identified 1.21 Gb of non-redundant sequences as repetitive elements and annotated 24,548 protein-coding genes in the assembly. BUSCO assessment indicated that more than 94% of the expected vertebrate genes were present in the genome assembly and the gene set. We anticipate that this annotated genome assembly will advance the future study of anuran phylogeny and amphibian genome evolution.

CSA: A high-throughput chromosome-scale assembly pipeline for vertebrate genomes

BACKGROUND

Easy-to-use and fast bioinformatics pipelines for long-read assembly that go beyond the contig level to generate highly continuous chromosome-scale genomes from raw data remain scarce.

RESULT

Chromosome-Scale Assembler (CSA) is a novel computationally highly efficient bioinformatics pipeline that fills this gap. CSA integrates information from scaffolded assemblies (e.g., Hi-C or 10X Genomics) or even from diverged reference genomes into the assembly process. As CSA performs automated assembly of chromosome-sized scaffolds, we benchmark its performance against state-of-the-art reference genomes, i.e., conventionally built in a laborious fashion using multiple separate assembly tools and manual curation. CSA increases the contig lengths using scaffolding, local re-assembly, and gap closing. On certain datasets, initial contig N50 may be increased up to 4.5-fold. For smaller vertebrate genomes, chromosome-scale assemblies can be achieved within 12 h using low-cost, high-end desktop computers. Mammalian genomes can be processed within 16 h on compute-servers. Using diverged reference genomes for fish, birds, and mammals, we demonstrate that CSA calculates chromosome-scale assemblies from long-read data and genome comparisons alone. Even contig-level draft assemblies of diverged genomes are helpful for reconstructing chromosome-scale sequences. CSA is also capable of assembling ultra-long reads.

CONCLUSIONS

CSA can speed up and simplify chromosome-level assembly and significantly lower costs of large-scale family-level vertebrate genome projects.