Dynamic Profiling of Protein Mole Synthesis Rates during C2C12 Myoblast Differentiation

Mole (MSR) and fractional (FSR) synthesis rates of proteins during C2C12 myoblast differentiation are investigated. Myoblast cultures supplemented with D₂O during 0–24 h or 72–96 h of differentiation are analyzed by LC‐MS/MS to calculate protein FSR and MSR after samples are spiked with yeast alcohol dehydrogenase (ADH1). Profiling of 153 proteins detected 70 significant (p ≤ 0.05, FDR ≤ 1%) differences in abundance between cell states. Early differentiation is enriched by clusters of ribosomal and heat shock proteins, whereas later differentiation is associated with actin filament binding. The median (first–third quartile) FSR (%/h) during early differentiation 4.1 (2.7–5.3) is approximately twofold greater than later differentiation 1.7 (1.0–2.2), equating to MSR of 0.64 (0.38–1.2) and 0.28 (0.1–0.5) fmol h^?1 µg^?1 total protein, respectively. MSR corresponds more closely with abundance data and highlights proteins associated with glycolytic processes and intermediate filament protein binding that are not evident among FSR data. Similarly, MSR during early differentiation accounts for 78% of the variation in protein abundance during later differentiation, whereas FSR accounts for 4%. Conclusively, the interpretation of protein synthesis data differs when reported in mole or fractional terms, which has consequences when studying the allocation of cellular resources.     

Horizontal, but not vertical canopy structure is related to stand functional diversity in a subtropical slope forest

The aim of this study was to analyse the relation of horizontal and vertical canopy structure to tree functional diversity of a highly diverse subtropical broad-leaved slope forest, stratified for different successional stages. This is of particular interest because many key ecosystem processes and functions are related to the arrangement of forest canopies. We assessed the effect of stand-related functional diversity (FD_Q, measured as Rao’s quadratic entropy of leaf traits), together with other environmental variables on horizontal [measured as relative crown projection areas (CPA_r)] and vertical [relative crown overlap, coefficients of variation (CV) of crown positioning variables] structure of the upper canopy at the local neighbourhood level. The analyses with mixed effects models revealed a negative relation (p = 0.025; estimate −0.07) between FD_Q and CPA_r. No significant effect of FD_Q on vertical canopy structure has been found (p > 0.05). The findings are discussed with regard to resource partitioning and niche differentiation of canopy and sub-canopy species. Successional stage positively impacted the CV of crown length (p = 0.019; estimate 0.03) but did not affect other response variables. The sloping terrain strongly influenced vertical canopy structure as revealed by the significant effect of slope inclination on CV of crown length (p = 0.004; estimate −0.05) and of slope aspect on CV of mean crown height (p = 0.036; estimate −0.03). The high complexity of vertical crown positioning depending on the heterogeneous sloping terrain of the study area may have obscured relations of FD_Q to vertical canopy structure.     

Gross primary production responses to warming,elevated CO2, and irrigation: quantifying the drivers of ecosystem physiology in a semiarid grassland

