Validation and development of eDNA metabarcoding primers for comprehensive assessment of Chinese amphibians.

Environmental DNA (eDNA) metabarcoding has emerged as a powerful, non-invasive tool for biodiversity assessments. However, the accuracy and limitations of these assessment techniques are highly dependent on the choice of primer pairs being used. Although several primer sets have been used in eDNA metabarcoding studies of amphibians, there are few comparisons of their reliability and efficiency. Here, we employed lab- and field-tested sets of publicly available and de novo-designed primers in amplifying 83 species of amphibian from all three orders (Anura, Caudata, and Gymnophiona) and 13 families present in China to evaluate the versatility and specificity of these primers sets in amphibian eDNA metabarcoding studies. Three pairs of primers were highly effective, as they could successfully amplify all the major clades of Chinese amphibians in our study. A few non-amphibian taxa were also amplified by these primers, which implies that further optimization of amphibian-specific primers is still needed. The simultaneous use of three primer sets can completely cover all the species obtained by conventional survey methods and has even effectively distinguished quite a number of species (n = 20) in the Wenshan National Nature Reserve. No single primer set could individually detect all of the species from the studied region, indicating that multiple primers might be necessary for a comprehensive survey of Chinese amphibians. Besides, seasonal variations in amphibian species composition were also revealed by eDNA metabarcoding, which was consistent with traditional survey methods. These results indicate that eDNA metabarcoding has the potential to be a powerful tool for studying spatial and temporal community changes in amphibian species richness.

Interactions between plants and bacterial communities for phytoremediation of petroleum-contaminated soil.

Plants can stimulate the microbes to degrade ubiquitous petroleum hydrocarbons (PHCs), which has prompted a novel view on rhizoremediation. In the present study, the degradation rate of PHCs was investigated and 16S rRNA gene analysis was performed to investigate the PHC-degrading bacteria in petroleum-contaminated soil with different plants. Mirabilis jalapa (M. jalapa) has a higher PHC degradation rate than Lolium perenne (L. perenne) under petroleum contamination. The bacterial diversity in rhizospheric soil was decreased but the relative abundance of Actinobacteriota, Proteobacteria, and Candidatus Saccharibacteria were significant increased on 45 days petroleum-contaminated rhizospheric soil. In addition, the relative expression of PHC degradation-related genes, the content of malic acid and citric acid of the root exudates in the two plants was significantly increased in response to petroleum stress. The content of citric acid increased 11.9 times in M. jalapa and 3.4 times in L. perenne, respectively, in response to petroleum stress. These results indicate that M. jalapa changes the hydrocarbon-degrading microbial community to enhance the degradation of PHCs by root exudates and phytostimulation.

Integrated transcriptomics and metabolomics analyses of the effects of bagging treatment on carotenoid biosynthesis and regulation of Areca catechu L.

Introduction: Fresh Aareca nut fruit for fresh fruit chewing commonly found in green or dark green hues. Despite its economic significance, there is currently insufficient research on the study of color and luster of areca. And the areca nut fruits after bagging showed obvious color change from green to tender yellow. In the study, we tried to explain this interesting variation in exocarp color. Methods: Fruits were bagged (with a double-layered black interior and yellow exterior) 45 days after pollination and subsequently harvested 120 days after pollination. In this study, we examined the the chlorophyll and carotenoid content of pericarp exocarp, integrated transcriptomics and metabolomics to study the effects of bagging on the carotenoid pathway at the molecular level. Results: It was found that the chlorophyll and carotenoid content of bagged areca nut (YP) exocarp was significantly reduced. A total of 21 differentially expressed metabolites (DEMs) and 1784 differentially expressed genes (DEGs) were screened by transcriptomics and metabolomics. Three key genes in the carotenoid biosynthesis pathway as candidate genes for qPCR validation by co-analysis, which suggested their role in the regulation of pathways related to crtB, crtZ and CYP707A. Discussion: We described that light intensity may appear as a main factor influencing the noted shift from green to yellow and the ensuing reduction in carotenoid content after bagging.

Hidden hotspots of amphibian biodiversity in China.

