Phylogenetic insight into ABCE gene subfamily in plants.

ATP-BINDING CASSETTE SUBFAMILY E MEMBER (ABCE) proteins are one of the most conserved proteins across eukaryotes and archaea. Yeast and most animals possess a single ABCE gene encoding the critical translational factor ABCE1. In several plant species, including Arabidopsis thaliana and Oryza sativa, two or more ABCE gene copies have been identified, however information related to plant ABCE gene family is still missing. In this study we retrieved ABCE gene sequences of 76 plant species from public genome databases and comprehensively analyzed them with the reference to A. thaliana ABCE2 gene (AtABCE2). Using bioinformatic approach we assessed the conservation and phylogeny of plant ABCEs. In addition, we performed haplotype analysis of AtABCE2 and its paralogue AtABCE1 using genomic sequences of 1,135 A. thaliana ecotypes. Plant ABCE proteins showed overall high sequence conservation, sharing at least 78% of amino acid sequence identity with AtABCE2. We found that over half of the selected species have two to eight ABCE genes, suggesting that in plants ABCE genes can be classified as a low-copy gene family, rather than a single-copy gene family. The phylogenetic trees of ABCE protein sequences and the corresponding coding sequences demonstrated that Brassicaceae and Poaceae families have independently undergone lineage-specific split of the ancestral ABCE gene. Other plant species have gained ABCE gene copies through more recent duplication events. We also noticed that ploidy level but not ancient whole genome duplications experienced by a species impacts ABCE gene family size. Deeper analysis of AtABCE2 and AtABCE1 from 1,135 A. thaliana ecotypes revealed four and 35 non-synonymous SNPs, respectively. The lower natural variation in AtABCE2 compared to AtABCE1 is in consistence with its crucial role for plant viability. Overall, while the sequence of the ABCE protein family is highly conserved in the plant kingdom, many plants have evolved to have more than one copy of this essential translational factor.

User-friendly analysis of droplet array images.

Water-in-oil droplets allow performing massive experimental parallelization and high-throughput studies, such as single-cell experiments. However, analyzing such vast arrays of droplets usually requires advanced expertise and sophisticated workflow tools, which limits accessibility for a wider user base in the fields of chemistry and biology. Thus, there is a need for more user-friendly tools for droplet analysis. In this article, we deliver a set of analytical pipelines for user-friendly analysis of typical scenarios in droplet experiments. We built pipelines that combine various open-source image-analysis software with a custom-developed data processing tool called "EasyFlow". Our pipelines are applicable to the typical experimental scenarios that users encounter when working with droplets: i) mono- and polydisperse droplets, ii) brightfield and fluorescent images, iii) droplet and object detection, iv) signal profile of droplets and objects (e.g., fluorescence).

Sex-dependent expression levels of VAV1 and P2X7 in PBMC of multiple sclerosis patients.

Multiple sclerosis (MS) is an inflammatory autoimmune disorder of the central nervous system and the leading cause of progressive neurological disability in young adults. It decreases the patient's lifespan by about 10 years and affects women more than men. No medication entirely restricts or reverses neurological degradation. However, early diagnosis and treatment increase the possibility of a better outcome. To identify new MS biomarkers, we tested the expression of six potential markers (P2X4, P2X7, CXCR4, RGS1, RGS16 and VAV1) using qPCR in peripheral blood mononuclear cells (PBMC) of MS patients treated with interferon ß (IFNß), with glatiramer acetate (GA) or untreated. We showed that P2X7 and VAV1 are significantly induced in MS patients. In contrast, the expression of P2X4, CXCR4, RGS1 and RGS16 was not significantly modified by MS in PBMC. P2X7 and VAV1 are essentially induced in female patients, suggesting these markers are connected to sex-specific mechanisms. Strikingly, VAV1 expression is higher in healthy women than healthy men and IFNß treatment of MS reduced VAV1 expression in female MS patients while it up-regulated VAV1 in male MS patients. Our data point to the differential, sex-dependent value of MS markers and treatment effects. Although rgs16 expression in PBMC was not a valid MS marker in patients, the strong upregulation of P2X4 and P2X7 induced in the spinal cord of WT mice by EAE was abrogated in rgs16KO mice suggesting that rgs16 is required for P2X4 and P2X7 induction by neurological diseases.