Determining whether the terrestrial biosphere will be a source or sink of carbon (C) under a future climate of elevated CO₂ (eCO₂) and warming requires accurate quantification of gross primary production (GPP), the largest flux of C in the global C cycle. We evaluated 6 years (2007–2012) of flux‐derived GPP data from the Prairie Heating and CO₂ Enrichment (PHACE) experiment, situated in a grassland in Wyoming, USA. The GPP data were used to calibrate a light response model whose basic formulation has been successfully used in a variety of ecosystems. The model was extended by modeling maximum photosynthetic rate (A_max) and light‐use efficiency (Q) as functions of soil water, air temperature, vapor pressure deficit, vegetation greenness, and nitrogen at current and antecedent (past) timescales. The model fits the observed GPP well (R² = 0.79), which was confirmed by other model performance checks that compared different variants of the model (e.g. with and without antecedent effects). Stimulation of cumulative 6‐year GPP by warming (29%, P = 0.02) and eCO₂ (26%, P = 0.07) was primarily driven by enhanced C uptake during spring (129%, P = 0.001) and fall (124%, P = 0.001), respectively, which was consistent across years. Antecedent air temperature (Tair_ant) and vapor pressure deficit (VPD_ant) effects on A_max (over the past 3–4 days and 1–3 days, respectively) were the most significant predictors of temporal variability in GPP among most treatments. The importance of VPD_ant suggests that atmospheric drought is important for predicting GPP under current and future climate; we highlight the need for experimental studies to identify the mechanisms underlying such antecedent effects. Finally, posterior estimates of cumulative GPP under control and eCO₂ treatments were tested as a benchmark against 12 terrestrial biosphere models (TBMs). The narrow uncertainties of these data‐driven GPP estimates suggest that they could be useful semi‐independent data streams for validating TBMs.     

Adaptation of soil microbial growth to temperature: Using a tropical elevation gradient to predict future changes

Terrestrial biogeochemical feedbacks to the climate are strongly modulated by the temperature response of soil microorganisms. Tropical forests, in particular, exert a major influence on global climate because they are the most productive terrestrial ecosystem. We used an elevation gradient across tropical forest in the Andes (a gradient of 20°C mean annual temperature, MAT), to test whether soil bacterial and fungal community growth responses are adapted to long‐term temperature differences. We evaluated the temperature dependency of soil bacterial and fungal growth using the leucine‐ and acetate‐incorporation methods, respectively, and determined indices for the temperature response of growth: Q₁₀ (temperature sensitivity over a given 10oC range) and T_min (the minimum temperature for growth). For both bacterial and fungal communities, increased MAT (decreased elevation) resulted in increases in Q₁₀ and T_min of growth. Across a MAT range from 6°C to 26°C, the Q₁₀ and T_min varied for bacterial growth (Q_10–20 = 2.4 to 3.5; T_min = ?8°C to ?1.5°C) and fungal growth (Q_10–20 = 2.6 to 3.6; T_min = ?6°C to ?1°C). Thus, bacteria and fungi did not differ significantly in their growth temperature responses with changes in MAT. Our findings indicate that across natural temperature gradients, each increase in MAT by 1°C results in increases in T_min of microbial growth by approximately 0.3°C and Q_10–20 by 0.05, consistent with long‐term temperature adaptation of soil microbial communities. A 2°C warming would increase microbial activity across a MAT gradient of 6°C to 26°C by 28% to 15%, respectively, and temperature adaptation of microbial communities would further increase activity by 1.2% to 0.3%. The impact of warming on microbial activity, and the related impact on soil carbon cycling, is thus greater in regions with lower MAT. These results can be used to predict future changes in the temperature response of microbial activity over different levels of warming and over large temperature ranges, extending to tropical regions.     

Cytotoxic peptides with insulin‐releasing activities from skin secretions of the Italian stream frog Rana italica (Ranidae)

Peptidomic analysis of norepinephrine‐stimulated skin secretions from Italian stream frog Rana italica led to the purification and characterization of two host‐defense peptides differing by a single amino acid residue belonging to the brevinin‐1 family (brevinin‐1ITa and ‐1ITb), a peptide belonging to the temporin family (temporin‐ITa) and a component identified as prokineticin Bv8. The secretions contained relatively high concentrations of the methionine‐sulphoxide forms of brevinin‐1ITa and ‐1ITb suggesting that these peptides may have a role as antioxidants in the skin of this montane frog. Brevinin‐1ITa (IVPFLLGMVPKLVCLITKKC) displayed potent cytotoxicity against non‐small cell lung adenocarcinoma A549 cells (LC₅₀ = 18 μM), breast adenocarcinoma MDA‐MB‐231 cells (LC₅₀ = 8 μM) and colorectal adenocarcinoma HT‐29 cells (LC₅₀ = 18 μM), but the peptide was also strongly hemolytic against mouse erythrocytes (LC₅₀ = 7 μM). Temporin‐ITa (VFLGAIAQALTSLLGKL.NH₂) was between three and fivefold less potent against these cells. Brevinin‐1ITa inhibited growth of both Gram‐positive Staphylococcus epidermidis and Gram‐negative Escherichia coli as well as a strain of the opportunist yeast pathogen Candida parapsilosis, whereas temporin‐ITa was active only against S. epidermidis and C. parapsilosis. Both peptides stimulated the release of insulin from BRIN‐BD11 clonal β‐cells at concentrations ≥1 nM, but brevinin‐1ITa was cytotoxic to the cells at concentrations ≥3 μM. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.     

