Large-scale snake genome analyses provide insights into vertebrate development.

Snakes are a remarkable squamate lineage with unique morphological adaptations, especially those related to the evolution of vertebrate skeletons, organs, and sensory systems. To clarify the genetic underpinnings of snake phenotypes, we assembled and analyzed 14 de novo genomes from 12 snake families. We also investigated the genetic basis of the morphological characteristics of snakes using functional experiments. We identified genes, regulatory elements, and structural variations that have potentially contributed to the evolution of limb loss, an elongated body plan, asymmetrical lungs, sensory systems, and digestive adaptations in snakes. We identified some of the genes and regulatory elements that might have shaped the evolution of vision, the skeletal system and diet in blind snakes, and thermoreception in infrared-sensitive snakes. Our study provides insights into the evolution and development of snakes and vertebrates.

Genetic Basis and Evolutionary Forces of Sexually Dimorphic Color Variation in a Toad-Headed Agamid Lizard.

In the animal kingdom, sexually dimorphic color variation is a widespread phenomenon that significantly influences survival and reproductive success. However, the genetic underpinnings of this variation remain inadequately understood. Our investigation into sexually dimorphic color variation in the desert-dwelling Guinan population of the toad-headed agamid lizard (Phrynocephalus putjatai) utilized a multidisciplinary approach, encompassing phenotypic, ultrastructural, biochemical, genomic analyses, and behavioral experiments. Our findings unveil the association between distinct skin colorations and varying levels of carotenoid and pteridine pigments. The red coloration in males is determined by a genomic region on chromosome 14, housing four pigmentation genes: BCO2 and three 6-pyruvoyltetrahydropterin synthases. A Guinan population-specific nonsynonymous single nucleotide polymorphism in BCO2 is predicted to alter the electrostatic potential within the binding domain of the BCO2-ß-carotene complex, influencing their interaction. Additionally, the gene MAP7 on chromosome 2 emerges as a potential contributor to the blue coloration in subadults and adult females. Sex-specific expression patterns point to steroid hormone-associated genes (SULT2B1 and SRD5A2) as potential upstream regulators influencing sexually dimorphic coloration. Visual modeling and field experiments support the potential selective advantages of vibrant coloration in desert environments. This implies that natural selection, potentially coupled with assortative mating, might have played a role in fixing color alleles, contributing to prevalence in the local desert habitat. This study provides novel insights into the genetic basis of carotenoid and pteridine-based color variation, shedding light on the evolution of sexually dimorphic coloration in animals. Moreover, it advances our understanding of the driving forces behind such intricate coloration patterns.

Gatorbulin-1, a distinct cyclodepsipeptide chemotype, targets a seventh tubulin pharmacological site.

Tubulin-targeted chemotherapy has proven to be a successful and wide spectrum strategy against solid and liquid malignancies. Therefore, new ways to modulate this essential protein could lead to new antitumoral pharmacological approaches. Currently known tubulin agents bind to six distinct sites at α/ß-tubulin either promoting microtubule stabilization or depolymerization. We have discovered a seventh binding site at the tubulin intradimer interface where a novel microtubule-destabilizing cyclodepsipeptide, termed gatorbulin-1 (GB1), binds. GB1 has a unique chemotype produced by a marine cyanobacterium. We have elucidated this dual, chemical and mechanistic, novelty through multidimensional characterization, starting with bioactivity-guided natural product isolation and multinuclei NMR-based structure determination, revealing the modified pentapeptide with a functionally critical hydroxamate group; and validation by total synthesis. We have investigated the pharmacology using isogenic cancer cell screening, cellular profiling, and complementary phenotypic assays, and unveiled the underlying molecular mechanism by in vitro biochemical studies and high-resolution structural determination of the α/ß-tubulin-GB1 complex.

Bavachin protects against diet-induced hepatic steatosis and obesity in mice.