Single-cell RNA-seq analysis and cell-cluster deconvolution of the human preovulatory follicular fluid cells provide insights into the pathophysiology of ovarian hyporesponse.

Reduction in responsiveness to gonadotropins or hyporesponsiveness may lead to the failure of in vitro fertilization (IVF), due to a low number of retrieved oocytes. The ovarian sensitivity index (OSI) is used to reflect the ovarian responsiveness to gonadotropin stimulation before IVF. Although introduced to clinical practice already years ago, its usefulness to predict clinical outcomes requires further research. Nevertheless, pathophysiological mechanisms of ovarian hyporesponse, along with advanced maternal age and in younger women, have not been fully elucidated. Follicles consist of multiple cell types responsible for a repertoire of biological processes including responding to pituitary gonadotropins necessary for follicle growth and oocyte maturation as well as ovulation. Encouraging evidence suggests that hyporesponse could be influenced by many contributing factors, therefore, investigating the variability of ovarian follicular cell types and their gene expression in hyporesponders is highly informative for increasing their prognosis for IVF live birth. Due to advancements in single-cell analysis technologies, the role of somatic cell populations in the development of infertility of ovarian etiology can be clarified. Here, somatic cells were collected from the fluid of preovulatory ovarian follicles of patients undergoing IVF, and RNA-seq was performed to study the associations between OSI and gene expression. We identified 12 molecular pathways differentially regulated between hypo- and normoresponder patient groups (FDR<0.05) from which extracellular matrix organization, post-translational protein phosphorylation, and regulation of Insulin-like Growth Factor (IGF) transport and uptake by IGF Binding Proteins were regulated age-independently. We then generated single-cell RNA-seq data from matching follicles revealing 14 distinct cell clusters. Using cell cluster-specific deconvolution from the bulk RNA-seq data of 18 IVF patients we integrated the datasets as a novel approach and discovered that the abundance of three cell clusters significantly varied between hypo- and normoresponder groups suggesting their role in contributing to the deviations from normal ovarian response to gonadotropin stimulation. Our work uncovers new information regarding the differences in the follicular gene expression between hypo- and normoresponders. In addition, the current study fills the gap in understanding the inter-patient variability of cell types in human preovulatory follicles, as revealed by single-cell analysis of follicular fluid cells.

The transcription factor network of E. coli steers global responses to shifts in RNAP concentration.

The robustness and sensitivity of gene networks to environmental changes is critical for cell survival. How gene networks produce specific, chronologically ordered responses to genome-wide perturbations, while robustly maintaining homeostasis, remains an open question. We analysed if short- and mid-term genome-wide responses to shifts in RNA polymerase (RNAP) concentration are influenced by the known topology and logic of the transcription factor network (TFN) of Escherichia coli. We found that, at the gene cohort level, the magnitude of the single-gene, mid-term transcriptional responses to changes in RNAP concentration can be explained by the absolute difference between the gene's numbers of activating and repressing input transcription factors (TFs). Interestingly, this difference is strongly positively correlated with the number of input TFs of the gene. Meanwhile, short-term responses showed only weak influence from the TFN. Our results suggest that the global topological traits of the TFN of E. coli shape which gene cohorts respond to genome-wide stresses.

Improved chromosome-level genome assembly of the Glanville fritillary butterfly (Melitaea cinxia) integrating Pacific Biosciences long reads and a high-density linkage map.