Molecular design and optimization of hepatic cancer SLP76‐derived PLCγ1 SH3‐binding peptide with the systematic N‐substitution of peptide PXXP motif

The phospholipase Cγ1 (PLCγ1) is essential for T‐cell signaling and activation in hepatic cancer immune response, which has a regulatory Src homology 3 (SH3) domain that can specifically recognize and interact with the PXXP‐containing decapeptide segment (¹⁸⁵QP P VP P QRPM¹⁹⁴, termed as SLP76^185–194 peptide) of adaptor protein SLP76 following T‐cell receptor ligation. The isolated peptide can only bind to the PLCγ1 SH3 domain with a moderate affinity due to lack of protein context support. Instead of the traditional natural residue mutagenesis that is limited by low structural diversity and shifted target specificity, we herein attempt to improve the peptide affinity by replacing the two key proline residues Pro187 and Pro190 of SLP76^185–194 PXXP motif with nonnatural N‐substituted amino acids, as the proline is the only endogenous N‐substituted amino acid. The replacement would increase peptide flexibility but can restore peptide activity by establishing additional interactions with the domain. Structural analysis reveals that the domain pocket can be divided into a large amphipathic region and a small negatively charged region; they accommodate hydrophobic, aromatic, polar, and moderate‐sized N‐substituted amino acid types. A systematic replacement combination profile between the peptide residues Pro187 and Pro190 is created by structural modeling, dynamics simulation, and energetics analysis, from which six improved and two reduced N‐substituted peptides as well as native SLP76^185–194 peptide are identified and tested for their binding affinity to the recombinant protein of the human PLCγ1 SH3 domain using fluorescence‐based assays. Two N‐substituted peptides, SLP76^185–194(N‐Leu187/N‐Gln190) and SLP76^185–194(N‐Thr187/N‐Gln190), are designed to have high potency (K_d = 0.67 ± 0.18 and 1.7 ± 0.3 μM, respectively), with affinity improvement by, respectively, 8.5‐fold and 3.4‐fold relative to native peptide (K_d = 5.7 ± 1.2 μM).     

Hypercapnia is a Possible Determinant of the Function of the Blood-Cerebrospinal Fluid Barrier in Amyotrophic Lateral Sclerosis

Elevated cerebrospinal fluid (CSF)/serum quotients of albumin (Q_Alb) may occur in motor neuron diseases (MND) including amyotrophic lateral sclerosis (ALS), but the pathophysiologic mechanisms underlying these alterations are unclear. Evidence from animal experiments suggests that the arterial carbon dioxide level might affect the Q_Alb, i.e. the function of the blood-CSF barrier (BCB). We therefore compared basic CSF parameters in different forms of MND (ALS, n = 105; lower motor neuron diseases, n = 12; and upper motor neuron diseases, n = 7) and investigated the relationship between elevated Q_Alb and the arterial partial pressure of carbon dioxide (pCO₂) in ALS where respiratory insufficiency leads to hypercapnia in the course of the disease. Pathologic elevations of Q_Alb occurred in 32 of 124 MND patients. In ALS, Q_Alb significantly correlated with the arterial pCO₂ (r = 0.454; P = 0.001; n = 45). These data indicate that BCB dysfunction is a frequent finding in different forms of MND and may reflect distinct pathophysiological mechanisms. In ALS, an important underlying mechanism might be the influence of the arterial pCO₂ which may alter the CSF flow.     