Non-alcoholic fatty liver disease (NAFLD) is a major health concern worldwide, and the incidence of metabolic disorders associated with NAFLD is rapidly increasing because of the obesity epidemic. There are currently no approved drugs that prevent or treat NAFLD. Recent evidence shows that bavachin, a flavonoid isolated from the seeds and fruits of Psoralea corylifolia L., increases the transcriptional activity of PPARγ and insulin sensitivity during preadipocyte differentiation, but the effect of bavachin on glucose and lipid metabolism remains unclear. In the current study we investigated the effects of bavachin on obesity-associated NAFLD in vivo and in vitro. In mouse primary hepatocytes and Huh7 cells, treatment with bavachin (20 µM) significantly suppressed PA/OA or high glucose/high insulin-induced increases in the expression of fatty acid synthesis-related genes and the number and size of lipid droplets. Furthermore, bavachin treatment markedly elevated the phosphorylation levels of AKT and GSK-3ß, improving the insulin signaling activity in the cells. In HFD-induced obese mice, administration of bavachin (30 mg/kg, i.p. every other day for 8 weeks) efficiently attenuated the increases in body weight, liver weight, blood glucose, and liver and serum triglyceride contents. Moreover, bavachin administration significantly alleviated hepatic inflammation and ameliorated HFD-induced glucose intolerance and insulin resistance. We demonstrated that bavachin protected against HFD-induced obesity by inducing fat thermogenesis and browning subcutaneous white adipose tissue (subWAT). We revealed that bavachin repressed the expression of lipid synthesis genes in the liver of obese mice, while promoting the expression of thermogenesis, browning, and mitochondrial respiration-related genes in subWAT and brown adipose tissue (BAT) in the mice. In conclusion, bavachin attenuates hepatic steatosis and obesity by repressing de novo lipogenesis, inducing fat thermogenesis and browning subWAT, suggesting that bavachin is a potential drug for NAFLD therapy.

Discovery, Synthesis, and Biological Evaluation of Anaenamides C and D from a New Marine Cyanobacterium, Hormoscilla sp.

Our ongoing efforts to explore the chemical space associated with marine cyanobacteria from coral reefs of Guam have yielded two new members of the anaenamide family of natural products, anaenamides C (3) and D (4). These compounds were isolated from a novel Hormoscilla sp. (VPG16-58). Our phylogenetic profiling (16S rDNA) of this cyanobacterium indicated that VPG16-58 is taxonomically distinct from the previously reported producer of the anaephenes, VPG16-59 (Hormoscilla sp.), and other previously documented species of the genus Hormoscilla. The planar structures of 3 and 4 were determined via spectroscopic methods, and absolute configurations of the α-hydroxy acids were assigned by enantioselective HPLC analysis. To address the requirement for sufficient material for testing, we first adapted our published linear synthetic approach for 1 and 2 to generate anaenoic acid (7), which served as a point for diversification, providing the primary amides 3 and 4 from synthetic intermediates 5 and 6, respectively. The compounds were then tested for effects on HCT116 colon cancer cell viability and in an ARE-luciferase reporter gene assay for Nrf2 modulation using HEK293 human embryonic kidney cells. Our findings indicate that, in contrast to cytotoxic methyl esters 1 and 2, the primary amides 3 and 4 activate the Nrf2 pathway at noncytotoxic concentrations. Overall, our data suggest that the anaenamide scaffold is tunable to produce differential biological outcomes.

Bifunctional Doscadenamides Activate Quorum Sensing in Gram-Negative Bacteria and Synergize with TRAIL to Induce Apoptosis in Cancer Cells.

New cyanobacteria-derived bifunctional analogues of doscadenamide A, a LasR-dependent quorum sensing (QS) activator in Pseudomonas aeruginosa, characterized by dual acylation of the pyrrolinone core structure and the pendant side chain primary amine to form an imide/amide hybrid are reported. The identities of doscadenamides B-J were confirmed through total synthesis and a strategic focused library with different acylation and unsaturation patterns was created. Key molecular interactions for binding with LasR and a functional response through mutation studies coupled with molecular docking were identified. The structure-activity relationships (SARs) were probed in various Gram-negative bacteria, including P. aeruginosa and Vibrio harveyi, indicating that the pyrrolinone-N acyl chain is critical for full agonist activity, while the other acyl chain is dispensable or can result in antagonist activity, depending on the bacterial system. Since homoserine lactone (HSL) quorum sensing activators have been shown to act in synergy with TRAIL to induce apoptosis in cancer cells, selected doscadenamides were tested in orthogonal eukaryotic screening systems. The most potent QS agonists, doscadenamides S10-S12, along with doscadenamides F and S4 with partial or complete saturation of the acyl side chains, exhibited the most pronounced synergistic effects with TRAIL in triple negative MDA-MB-231 breast cancer cells. The overall correlation of the SAR with respect to prokaryotic and eukaryotic targets may hint at coevolutionary processes and intriguing host-bacteria relationships. The doscadenamide scaffold represents a non-HSL template for combination therapy with TRAIL pathway stimulators.