BACKGROUND: The Glanville fritillary (Melitaea cinxia) butterfly is a model system for metapopulation dynamics research in fragmented landscapes. Here, we provide a chromosome-level assembly of the butterfly's genome produced from Pacific Biosciences sequencing of a pool of males, combined with a linkage map from population crosses. RESULTS: The final assembly size of 484 Mb is an increase of 94 Mb on the previously published genome. Estimation of the completeness of the genome with BUSCO indicates that the genome contains 92-94% of the BUSCO genes in complete and single copies. We predicted 14,810 genes using the MAKER pipeline and manually curated 1,232 of these gene models. CONCLUSIONS: The genome and its annotated gene models are a valuable resource for future comparative genomics, molecular biology, transcriptome, and genetics studies on this species.

Hsa-mir-548 family expression in human reproductive tissues.

BACKGROUND: Hsa-miR-548ba expressed in ovarian granulosa cells targets PTEN and LIFR, which are essential for ovarian follicle activation and growth. The expression pattern of hsa-miR-548ba correlates with its host gene follicle-stimulating hormone receptor (FSHR), and FSH has a positive influence on hsa-miR-548ba expression. However, hsa-miR-548ba is a member of a large hsa-mir-548 family with potentially overlapping targets. The current study aims to investigate the co-expression of hsa-mir-548 family members in FSHR-positive reproductive tissues and to explore the potential co-regulation of pathways. RESULTS: For the above-described analysis, small RNA sequencing data from public data repositories were used. Sequencing results revealed that hsa-miR-548ba was expressed at the highest level in the ovarian granulosa cells and uterine myometrial samples together with another twelve and one hsa-miR-548 family members, respectively. Pathway enrichment analysis of microRNA targets in the ovarian samples revealed the hsa-miR-548ba and hsa-miR-548b-5p co-regulation of RAB geranylgeranylation in mural granulosa cells. Moreover, other hsa-mir-548 family members co-regulate pathways essential for ovarian functions (PIP3 activates AKT signalling and signalling by ERBB4). In addition to hsa-miR-548ba, hsa-miR-548o-3p is expressed in the myometrium, which separately targets the peroxisome proliferator-activated receptor alpha (PPARA) pathway. CONCLUSION: This study reveals that hsa-mir-548 family members are expressed in variable combinations in the reproductive tract, where they potentially fulfil different regulatory roles. The results provide a reference for further studies of the hsa-mir-548 family role in the reproductive tract.

Investigation of Different Free Image Analysis Software for High-Throughput Droplet Detection.

Droplet microfluidics has revealed innovative strategies in biology and chemistry. This advancement has delivered novel quantification methods, such as droplet digital polymerase chain reaction (ddPCR) and an antibiotic heteroresistance analysis tool. For droplet analysis, researchers often use image-based detection techniques. Unfortunately, the analysis of images may require specific tools or programming skills to produce the expected results. In order to address the issue, we explore the potential use of standalone freely available software to perform image-based droplet detection. We select the four most popular software and classify them into rule-based and machine learning-based types after assessing the software's modules. We test and evaluate the software's (i) ability to detect droplets, (ii) accuracy and precision, and (iii) overall components and supporting material. In our experimental setting, we find that the rule-based type of software is better suited for image-based droplet detection. The rule-based type of software also has a simpler workflow or pipeline, especially aimed for non-experienced users. In our case, CellProfiler (CP) offers the most user-friendly experience for both single image and batch processing analyses.

Cellular, Extracellular and Extracellular Vesicular miRNA Profiles of Pre-Ovulatory Follicles Indicate Signaling Disturbances in Polycystic Ovaries.