Molecular design of Stat3‐derived peptide selectivity between BET proteins Brd2 and Brd4 in ovarian cancer

Stat3 signaling has been recognized as a potential therapeutic target of human ovarian cancer. The signaling is transducted through the peptide‐medicated interaction of Stat3 with BET family members Brd2 and Brd4 –– 2 highly homologous proteins involved in differential downstream pathways. Here, we reported a successful design of peptide selectivity between the Brd2 and Brd4. The design resulted in 3 linear peptides SMSLQCX YLGVA, QSKVLTX SYWGA, and RQCNLGX LYMNY with high or moderate selectivity for Brd2 over Brd4 (S = 3.3‐fold, 6.8‐fold, and 4.2‐fold, respectively) as compared with the native Stat3 peptide ²⁸¹HNLLRIX QFLQS²⁹² (S = 2.5‐fold). Structural analysis revealed that peptide N‐terminus and hydrogen bonds play important roles in the peptide interaction stability and specificity with Brd2 and Brd4. This study would help to establish an integrated in silico‐in vitro method for rational molecular design of peptide ligand selectivity between homologous protein receptors.     

Effect of menstrual cycle phase and hormonal treatments on evaluation of tubal patency in baboons

Background

We evaluated whether menstrual cycle phase influences the assessment of tubal patency by hysterosalpingography (HSG) in baboons.

Methods

Retrospective analysis of baseline tubal patency studies and serum estradiol (E₂) and progesterone (P4) values obtained from female baboons used as models for development of non‐surgical permanent contraception in women. The main outcome measure was bilateral tubal patency (BTP) in relationship with estradiol level.

Results

Female baboons (n = 110) underwent a single (n = 81), two (n = 26), or three (n = 3) HSG examinations. In 33/142 (23%) HSG examinations, one or both tubes showed functional occlusion (FO). The median E₂ in studies with BTP (49 pg/mL) was significantly higher than in those studies with FO (32 pg/mL, P = .005). Among 18 animals with repeat examinations where serum E₂ changed from <60 to ≥ 60 pg/mL, 13 results changed from FO to BTP (P = .0001). No sets showed a change from BTP to FO with an increase in estradiol.

Conclusion

In baboons, functional occlusion of the fallopian tube is associated with low estradiol levels, supporting a role for estrogen‐mediated relaxation of the utero‐tubal junction.     

Using field photography to study avian moult

Methods to obtain moult data from wild birds have not changed much over the last century and most studies still depend on checking museum specimens or capturing birds. Here we assess the applicability of systematic field photography for detecting and scoring moult in adult Black Skimmers Rynchops niger from southern Brazil. Moult data extracted from photographs have a high within‐ (R_GLMM = 0.98) and between‐observer repeatability (R_GLMM = 0.97) and show very good fit to current Underhill–Zucchini moult models (R² = 0.75). Photography offers the advantages of being less invasive, requiring less equipment and human effort, being feasible in areas where captures may not be possible, and causing less disturbance, so enhancing the number of sampled individuals.     

Anti‐HSV‐1 and HSV‐2 Flavonoids and a New Kaempferol Triglycoside from the Medicinal Plant Kalanchoe daigremontiana