Largazole Inhibits Ocular Angiogenesis by Modulating the Expression of VEGFR2 and p21.

Ocular angiogenic diseases, characterized by abnormal blood vessel formation in the eye, are the leading cause of blindness. Although Anti-VEGF therapy is the first-line treatment in the market, a substantial number of patients are refractory to it or may develop resistance over time. As uncontrolled proliferation of vascular endothelial cells is one of the characteristic features of pathological neovascularization, we aimed to investigate the role of the class I histone deacetylase (HDAC) inhibitor Largazole, a cyclodepsipeptide from a marine cyanobacterium, in ocular angiogenesis. Our study showed that Largazole strongly inhibits retinal vascular endothelial cell viability, proliferation, and the ability to form tube-like structures. Largazole strongly inhibits the vessel outgrowth from choroidal explants in choroid sprouting assay while it does not affect the quiescent choroidal vasculature. Largazole also inhibits vessel outgrowth from metatarsal bones in metatarsal sprouting assay without affecting pericytes coverage. We further demonstrated a cooperative effect between Largazole and an approved anti-VEGF drug, Alflibercept. Mechanistically, Largazole strongly inhibits the expression of VEGFR2 and leads to an increased expression of cell cycle inhibitor, p21. Taken together, our study provides compelling evidence on the anti-angiogenic role of Largazole that exerts its function through mediating different signaling pathways.

A Rare Case of Hb H Disease and Systemic Lupus Erythematosus.

Anemia is common in patients with systemic lupus erythematosus (SLE). The association between thalassemia and SLE is rare. In this study, we report the first patient who was found to have a severe hemolytic anemia caused by combination of SLE and Hb H disease. The patient had a more severe presentation in the hematological system. Our case indicates that for a patient who was diagnosed with SLE and developed deterioration in her hematological cell lines, investigation of other possible coexisting causes would be warranted.

Ahp-Cyclodepsipeptides as tunable inhibitors of human neutrophil elastase and kallikrein 7: Total synthesis of tutuilamide A, serine protease selectivity profile and comparison with lyngbyastatin 7.

We describe the total synthesis of tutuilamide A, a potent porcine pancreatic elastase (PPE) inhibitor and a representative member of the 3-amino-6-hydroxy-2-piperidone (Ahp) cyclodepsipeptide family, isolated from marine cyanobacteria. The Ahp unit serves as a pharmacophore and the adjacent 2-amino-2-butenoic acid (Abu) is a main driver of the selectivity among serine proteases. We adapted our previous convergent strategy to generate the macrocycle, common with lyngbyastatin 7 and related elastase inhibitors, and then appended the tutuilamide A-specific side chain bearing a vinyl chloride. Tutuilamide A and lyngbyastatin 7 were evaluated side by side for the inhibition of the disease-relevant human neutrophil elastase (HNE). Tutuilamide A and lyngbyastatin 7 were approximately equipotent against HNE, while tutuilamide A was previously shown to be more active against PPE compared with lyngbyastatin 7, further demonstrating that the side chain provides opportunities to not only modulate potency but also selectivity among proteases of the same function from different organisms. Profiling of tutuilamide A against mainly human serine proteases revealed high selectivity for HNE (IC50 0.73 nM) and pleiotropic activity against kallikrein 7 (KLK7, IC50 5.0 nM), without affecting other kallikreins, similarly to lyngbyastatin 7 (IC50 0.85 nM for HNE and 3.1 nM for KLK7). A comprehensive molecular docking study for elastases and KLK7 afforded deeper insight into the intricate differences between inhibitor interactions with HNE and PPE, accounting for the differential activities for both compounds. The synthesis and molecular studies serve as a proof-of-concept that the macrocyclic scaffold can be diversified to fine-tune the activity of serine protease inhibitors.

Development of apratoxin S10 (Apra S10) as an anti-pancreatic cancer agent and its preliminary evaluation in an orthotopic patient-derived xenograft (PDX) model.