Identifying and protecting hotspots of endemism and species richness is crucial for mitigating the global biodiversity crisis. However, our understanding of spatial diversity patterns is far from complete, which severely limits our ability to conserve biodiversity hotspots. Here, we report a comprehensive analysis of amphibian species diversity in China, one of the most species-rich countries on Earth. Our study combines 20 y of field surveys with new molecular analyses of 521 described species and also identifies 100 potential cryptic species. We identify 10 hotspots of amphibian diversity in China, each with exceptional species richness and endemism and with exceptional phylogenetic diversity and phylogenetic endemism (based on a new time-calibrated, species-level phylogeny for Chinese amphibians). These 10 hotspots encompass 59.6% of China's described amphibian species, 49.0% of cryptic species, and 55.6% of species endemic to China. Only four of these 10 hotspots correspond to previously recognized biodiversity hotspots. The six new hotspots include the Nanling Mountains and other mountain ranges in South China. Among the 186 species in the six new hotspots, only 9.7% are well covered by protected areas and most (88.2%) are exposed to high human impacts. Five of the six new hotspots are under very high human pressure and are in urgent need of protection. We also find that patterns of richness in cryptic species are significantly related to those in described species but are not identical.

Genome Sequencing of a Gray Wolf from Peninsular India Provides New Insights into the Evolution and Hybridization of Gray Wolves.

The gray wolf (Canis lupus) is among the few large carnivores that survived the Late Pleistocene megafaunal extinctions. Thanks to their complex history of admixture and extensive geographic range, the number of gray wolf subspecies and their phylogenetic relationships remain poorly understood. Here, we perform whole-genome sequencing of a gray wolf collected from peninsular India that was phenotypically distinct from gray wolves outside India. Genomic analyses reveal that the Indian gray wolf is an evolutionarily distinct lineage that diverged from other extant gray wolf lineages ∼110 thousand years ago. Demographic analyses suggest that the Indian wolf population declined continuously decline since separating from other gray wolves and, today, has exceptionally low genetic diversity. We also find evidence for pervasive and mosaic gene flow between the Indian wolf and African canids including African wolf, Ethiopian wolf, and African wild dog despite their current geographical separation. Our results support the hypothesis that the Indian subcontinent was a Pleistocene refugium and center of diversification and further highlight the complex history of gene flow that characterized the evolution of gray wolves.

Herpetological phylogeographic analyses support a Miocene focal point of Himalayan uplift and biological diversification.

The Himalaya are among the youngest and highest mountains in the world, but the exact timing of their uplift and origins of their biodiversity are still in debate. The Himalayan region is a relatively small area but with exceptional diversity and endemism. One common hypothesis to explain the rich montane diversity is uplift-driven diversification-that orogeny creates conditions favoring rapid in situ speciation of resident lineages. We test this hypothesis in the Himalayan region using amphibians and reptiles, two environmentally sensitive vertebrate groups. In addition, analysis of diversification of the herpetofauna provides an independent source of information to test competing geological hypotheses of Himalayan orogenesis. We conclude that the origins of the Himalayan herpetofauna date to the early Paleocene, but that diversification of most groups was concentrated in the Miocene. There was an increase in both rates and modes of diversification during the early to middle Miocene, together with regional interchange (dispersal) between the Himalaya and adjacent regions. Our analyses support a recently proposed stepwise geological model of Himalayan uplift beginning in the Paleocene, with a subsequent rapid increase of uplifting during the Miocene, finally giving rise to the intensification of the modern South Asian Monsoon.

Characterization of an Acidic Polysaccharides from Carrot and Its Hepatoprotective Effect on Alcoholic Liver Injury in Mice.