Kalanchoe daigremontiana (Crassulaceae) is a medicinal plant native to Madagascar. The aim of this study was to investigate the flavonoid content of an aqueous leaf extract from K. daigremontiana (Kd), and assess its antiherpetic potential. The major flavonoid, kaempferol 3‐O‐β‐d ‐xylopyranosyl‐(1 → 2)‐α‐l ‐rhamnopyranoside ( 1 ), was isolated from the AcOEt fraction (Kd‐AC). The BuOH‐soluble fraction afforded quercetin 3‐O‐β‐d ‐xylopyranosyl‐(1 → 2)‐α‐l ‐rhamnopyranoside ( 2 ) and the new kaempferol 3‐O‐β‐d ‐xylopyranosyl‐(1 → 2)‐α‐l ‐rhamnopyranoside‐7‐O‐β‐d ‐glucopyranoside ( 3 ), named daigremontrioside. The crude extract, Kd‐AC fraction, flavonoids 1 and 2 were evaluated using acyclovir‐sensitive strains of HSV‐1 and HSV‐2. Kd‐AC was highly active against HSV‐1 (EC₅₀ = 0.97 μg/ml, SI > 206.1) and HSV‐2 (EC₅₀ = 0.72 μg/ml, SI > 277.7). Flavonoids 1 and 2 showed anti‐HSV‐1 (EC₅₀ = 7.4 μg/ml; SI > 27 and EC₅₀ = 5.8 μg/ml; SI > 8.6, respectively) and anti‐HSV‐2 (EC₅₀ = 9.0 μg/ml; SI > 22.2 and EC₅₀ = 36.2 μg/ml; SI > 5.5, respectively) activities, suggesting the contribution of additional substances to the antiviral activity.     

Repeatability of Quantitative Sodium Magnetic Resonance Imaging for Estimating Pseudo-Intracellular Sodium Concentration and Pseudo-Extracellular Volume Fraction in Brain at 3 T

The purpose of this study is to assess the repeatability of the quantification of pseudo-intracellular sodium concentration (C₁) and pseudo-extracellular volume fraction (α) estimated in brain in vivo using sodium magnetic resonance (MRI) at 3 T. Eleven healthy subjects were scanned twice, with two sodium MRI acquisitions (with and without fluid suppression by inversion recovery), and two double inversion recovery (DIR) proton MRI. DIR MRIs were used to create masks of gray and white matter (GM, WM), that were subsequently applied to the C₁ and α maps calculated from sodium MRI and a tissue three-compartment model, in order to measure the distributions of these two parameters in GM, WM or full brain (GM+WM) separately. The mean, median, mode, standard deviation (std), skewness and kurtosis of the C₁ and α distributions in whole GM, WM and full brain were calculated for each subject, averaged over all data, and used as parameters for the repeatability assessment. The coefficient of variation (CV) was calculated as a measure of reliability for the detection of intra-subject changes in C₁ and αfor each parameter, while intraclass correlation (ICC) was used as a measure of repeatability. It was found that the CV of most of the parameters was around 10–20% (except for C₁ kurtosis which is about 40%) for C₁ and α measurements, and that ICC was moderate to very good (0.4 to 0.9) for C₁ parameters and for some of the α parameters (mainly skewness and kurtosis). In conclusion, the proposed method could allow to reliably detect changes of 50% and above of the different measurement parameters of C₁ and αin neuropathologies (multiple sclerosis, tumor, stroke, Alzheimer’s disease) compared to healthy subjects, and that skewness and kurtosis of the distributions of C₁ and αseem to be the more sensitive parameters to these changes.     

Changes in smolt traits of Atlantic salmon (Salmo salar Linnaeus, 1758) and linkages to parr density and water temperature

Smolt traits (length, age) and timing of smolt migration of wild Atlantic salmon, Salmo salar L., were investigated in the Simojoki River, northern Baltic Sea. The aim was to determine whether they responded to changes in parr length, parr density and temperature from 2000 to 2014. Annual electrofishing surveys and smolt numbers determined parr densities by springtime trapping in the river mouth. During the smolt trapping period captured parr and smolts were aged from scales. Water temperature was measured daily. Mean length decreased from 137 mm (TL) to 129 mm among 2‐year‐old smolts, and from 150 mm to 139 mm among 3‐year‐olds. Median date of the smolt migration was 10 days earlier, from early June to late May during the study period, linked to the rise in air temperature in May at the nearby Kemi‐Tornio airport. However, the median day temperature and the mean daily water temperatures during the second (Q₂) and third (Q₃) migration quartiles did not change. This implied that migration began when a suitable water temperature was reached, independent of the date.     