Despite the significant progress in the field of cancer therapeutics, the incidence of pancreatic cancer (PC) has continuously increased. One possible mechanism for this increasing burden is impaired drug delivery and drug resistance resulting from a unique tumor microenvironment and genetic mutations. Apratoxins are potent anticancer agents and cotranslational translocation inhibitors with potential therapeutic applications to treat cancers with active secretory pathways. Here, we developed apratoxin S10 (Apra S10) as an anti-pancreatic cancer agent which potently inhibited the growth of both established and patient-derived primary pancreatic cancer cells. We validated its mechanism of action on pancreatic cancer cells by demonstrating the downregulation of multiple receptor tyrosine kinases and inhibition of growth factor and cytokine secretion. Apra S10 also inhibited a number of cytokines secreted by stromal cells, suggesting that Apra S10 not only inhibited pancreatic cancer cell secretion, but also reduced the level of factors secreted by other cell types active within the tumor microenvironment. As Apra S10 tissue distribution indicated its high enrichment in pancreas tissue, an orthotopic pancreatic patient-derived xenograft mouse model that closely mimics the human pancreatic tumor microenvironment was for the first time used in apratoxin studies. Apra S10 showed promising antitumor effect in this pancreatic cancer model and this effect was mediated through anti-proliferation properties.

Evolutionary melting pots and reproductive isolation: A ring-shaped diversification of an odorous frog (Odorrana margaratea) around the Sichuan Basin.

Quaternary climatic oscillations and geography are of primary importance in shaping intraspecific genetic diversity. We examined the diversification patterns and inferred processes for the green odorous frog (Odorrana margaretae) of western China. Species distribution modelling showed that the species has a continuous circular distribution around the Sichuan Basin while the basin itself is largely uninhabitable. Population genetic and phylogenetic analyses revealed that the species has a ring-shaped divergent pattern. While the chain of populations around the Basin maintains a mostly gradual and continuous genetic variation, populations between the north and west showed little gene exchange. Two processes, glacial refugial history and geography, likely contributed to the observed patterns. Our genetic clustering analysis revealed two clusters, suggesting two refugial groups among the populations, one from the west and the other from the east. Postglacial expansion may have created two contact zones. One at the south had extensive population admixture and produced a gradual transition between the western and eastern populations. Consequently, this region has the highest genetic diversity and represents an evolutionary "melting pot." In contrast, the second contact zone at the northwestern side of the Basin has limited admixture, suggesting partial reproductive isolation between the northern and western populations. Furthermore, an isolation-by-distance analysis revealed a strong correlation (Mantel r = 0.745) between the genetic and geographic distances, when we assumed that populations were connected following the circular distribution without gene flow across the NW contact zone. We also explored alternative explanations, such as a one-refugium scenario. With its micro-ring, the green odorous frog is poised to make an excellent model system for understanding the process of speciation.

Testing domain general learning in an Australian lizard.

A key question in cognition is whether animals that are proficient in a specific cognitive domain (domain specific hypothesis), such as spatial learning, are also proficient in other domains (domain general hypothesis) or whether there is a trade-off. Studies testing among these hypotheses are biased towards mammals and birds. To understand constraints on the evolution of cognition more generally, we need broader taxonomic and phylogenetic coverage. We used Australian eastern water skinks (Eulamprus quoyii) with known spatial learning ability in three additional tasks: an instrumental and two discrimination tasks. Under domain specific learning we predicted that lizards that were good at spatial learning would perform less well in the discrimination tasks. Conversely, we predicted that lizards that did not meet our criterion for spatial learning would likewise perform better in discrimination tasks. Lizards with domain general learning should perform approximately equally well (or poorly) in these tasks. Lizards classified as spatial learners performed no differently to non-spatial learners in both the instrumental and discrimination learning tasks. Nevertheless, lizards were proficient in all tasks. Our results reveal two patterns: domain general learning in spatial learners and domain specific learning in non-spatial learners. We suggest that delineating learning into domain general and domain specific may be overly simplistic and we need to instead focus on individual variation in learning ability, which ultimately, is likely to play a key role in fitness. These results, in combination with previously published work on this species, suggests that this species has behavioral flexibility because they are competent across multiple cognitive domains and are capable of reversal learning.

Discovery of new A- and B-type laxaphycins with synergistic anticancer activity.