The characteristics of acidic polysaccharides extracted from Daucus carota L. var. sativa Hoffm were investigated and its hepatoprotective effects on alcoholic liver injury were determined in the mice model. A carrot polysaccharide (CPS-I: Carrot polysaccharide-I) with the molecular weight of 3.40×104  kDa was isolated from Daucus carota L. and purified by diethylaminoethyl-52 and Sephadex G-150 column chromatography. The components were analyzed by HPLC, which revealed that CPS-I consisted of galacturonic acid, rhamnose, xylose, arabinose, fructose, and galactose at a relative ratio of 1 : 3.16 : 1.13 : 5.53 : 3.45 : 7.76. Structural characterization analysis suggested that CPS-I was mainly composed of →6)-ß-D-Galp-(1→ and →5)-α-L-Araf-(1→. The hepatoprotective effect of CPS-I was evaluated by alcoholic liver injury mice model. The results showed that the administration of CPS-I (300 mg/kg/day) alleviated the alcoholic liver injury in mice by increasing the levels of ADH and ALDH and reducing oxidative stress. CPS-I ameliorated the pathological changes of liver characterized by lipid accumulation, and reduced the number of lipid droplets.

Multi-locus genetic analyses of Quasipaa from throughout its distribution.

Montane frogs of the genus Quasipaa Dubois, 1992 occur from southern China to Southeast Asia (Frost 2021). Analyses of mtDNA (Cytb) and nuDNA data (Rag1, Rag2, Rhod, Tyr) for samples from 93 localities throughout its distribution yield a phylogeny. Clades A and B occur in Southeast Asia, clade C in northern Yangtze River, China, clade D in southwestern China, and clades E and F in southeastern China. Results place Q. yei within monophyletic Quasipaa and identify two new species. Based on nuDNA data, the basal split of clade A and B indicates an Indochinese origin of Quasipaa. The west-east diversification of five species across South China (Q. spinosa, Q. exilispinosa, Q. jiulongensis, Q. shini, Q. boulengeri) corresponds to topographic terrains II and III of China. Divergence of species from southeastern China (Q. shini, Q. jiulongensis, Q. spinosa, Q. exilispinosa) and southwestern China (Q. boulengeri) dates to 15.30-16.56 Ma (million years ago). A principal component analysis (PCA) and t-test involving 19 bioclimatic variables identifies significantly different environmental conditions between the two regions. Species' distribution models (SDM) for Q. spinosa and Q. boulengeri identify the best areas to be eastern and western South China, respectively. Thus, environmental variation appears to have influenced the genetic divergence and distributions of Quasipaa in South China. Mito-nuclear discordance indicates that some individuals of Q. exilispinosa and Q. spinosa hybridized historically.

Author Correction: 863 genomes reveal the origin and domestication of chicken.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

863 genomes reveal the origin and domestication of chicken.

Despite the substantial role that chickens have played in human societies across the world, both the geographic and temporal origins of their domestication remain controversial. To address this issue, we analyzed 863 genomes from a worldwide sampling of chickens and representatives of all four species of wild jungle fowl and each of the five subspecies of red jungle fowl (RJF). Our study suggests that domestic chickens were initially derived from the RJF subspecies Gallus gallus spadiceus whose present-day distribution is predominantly in southwestern China, northern Thailand and Myanmar. Following their domestication, chickens were translocated across Southeast and South Asia where they interbred locally with both RJF subspecies and other jungle fowl species. In addition, our results show that the White Leghorn chicken breed possesses a mosaic of divergent ancestries inherited from other subspecies of RJF. Despite the strong episodic gene flow from geographically divergent lineages of jungle fowls, our analyses show that domestic chickens undergo genetic adaptations that underlie their unique behavioral, morphological and reproductive traits. Our study provides novel insights into the evolutionary history of domestic chickens and a valuable resource to facilitate ongoing genetic and functional investigations of the world's most numerous domestic animal.

Ancient Hybridization with an Unknown Population Facilitated High-Altitude Adaptation of Canids.

Genetic introgression not only provides material for adaptive evolution but also confounds our understanding of evolutionary history. This is particularly true for canids, a species complex in which genome sequencing and analysis has revealed a complex history of admixture and introgression. Here, we sequence 19 new whole genomes from high-altitude Tibetan and Himalayan wolves and dogs and combine these into a larger data set of 166 whole canid genomes. Using these data, we explore the evolutionary history and adaptation of these and other canid lineages. We find that Tibetan and Himalayan wolves are closely related to each other, and that ∼39% of their nuclear genome is derived from an as-yet-unrecognized wolf-like lineage that is deeply diverged from living Holarctic wolves and dogs. The EPAS1 haplotype, which is present at high frequencies in Tibetan dog breeds and wolves and confers an adaptive advantage to animals living at high altitudes, was probably derived from this ancient lineage. Our study underscores the complexity of canid evolution and demonstrates how admixture and introgression can shape the evolutionary trajectories of species.

New country records for three species of frog from Myanmar including two genera (Nasutixalus and Oreolalax).

Myanmar, a biodiversity hotspot, harbors a striking diversity and endemism of species. Despite this, its herpetofauna remains one of the least explored in continental Asia due to restrictions of crossing political boundaries and infrastructure in remote regions. Many species in adjacent China and India are hypothesized to occur in Myanmar but records are wanting. Recent fieldwork found the frogs Polypedates braueri, Nasutixalus jerdonii and Oreolalax jingdongensis there, and the latter two species represent new generic records for Myanmar. All major morphological characters of these populations match the original descriptions. In addition, our matrilineal genealogy based on DNA barcoding confirms their identities. Overall, these findings confirm that the amphibian diversity is underestimated and this has important implications for conservation. Analyses indicate that northern Myanmar is a biogeographic corridor for the Himalayas, southern China, and northeastern India.

Synthesis of selenylated polysaccharides from Momordica charantia L. and its hypoglycemic activity in streptozotocin-induced diabetic mice.

In order to obtain the effective antidiabetic polysaccharide derivative, a selenylated polysaccharide (Se-MCPIIa-1), with an average molecular weight (MW) of 4.0038 × 104 Da, was synthesized by reduction of sodium selenite with ascorbic acid in the presence of Momordica polysaccharides (MCPIIa). The selenium (Se) content of Se-MCPIIa-1 was up to 445.0 µg/g, and its average diameter of monodisperse spherical particle size was around 63.78 nm. The morphology and physicochemical properties of Se-MCPIIa-1 were characterized by scanning electron microscope (SEM), atomic force microscope (AFM), Fourier transform infrared spectrometry (FT-IR) and Raman spectroscopy, respectively. The results indicated that Se was conjunct with MCPIIa by esterification. Moreover, the oral administration of Se-MCPIIa-1 showed a gradual normalization in the levels of hypoglyemic test in the STZ-induced diabetic mice. The anti-diabetic effects of Se-MCPIIa-1 in vivo showed that Se-MCPIIa-1 can significantly reduce fasting blood glucose levels and enhance insulin levels as well as antioxidant enzyme activities in diabetic mice with an optimal dosage of 20 mg/kg/body weight. In addition, it was found from histopathological data that Se-MCPIIa-1 could prevent pancreatic islets, liver and kidney damage from diabetes, which suggested that Se-MCPIIa-1 is a promising novel Se supplement and could be applied in the field of food and medicine.

A combined approach of mitochondrial DNA and anchored nuclear phylogenomics sheds light on unrecognized diversity, phylogeny, and historical biogeography of the torrent frogs, genus Amolops (Anura: Ranidae).

The genus Amolops ("torrent frogs") is one of the most species-rich genera in Ranidae, with 59 recognized species. This genus currently includes six species groups diagnosed mainly by morphology. Several recent molecular studies indicated that the classification of species groups within Amolops remains controversial, and key nodes in the phylogeny have been inadequately resolved. In addition, the diversity of Amolops remains poorly understood, especially for those from incompletely sampled regions. Herein, we investigate species-level diversity within the genus Amolops throughout southern China and Southeast Asia, and infer evolutionary relationships among the species using mtDNA data (16S, COI, and ND2). Molecular analyses indicate nine unnamed species, mostly distributed in the Himalayas. We then utilized anchored hybrid enrichment to generate a dataset representing the major mitochondrial lineages to resolve phylogenetic relationships, biogeography, and pattern of species diversification. Our resulting phylogeny strongly supports the monophyly of four previously identified species groups (the A. ricketti, A. daiyunensis, A. hainanensis, and A. monticola groups), but paraphyly for the A. mantzorum and A. marmoratus groups, as previously defined. We erect one new species group, the A. viridimaculatus group, and recognize Dubois' (1992) subgenus Amo as the A. larutensis species group. Biogeographic analysis suggests that Amolops originated on the Indo-Burma/Thai-Malay Peninsula at the Eocene/Oligocene boundary, and dispersed outward, exemplifying a common pattern observed for the origin of Asian biodiversity. The early divergence within Amolops coincides with the Himalayan uplift and the lateral extrusion of Indochina at the Oligocene/Miocene boundary. Our results show that paleoclimatic and geomorphological events have profoundly influenced the patterns of lineage diversification within Amolops.

An integrative phylogenomic approach illuminates the evolutionary history of Old World tree frogs (Anura: Rhacophoridae).

Rhacophoridae are one of the most speciose and ecologically diverse families of amphibians. Resolution of their evolutionary relationships is key to understanding the accumulation of biodiversity, yet previous hypotheses based on Sanger sequencing exhibit much discordance amongst generic relationships. This conflict precludes the making of sound macroevolutionary conclusions. Herein, we conduct the first phylogenomic study using broad-scale sampling and sequences of 352 nuclear DNA loci obtained using anchored hybrid enrichment targeted sequencing. The robust time-calibrated phylogenetic hypothesis clarifies several long-disputed relationships and facilitates the testing of evolutionary hypotheses on spatiotemporal diversification and reproductive modes. The major extant lineages of Rhacophoridae appear to have radiated in mainland Asia, and the spatiotemporal process corresponds with several common accumulations of biodiversity in Asia. Analyses do not detect any case of "Out of Himalaya" in Rhacophoridae. All transitions of reproductive modes appear to have evolved in an ordered, gradual sequence associated with gaining independence of standing water for larval development. The different reproductive modes are phylogenetically conserved and the completion of their transitions appear to have occurred over a period of ~30 Ma, which does not fit a pattern of a rapid burst of diversification. Innovations in reproductive modes associate statistically with the uneven distribution of species-richness between clades, where higher diversification is linked to increased terrestrial modes of reproduction. These results strengthen the hypothesis that breeding innovations drive diversification by providing new opportunities for ecological release and dispersion.

A new species of the genus Raorchestes (Anura: Rhacophoridae) from Yunnan Province, China.

A new bush frog species, Raorchestes cangyuanensis sp. nov., from Cangyuan, Yunnan Province, China, is described based on morphological and molecular analyses. It differs from all known congeners by a combination of the following characters: body size small, adult snout-vent length (SVL) 16.1-20.0 mm in males (n=3); tympanum indistinct; tips of all fingers and toes expanded into discs with circummarginal grooves; rudimentary webbing between toes; fingers and toes with lateral dermal fringes; inner and outer metacarpal tubercles present; heels meeting when limbs held at right angles to body; crotch with a distinct black patch; discs of fingers and toes orange; male with external single subgular vocal sac and reddish nuptial pad at the base of first finger.

Genomic and transcriptomic investigations of the evolutionary transition from oviparity to viviparity.

Viviparous (live-bearing) vertebrates have evolved repeatedly within otherwise oviparous (egg-laying) clades. Over two-thirds of these changes in vertebrate reproductive parity mode happened in squamate reptiles, where the transition has happened between 98 and 129 times. The transition from oviparity to viviparity requires numerous physiological, morphological, and immunological changes to the female reproductive tract, including eggshell reduction, delayed oviposition, placental development for supply of water and nutrition to the embryo by the mother, enhanced gas exchange, and suppression of maternal immune rejection of the embryo. We performed genomic and transcriptomic analyses of a closely related oviparous-viviparous pair of lizards (Phrynocephalus przewalskii and Phrynocephalus vlangalii) to examine these transitions. Expression patterns of maternal oviduct through reproductive development of the egg and embryo differ markedly between the two species. We found changes in expression patterns of appropriate genes that account for each of the major aspects of the oviparity to viviparity transition. In addition, we compared the gene sequences in transcriptomes of four oviparous-viviparous pairs of lizards in different genera (Phrynocephalus, Eremias, Scincella, and Sphenomorphus) to look for possible gene convergence at the sequence level. We discovered low levels of convergence in both amino acid replacement and evolutionary rate shift. This suggests that most of the changes that produce the oviparity-viviparity transition are changes in gene expression, so occasional reversals to oviparity from viviparity may not be as difficult to achieve as has been previously suggested.

Preparation of a Momordica charantia L. polysaccharide­chromium (III) complex and its anti-hyperglycemic activity in mice with streptozotocin-induced diabetes.

Polysaccharides comprise the major bioactive components in Momordica charantia L. We here synthesized and characterized a novel M. charantia polysaccharide­chromium (III) complex (MCPIIaC) and assessed its anti-diabetic effects in mice with streptozotocin (STZ)-induced diabetes mellitus (DM) and its mechanism underlying hypoglycemia. MCPIIaC is a novel polysaccharide­chromium (III) complex containing 14.68% elemental chromium. A Fourier transform infrared (FTIR) spectrogram experiment showed that the chromium ions are linked to the polysaccharide's hydroxyl groups. Combined circular dichroism (CD) and atomic force microscopy (AFM) analyses indicated that after linking with chromium ions, the flexibility of the three-dimensional structure of the polysaccharide increased. After the treatment of MCPIIaC for four weeks, the mice with STZ-induced DM exhibited significantly lower fasting blood glucose levels and body weight, whereas higher insulin levels and antioxidant enzyme activity than in the diabetic group. Optimal effects were obtained with a dosage of 30 mg MCPIIaC/kg body weight. Histological analysis indicated that MCPIIaC alleviated the oxidative tissue damage in STZ-lesioned mice. An acute toxicity experiment indicated that MCPIIaC was safe at a dose of 1500 mg/kg. These results suggest that MCPIIaC might be an excellent candidate hypoglycemic agent for the prevention of diabetes.

Transcriptomic analysis of Momordica charantia polysaccharide on streptozotocin-induced diabetic rats.

Polysaccharides obtained from Momordica charantia have been shown to exert a hypoglycemic effect, thereby improving glucose metabolism. However, little is known regarding the molecular mechanisms underlying this process. In recent years, RNA sequencing (RNA-seq) has rapidly emerged as a major transcriptome profiling system. Herein, RNA-seq was utilized to obtain whole transcriptomes from the livers of SD rats with type 2 diabetes (DM) and type 2 diabetic rats treated with M. charantia polysaccharides (DMCPIIa). In total, 44,752,508-52,585,920 uniquely mapped reads were obtained, covering 85.54% of the annotated transcripts and representing 12,857 mRNA transcripts. Following analytical analysis, 17 differentially expressed genes (p < 0.05, false discovery rate q < 0.05) were identified in the DMCPIIa group relative to the DM group. Gene ontology and pathway analysis demonstrated that 17 differently expressed genes were enriched in specific biological processes mainly relating to glucose metabolism and fat metabolism (p < 0.05). Furthermore, a subset of the differential genes was selected for validation via real-time quantitative reverse-transcription PCR. Collectively, the integrated analysis of differential gene expression, data obtained herein and the examination of previously reported quantitative trait loci, and genome-wide association studies has suggested that pdk4, pkl, hsd11ß1, msmo1, rbp4 and fads2 may serve as promising candidates for the modulation of type 2 diabetes.

Selection and environmental adaptation along a path to speciation in the Tibetan frog Nanorana parkeri.

Tibetan frogs, Nanorana parkeri, are differentiated genetically but not morphologically along geographical and elevational gradients in a challenging environment, presenting a unique opportunity to investigate processes leading to speciation. Analyses of whole genomes of 63 frogs reveal population structuring and historical demography, characterized by highly restricted gene flow in a narrow geographic zone lying between matrilines West (W) and East (E). A population found only along a single tributary of the Yalu Zangbu River has the mitogenome only of E, whereas nuclear genes of W comprise 89-95% of the nuclear genome. Selection accounts for 579 broadly scattered, highly divergent regions (HDRs) of the genome, which involve 365 genes. These genes fall into 51 gene ontology (GO) functional classes, 14 of which are likely to be important in driving reproductive isolation. GO enrichment analyses of E reveal many overrepresented functional categories associated with adaptation to high elevations, including blood circulation, response to hypoxia, and UV radiation. Four genes, including DNAJC8 in the brain, TNNC1 and ADORA1 in the heart, and LAMB3 in the lung, differ in levels of expression between low- and high-elevation populations. High-altitude adaptation plays an important role in maintaining and driving continuing divergence and reproductive isolation. Use of total genomes enabled recognition of selection and adaptation in and between populations, as well as documentation of evolution along a stepped cline toward speciation.