Cell-free RNAs have the potential to act as a means of gene expression regulation between cells and are therefore used as diagnostic markers describing the state of tissue environment. The origin and functions of such RNAs in human ovarian follicle, the environment of oocyte maturation, are unclear. The current study investigates the difference in the microRNA profiles of fertile women and polycystic ovary syndrome (PCOS) patients in three compartments from the same preovulatory follicle: mural granulosa cells (MGC), cell-free follicular fluid (FF), and extracellular vesicles (EV) of the FF by small RNA sequencing. In silico analysis was used for the prediction and over-representation of targeted pathways for the detected microRNAs. PCOS follicles were distinguished from normal tissue by the differential expression of 30 microRNAs in MGC and 10 microRNAs in FF (FDR < 0.1) that commonly regulate cytokine signaling pathways. The concentration of EV-s was higher in the FF of PCOS patients (p = 0.04) containing eight differentially expressed microRNAs (p < 0.05). In addition, we present the microRNA profiles of MGC, FF, and EV in the fertile follicle and demonstrate that microRNAs loaded into EVs target mRNAs of distinct signaling pathways in comparison to microRNAs in FF. To conclude, the three follicular compartments play distinct roles in the signaling disturbances associated with PCOS.

Molecular profile of the rat peri-infarct region four days after stroke: Study with MANF.

The peri-infarct region after ischemic stroke is the anatomical location for many of the endogenous recovery processes; however, -the molecular events in the peri-infarct region remain poorly characterized. In this study, we examine the molecular profile of the peri-infarct region on post-stroke day four, a time when reparative processes are ongoing. We used a multiomics approach, involving RNA sequencing, and mass spectrometry-based proteomics and metabolomics to characterize molecular changes in the peri-infarct region. We also took advantage of our previously developed method to express transgenes in the peri-infarct region where self-complementary adeno-associated virus (AAV) vectors were injected into the brain parenchyma on post-stroke day 2. We have previously used this method to show that mesencephalic astrocyte-derived neurotrophic factor (MANF) enhances functional recovery from stroke and recruits phagocytic cells to the peri-infarct region. Here, we first analyzed the effects of stroke to the peri-infarct region on post-stroke day 4 in comparison to sham-operated animals, finding that strokeinduced changes in 3345 transcripts, 341 proteins, and 88 metabolites. We found that after stroke, genes related to inflammation, proliferation, apoptosis, and regeneration were upregulated, whereas genes encoding neuroactive ligand receptors and calcium-binding proteins were downregulated. In proteomics, we detected upregulation of proteins related to protein synthesis and downregulation of neuronal proteins. Metabolomic studies indicated that in after stroke tissue there is an increase in saccharides, sugar phosphates, ceramides and free fatty acids and a decrease of adenine, hypoxantine, adenosine and guanosine. We then compared the effects of post-stroke delivery of AAV1-MANF to AAV1-eGFP (enhanced green fluorescent protein). MANF administration increased the expression of 77 genes, most of which were related to immune response. In proteomics, MANF administration reduced S100A8 and S100A9 protein levels. In metabolomics, no significant differences between MANF and eGFP treatment were detected, but relative to sham surgery group, most of the changes in lipids were significant in the AAV-eGFP group only. This work describes the molecular profile of the peri-infarct region during recovery from ischemic stroke, and establishes a resource for further stroke studies. These results provide further support for parenchymal MANF as a modulator of phagocytic function.

Evolutionary Origin of the P2X7 C-ter Region: Capture of an Ancient Ballast Domain by a P2X4-Like Gene in Ancient Jawed Vertebrates.

P2X purinergic receptors are extracellular ATP-gated ion channel receptors present on the cell plasma membrane. P2X receptors have been found in Metazoa, fungi, amoebas, and in plants. In mammals, P2X7 is expressed by a large number of cell types and is involved in inflammation and immunity. Remarkably, P2X7 does not desensitize as other P2X do, a feature linked to a "C-cysteine anchor" intra-cytoplasmic motif encoded by exon 11. Another specific feature of P2X7 is its C-terminal cytoplasmic ballast domain (exon 13) which contains a zinc (Zn) coordinating cysteine motif and a GDP-binding region. To determine the origin of P2X7, we analyzed and compared sequences and protein motifs of the C-terminal intra-cytoplasmic region across all main groups of Metazoa. We identified proteins with typical ballast domains, sharing a remarkably conserved Zn-coordinating cysteine motif. Apart from vertebrates, these ballast domains were not associated with a typical P2X architecture. These results strongly suggest that P2X7 resulted from the fusion of a P2X gene, highly similar to P2X4, with an exon encoding a ballast domain. Our work brings new evidence on the origin of the P2X7 purinergic receptor and identifies the Zn-coordinating cysteine domain as the fundamental feature of the ancient ballast fold.