Two new cyclic lipopeptides termed laxaphycins B4 (1) and A2 (2) were discovered from a collection of the marine cyanobacterium Hormothamnion enteromorphoides, along with the known compound laxaphycin A. The planar structures were solved based on a combined interpretation of 1D and 2D NMR data and mass spectral data. The absolute configurations of the subunits were determined by chiral LC-MS analysis of the hydrolysates, advanced Marfey's analysis and 1D and 2D ROESY experiments. Consistent with similar findings on other laxaphycin A- and B-type peptides, laxaphycin B4 (1) showed antiproliferative effects against human colon cancer HCT116 cells with IC50 of 1.7 µM, while laxaphycins A and A2 (2) exhibited weak activities. The two major compounds isolated from the sample, laxaphycins A and B4, were shown to act synergistically to inhibit the growth of HCT116 colorectal cancer cells.

Gene expression variations in high-altitude adaptation: a case study of the Asiatic toad (Bufo gargarizans).

BACKGROUND: Genome-wide investigation of molecular mechanisms for high-altitude adaptation has attracted great attention in the last few years. In order to understand the contribution of gene expression level variations to high-altitude adaptation in Asiatic toads (Bufo gargarizans), we implemented a reciprocal transplant experiment between low- and high-altitude sites and sequenced 12 transcriptomes from brain, heart, and liver tissues. RESULTS: A large number of genes with expression differences (DEGs) between high- and low-altitude individuals (193 fixed and 844 plastic) were identified, and the majority of them were tissue specific. Heart displayed the largest number of DEGs, both plastic and fixed. Fixed DEGs were particularly concentrated in functions associated with muscle contraction, and the majority of them were down-regulated in high-altitude individuals. Plastic DEGs were highly concentrated in several energy metabolism related functional categories, and the majority of them were also down-regulated at high-altitude environments. In liver samples, genes associated with nutrient metabolism experienced a broad-scale expression down-regulation in high-altitude toads. CONCLUSIONS: These broadly suppressed expression patterns at high altitudes are in strong contrast to those of endothermic homeotherms, suggesting poikilothermic vertebrates may have adopted different strategies at high altitudes. Our results strongly support that both genotypic specialization and phenotypic plasticity play crucial role in adaptation to high altitude for Asiatic toads. Poikilothermic vertebrates are among the most hypoxia-tolerant animals known, and many molecular mechanisms remain elusive. We hope that our results will provide useful directions for future research.

Genetic signals of high-altitude adaptation in amphibians: a comparative transcriptome analysis.

BACKGROUND: High-altitude adaptation provides an excellent system for studying how organisms cope with multiple environmental stressors and interacting genetic modifications. To explore the genetic basis of high-altitude adaptation in poikilothermic animals, we acquired transcriptome sequences from a high-altitude population and a low-altitude population of the Asiatic toad (Bufo gargarizans). Transcriptome data from another high-altitude amphibian, Rana kukunoris and its low-altitude relative R. chensiensis, which are from a previous study, were also incorporated into our comparative analysis. RESULTS: More than 40,000 transcripts were obtained from each transcriptome, and 5107 one-to-one orthologs were identified among the four taxa for comparative analysis. A total of 29 (Bufo) and 33 (Rana) putative positively selected genes were identified for the two high-altitude species, which were mainly concentrated in nutrient metabolism related functions. Using SNP-tagging and FST outlier analysis, we further tested 89 other nutrient metabolism related genes for signatures of natural selection, and found that two genes, CAPN2 and ITPR1, were likely under balancing selection. We did not detect any positively selected genes associated with response to hypoxia. CONCLUSIONS: Amphibians clearly employ different genetic mechanisms for high-altitude adaptation compared to endotherms. Modifications of genes associated with nutrient metabolism feature prominently while genes related to hypoxia tolerance appear to be insignificant. Poikilotherms represent the majority of animal diversity, and we hope that our results will provide useful directions for future studies of amphibians as well as other poikilotherms.

Total Synthesis of the Potent Marine-Derived Elastase Inhibitor Lyngbyastatin 7 and in Vitro Biological Evaluation in Model Systems for Pulmonary Diseases.