Limited mass‐independent individual variation in resting metabolic rate in a wild population of snow voles (Chionomys nivalis)

Resting metabolic rate (RMR) is a potentially important axis of physiological adaptation to the thermal environment. However, our understanding of the causes and consequences of individual variation in RMR in the wild is hampered by a lack of data, as well as analytical challenges. RMR measurements in the wild are generally characterized by large measurement errors and a strong dependency on mass. The latter is problematic when assessing the ability of RMR to evolve independently of mass. Mixed models provide a powerful and flexible tool to tackle these challenges, but they have rarely been used to estimate repeatability of mass‐independent RMR from field data. We used respirometry to obtain repeated measurements of RMR in a long‐term study population of snow voles (Chionomys nivalis) inhabiting an environment subject to large circadian and seasonal fluctuations in temperature. Using both uni‐ and bivariate mixed models, we quantify individual repeatability in RMR and decompose repeatability into mass‐dependent and mass‐independent components, while accounting for measurement error. RMR varies among individuals, that is, is repeatable (R = .46) and strongly co‐varies with BM. Indeed, much of the repeatability of RMR is attributable to individual variation in BM, and the repeatability of mass‐independent RMR is reduced by 41% to R = .27. These empirical results suggest that the evolutionary potential of RMR independent of mass may be severely constrained. This study illustrates how to leverage bivariate mixed models to model field data for metabolic traits, correct for measurement error and decompose the relative importance of mass‐dependent and mass‐independent physiological variation.     

Spatial and temporal variation in nitrogen fixation and its importance to phytoplankton in phosphorus‐rich lakes

相似文献   

High-level SUMO-mediated fusion expression of ABP-dHC-cecropin A from multiple joined genes in Escherichia coli

The antimicrobial peptide ABP-dHC-cecropin A is a small cationic peptide with potent activity against a wide range of bacterial species. Evidence of antifungal activity has also been suggested; however, evaluation of this peptide has been limited due to the low expression of cecropin proteins in Escherichia coli. To improve the expression level of ABP-dHC-cecropin A in E. coli, tandem repeats of the ABP-dHC-cecropin A gene were constructed and expressed as fusion proteins (SUMO-nABP-dHC-cecropin, n = 1, 2, 3, 4) via pSUMO-nABP-dHC-cecropin A vectors (n = 1, 2, 3, 4). Comparison of the expression levels of soluble SUMO-nABP-dHC-cecropin A fusion proteins (n = 1, 2, 3, 4) suggested that BL21 (DE3)/pSUMO-3ABP-dHC-cecropin A is an ideal recombinant strain for ABP-dHC-cecropin A production. Under the selected conditions of cultivation and isopropylthiogalactoside (IPTG) induction, the expression level of ABP-dHC-cecropin A was as high as 65 mg/L, with ∼21.3% of the fusion protein in soluble form. By large-scale fermentation, protein production reached nearly 300 mg/L, which is the highest yield of ABP-dHC-cecropin A reported to date. In antibacterial experiments, the efficacy was approximately the same as that of synthetic ABP-dHC-cecropin A. This method provides a novel and effective means of producing large amounts of ABP-dHC-cecropin A.     

Synthesis and Anticancer Activity of 3‐(Substituted Aroyl)‐4‐(3,4,5‐trimethoxyphenyl)‐1H‐pyrrole Derivatives