Droplet-based digital antibiotic susceptibility screen reveals single-cell clonal heteroresistance in an isogenic bacterial population.

Since antibiotic resistance is a major threat to global health, recent observations that the traditional test of minimum inhibitory concentration (MIC) is not informative enough to guide effective antibiotic treatment are alarming. Bacterial heteroresistance, in which seemingly susceptible isogenic bacterial populations contain resistant sub-populations, underlies much of this challenge. To close this gap, here we developed a droplet-based digital MIC screen that constitutes a practical analytical platform for quantifying the single-cell distribution of phenotypic responses to antibiotics, as well as for measuring inoculum effect with high accuracy. We found that antibiotic efficacy is determined by the amount of antibiotic used per bacterial colony forming unit (CFU), not by the absolute antibiotic concentration, as shown by the treatment of beta-lactamase-carrying Escherichia coli with cefotaxime. We also noted that cells exhibited a pronounced clustering phenotype when exposed to near-inhibitory amounts of cefotaxime. Overall, our method facilitates research into the interplay between heteroresistance and antibiotic efficacy, as well as research into the origin and stimulation of heterogeneity by exposure to antibiotics. Due to the absolute bacteria quantification in this digital assay, our method provides a platform for developing reference MIC assays that are robust against inoculum-density variations.

ELIMÄKI Locus Is Required for Vertical Proprioceptive Response in Birch Trees.

Tree architecture has evolved to support a top-heavy above-ground biomass, but this integral feature poses a weight-induced challenge to trunk stability. Maintaining an upright stem is expected to require vertical proprioception through feedback between sensing stem weight and responding with radial growth. Despite its apparent importance, the principle by which plant stems respond to vertical loading forces remains largely unknown. Here, by manipulating the stem weight of downy birch (Betula pubescens) trees, we show that cambial development is modulated systemically along the stem. We carried out a genetic study on the underlying regulation by combining an accelerated birch flowering program with a recessive mutation at the ELIMÄKI locus (EKI), which causes a mechanically defective response to weight stimulus resulting in stem collapse after just 3 months. We observed delayed wood morphogenesis in eki compared with WT, along with a more mechanically elastic cambial zone and radial compression of xylem cell size, indicating that rapid tissue differentiation is critical for cambial growth under mechanical stress. Furthermore, the touch-induced mechanosensory pathway was transcriptionally misregulated in eki, indicating that the ELIMÄKI locus is required to integrate the weight-growth feedback regulation. By studying this birch mutant, we were able to dissect vertical proprioception from the gravitropic response associated with reaction wood formation. Our study provides evidence for both local and systemic responses to mechanical stimuli during secondary plant development.

Red-Brown Pigmentation of Acidipropionibacterium jensenii Is Tied to Haemolytic Activity and cyl-Like Gene Cluster.