Lyngbyastatin 7 (1) is a marine cyanobacteria-derived lariat-type cyclic depsipeptide of which the macrocyclic core possesses modified amino acids, including a featured 3-amino-6-hydroxy-2-piperidone (Ahp) moiety and a (Z)-2-amino-2-butenoic acid (Abu) moiety. The first total synthesis of 1 was successfully established via 31 steps, and the conditions of several crucial steps were optimized to ensure smooth operations. The previously reported structural assignment and elastase inhibitory activity of the isolated natural product were confirmed. According to the extensive in vitro biological evaluation, compound 1 displayed low nanomolar IC50 in blocking elastase activity and strong ability in protecting bronchial epithelial cells against elastase-induced antiproliferation and abrogating the elastase-triggered induction of pro-inflammatory cytokine expression. Its overall performance was superior over sivelestat, the only approved small molecule drug targeting elastase, which indicated its potential in developing as a pharmacotherapeutic against elastase-mediated pathologies. The success in total synthesis, designed with a novel convergent strategy, not only overcame the supply issue for thorough preclinical studies but also paved the way for convenient synthesis of analogues with improved potency and druglike properties.

Application of Serratia marcescens RZ-21 significantly enhances peanut yield and remediates continuously cropped peanut soil.

BACKGROUND: Continuous cropping practices cause a severe decline in peanut yield. The aim of this study was to investigate the remediation effect of Serratia marcescens on continuously cropped peanut soil. A pot experiment was conducted under natural conditions to determine peanut agronomic indices, soil microorganism characteristics, soil enzyme activities and antagonism ability to typical pathogens at different growth stages. Four treatments were applied to red soil as follows: an active fermentation liquor of S. marcescens (RZ-21), an equivalent sterilized fermentation liquor (M), an equivalent fermentation medium (P) and distilled water (CK). RESULTS: S. marcescens significantly inhibited the two typical plant pathogens Fusarium oxysporum A1 and Ralstonia solanacearum B1 and reduced their populations in rhizosphere soil. The RZ-21 treatment significantly increased peanut yield, vine dry weight, root nodules and taproot length by 62.3, 33, 72 and 61.4% respectively, followed by the M treatment. The P treatment also increased root nodules and root length slightly. RZ-21 also enhanced the activities of soil urease, sucrase and hydrogen peroxidase at various stages. In addition, RZ-21 and M treatments increased the average population of soil bacteria and decreased the average population of fungi in the three critical peanut growth stages, except for M in the case of the fungal population at flowering, thus balancing the structure of the soil microorganism community. CONCLUSION: This is the first report of S. marcescens being applied to continuously cropped peanut soil. The results suggest that S. marcescens RZ-21 has the potential to improve the soil environment and agricultural products and thus allow the development of sustainable management practices.

Randomized, controlled pharmacokinetic and pharmacodynamic evaluation of albinterferon in patients with chronic hepatitis B infection.

BACKGROUND AND AIMS: Albinterferon is a fusion of albumin and interferon-α2b developed to improve the pharmacokinetics, convenience, and potential efficacy of interferon-α for the treatment of chronic hepatitis infections. METHODS: This open-label, randomized, active-controlled, multicenter study investigated the safety and efficacy of albinterferon in patients with chronic hepatitis B virus (HBV) infection who were e-antigen (HBeAg) positive. One hundred and forty-one patients received one of four albinterferon doses/regimens or pegylated-interferon-α2a. Primary efficacy outcomes were changes in serum HBeAg and antibody, HBV-DNA, and alanine aminotransferase. Principal safety outcomes were changes in laboratory values, pulmonary function, and adverse events. RESULTS: The study was prematurely terminated as phase III trials in hepatitis C infection indicated noninferior efficacy but inferior safety compared with pegylated-interferon-α2a. Here, all treatment groups had a significant reduction in HBV-DNA from baseline. Reductions in HBV-DNA were not significantly different, except the 1200 µg every 4 weeks albinterferon dose which was inferior compared with pegylated-interferon-α2a. The serum alanine aminotransferase levels decreased in all arms. The per-patient incidence of adverse events was not significantly different for albinterferon (96.4-100%) and pegylated-interferon-α2a (93.1%). Total adverse events, however, were higher for albinterferon and correlated to dose. Decreased lung function was found in all arms (∼93% of patients), and was more common in some albinterferon groups. CONCLUSIONS: Albinterferon doses with similar anti-HBV efficacy to pegylated-interferon-α2a had higher rates of certain adverse events, particularly changes in lung diffusion capacity (http://www.clinicaltrials.gov number NCT00964665).

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris.

Mitochondrial genome (mito-genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito-genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex-biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt-b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt-b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial-nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt-b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito-nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex-biased dispersal. A weak pattern of female-biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex-biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris-type mito-genome in thermal adaptation may also contribute to the introgression between the two species.