A series of 3‐(substituted aroyl)‐4‐(3,4,5‐trimethoxyphenyl)‐1H‐pyrrole derivatives were synthesized and determined for their anticancer activity against eleven cancer cell lines and two normal tissue cell lines using MTT assay. Among the synthesized compounds, compound 3f was the most potent compound against A375, CT‐26, HeLa, MGC80‐3, NCI‐H460 and SGC‐7901 cells (IC₅₀ = 8.2 – 31.7 μm ); 3g , 3n and 3a were the most potent compounds against CHO (IC₅₀ = 8.2 μm ), HCT‐15 (IC₅₀ = 21 μm ) and MCF‐7 cells (IC₅₀ = 18.7 μm ), respectively. Importantly, all the target compounds showed no cytotoxicity towards the normal tissue cell (IC₅₀ > 100 μm ). Thus, these compounds with the potent anticancer activity and low toxicity have potential for the development of new anticancer chemotherapy agents.     

Safety and efficacy of cryoballoon versus radiofrequency ablation for atrial fibrillation in elderly patients: A real-world evidence

BackgroundThere are limited data describing the experience of radiofrequency (RF) vs. cryoballoon (CB) ablation for atrial fibrillation (AF) among elderly patients in the United States.MethodsWe conducted a retrospective analysis of patients ≥75 years of age undergoing index RF vs. CB ablation between January 2014 and May 2020 at our center. The choice of ablation technique was left to the operator's discretion. Major complications and efficacy, defined as freedom from any atrial tachyarrhythmia (ATA) lasting ≥30 s after one year of follow-up, were assessed in patients with index RF vs. CB ablation.ResultsIn our cohort of 186 patients, the median age was 78 (76–81) years, 54.8% were men, and 39.2% had persistent AF. The median CHA₂DS₂-VASc score was 4 (3–4), while the median duration of AF was 3 (1–7) years. The majority (n = 112, 60.2%) underwent RF ablation. The median procedure time was significantly lower in CB group (197 vs 226.5 min, p=<0.01). The incidence of complications was similar in the two sub-groups (RF: 1.8% vs. CB: 2.7%, p = 0.67). Similarly, arrhythmia-free survival rate on antiarrhythmic drugs at 1-year follow-up remained statistically comparable (63.4% vs. 68.9%, p = 0.33) between patients receiving RF vs. CB ablation.ConclusionThe safety and efficacy of RF vs. CB ablation for AF remained comparable in our cohort of patients older than 75 years. CB ablation was associated with a shorter procedure time.     

Speciation genomics and a role for the Z chromosome in the early stages of divergence between Mexican ducks and mallards

Speciation is a continuous and dynamic process, and studying organisms during the early stages of this process can aid in identifying speciation mechanisms. The mallard (Anas platyrhynchos) and Mexican duck (A. [p.] diazi) are two recently diverged taxa with a history of hybridization and controversial taxonomy. To understand their evolutionary history, we conducted genomic scans to characterize patterns of genetic diversity and divergence across the mitochondrial DNA (mtDNA) control region, 3523 autosomal loci and 172 Z‐linked sex chromosome loci. Between the two taxa, Z‐linked loci (Φ_ST = 0.088) were 5.2 times more differentiated than autosomal DNA (Φ_ST = 0.017) but comparable to mtDNA (Φ_ST = 0.092). This elevated Z differentiation deviated from neutral expectations inferred from simulated data that incorporated demographic history and differences in effective population sizes between marker types. Furthermore, 3% of Z‐linked loci, compared to <0.1% of autosomal loci, were detected as outlier loci under divergent selection with elevated relative (Φ_ST) and absolute (d_XY) estimates of divergence. In contrast, the ratio of Z‐linked and autosomal differentiation among the seven Mexican duck sampling locations was close to 1:1 (Φ_ST = 0.018 for both markers). We conclude that between mallards and Mexican ducks, divergence at autosomal markers is largely neutral, whereas greater divergence on the Z chromosome (or some portions thereof) is likely the product of selection that has been important in speciation. Our results contribute to a growing body of literature indicating elevated divergence on the Z chromosome and its likely importance in avian speciation.