The novel Acidipropionibacterium genus encompasses species of industrial importance but also those associated with food spoilage. In particular, Acidipropionibacterium acidipropionici, Acidipropionibacterium thoenii, and Acidipropionibacterium jensenii play an important role in food fermentation, as biopreservatives, or as potential probiotics. Notably, A. jensenii and A. thoenii can cause brown spot defects in Swiss-type cheeses, which have been tied to the rhamnolipid pigment granadaene. In the pathogenic bacterium Streptococcus agalactiae, production of granadaene depends on the presence of a cyl gene cluster, an important virulence factor linked with haemolytic activity. Here, we show that the production of granadaene in pigmented Acidipropionibacterium, including A. jensenii, A. thoenii, and Acidipropionibacterium virtanenii, is tied to haemolytic activity and the presence of a cyl-like gene cluster. Furthermore, we propose a PCR-based test, which allows pinpointing acidipropionibacteria with the cyl-like gene cluster. Finally, we present the first two whole genome sequence analyses of the A. jensenii strains as well as testing phenotypic characteristics important for industrial applications. In conclusion, the present study sheds light on potential risks associated with the presence of pigmented Acidipropionibacterium strains in food fermentation. In addition, the results presented here provide ground for development of a quick and simple diagnostic test instrumental in avoiding potential negative effects of Acidipropionibacterium strains with haemolytic activity on food quality.

gapFinisher: A reliable gap filling pipeline for SSPACE-LongRead scaffolder output.

Unknown sequences, or gaps, are present in many published genomes across public databases. Gap filling is an important finishing step in de novo genome assembly, especially in large genomes. The gap filling problem is nontrivial and while there are many computational tools partially solving the problem, several have shortcomings as to the reliability and correctness of the output, i.e. the gap filled draft genome. SSPACE-LongRead is a scaffolding tool that utilizes long reads from multiple third-generation sequencing platforms in finding links between contigs and combining them. The long reads potentially contain sequence information to fill the gaps created in the scaffolding, but SSPACE-LongRead currently lacks this functionality. We present an automated pipeline called gapFinisher to process SSPACE-LongRead output to fill gaps after the scaffolding. gapFinisher is based on the controlled use of a previously published gap filling tool FGAP and works on all standard Linux/UNIX command lines. We compare the performance of gapFinisher against two other published gap filling tools PBJelly and GMcloser. We conclude that gapFinisher can fill gaps in draft genomes quickly and reliably. In addition, the serial design of gapFinisher makes it scale well from prokaryote genomes to larger genomes with no increase in the computational footprint.

Notum produced by Paneth cells attenuates regeneration of aged intestinal epithelium.

A decline in stem cell function impairs tissue regeneration during ageing, but the role of the stem-cell-supporting niche in ageing is not well understood. The small intestine is maintained by actively cycling intestinal stem cells that are regulated by the Paneth cell niche1,2. Here we show that the regenerative potential of human and mouse intestinal epithelium diminishes with age owing to defects in both stem cells and their niche. The functional decline was caused by a decrease in stemness-maintaining Wnt signalling due to production of Notum, an extracellular Wnt inhibitor, in aged Paneth cells. Mechanistically, high activity of mammalian target of rapamycin complex 1 (mTORC1) in aged Paneth cells inhibits activity of peroxisome proliferator activated receptor α (PPAR-α)3, and lowered PPAR-α activity increased Notum expression. Genetic targeting of Notum or Wnt supplementation restored function of aged intestinal organoids. Moreover, pharmacological inhibition of Notum in mice enhanced the regenerative capacity of aged stem cells and promoted recovery from chemotherapy-induced damage. Our results reveal a role of the stem cell niche in ageing and demonstrate that targeting of Notum can promote regeneration of aged tissues.

Genome description of Phlebia radiata 79 with comparative genomics analysis on lignocellulose decomposition machinery of phlebioid fungi.

BACKGROUND: The white rot fungus Phlebia radiata, a type species of the genus Phlebia, is an efficient decomposer of plant cell wall polysaccharides, modifier of softwood and hardwood lignin, and is able to produce ethanol from various waste lignocellulose substrates. Thus, P. radiata is a promising organism for biotechnological applications aiming at sustainable utilization of plant biomass. Here we report the genome sequence of P. radiata isolate 79 originally isolated from decayed alder wood in South Finland. To better understand the evolution of wood decay mechanisms in this fungus and the Polyporales phlebioid clade, gene content and clustering of genes encoding specific carbohydrate-active enzymes (CAZymes) in seven closely related fungal species was investigated. In addition, other genes encoding proteins reflecting the fungal lifestyle including peptidases, transporters, small secreted proteins and genes involved in secondary metabolism were identified in the genome assembly of P. radiata. RESULTS: The PACBio sequenced nuclear genome of P. radiata was assembled to 93 contigs with 72X sequencing coverage and annotated, revealing a dense genome of 40.4 Mbp with approximately 14 082 predicted protein-coding genes. According to functional annotation, the genome harbors 209 glycoside hydrolase, 27 carbohydrate esterase, 8 polysaccharide lyase, and over 70 auxiliary redox enzyme-encoding genes. Comparisons with the genomes of other phlebioid fungi revealed shared and specific properties among the species with seemingly similar saprobic wood-decay lifestyles. Clustering of especially GH10 and AA9 enzyme-encoding genes according to genomic localization was discovered to be conserved among the phlebioid species. In P. radiata genome, a rich repertoire of genes involved in the production of secondary metabolites was recognized. In addition, 49 genes encoding predicted ABC proteins were identified in P. radiata genome together with 336 genes encoding peptidases, and 430 genes encoding small secreted proteins. CONCLUSIONS: The genome assembly of P. radiata contains wide array of carbohydrate polymer attacking CAZyme and oxidoreductase genes in a composition identifiable for phlebioid white rot lifestyle in wood decomposition, and may thus serve as reference for further studies. Comparative genomics also contributed to enlightening fungal decay mechanisms in conversion and cycling of recalcitrant organic carbon in the forest ecosystems.

Alternative oxidase-mediated respiration prevents lethal mitochondrial cardiomyopathy.

Alternative oxidase (AOX) is a non-mammalian enzyme that can bypass blockade of the complex III-IV segment of the respiratory chain (RC). We crossed a Ciona intestinalis AOX transgene into RC complex III (cIII)-deficient Bcs1lp.S78G knock-in mice, displaying multiple visceral manifestations and premature death. The homozygotes expressing AOX were viable, and their median survival was extended from 210 to 590 days due to permanent prevention of lethal cardiomyopathy. AOX also prevented renal tubular atrophy and cerebral astrogliosis, but not liver disease, growth restriction, or lipodystrophy, suggesting distinct tissue-specific pathogenetic mechanisms. Assessment of reactive oxygen species (ROS) production and damage suggested that ROS were not instrumental in the rescue. Cardiac mitochondrial ultrastructure, mitochondrial respiration, and pathological transcriptome and metabolome alterations were essentially normalized by AOX, showing that the restored electron flow upstream of cIII was sufficient to prevent cardiac energetic crisis and detrimental decompensation. These findings demonstrate the value of AOX, both as a mechanistic tool and a potential therapeutic strategy, for cIII deficiencies.

Bracketing phenogenotypic limits of mammalian hybridization.

An increasing number of mammalian species have been shown to have a history of hybridization and introgression based on genetic analyses. Only relatively few fossils, however, preserve genetic material, and morphology must be used to identify the species and determine whether morphologically intermediate fossils could represent hybrids. Because dental and cranial fossils are typically the key body parts studied in mammalian palaeontology, here we bracket the potential for phenotypically extreme hybridizations by examining uniquely preserved cranio-dental material of a captive hybrid between grey and ringed seals. We analysed how distinct these species are genetically and morphologically, how easy it is to identify the hybrids using morphology and whether comparable hybridizations happen in the wild. We show that the genetic distance between these species is more than twice the modern human-Neanderthal distance, but still within that of morphologically similar species pairs known to hybridize. By contrast, morphological and developmental analyses show grey and ringed seals to be highly disparate, and that the hybrid is a predictable intermediate. Genetic analyses of the parent populations reveal introgression in the wild, suggesting that grey-ringed seal hybridization is not limited to captivity. Taken together, we postulate that there is considerable potential for mammalian hybridization between phenotypically disparate taxa.