Corrigendum: Extension of mitogenome enrichment based on single long-range PCR: mtDNAs and putative mitochondrial-derived peptides of five rodent hibernators.

[This corrects the article DOI: 10.3389/fgene.2021.685806.].

Polarity reversal of canine intestinal organoids reduces proliferation and increases cell death.

Apical-out intestinal organoids are a relatively simple method of gaining access to the apical cell surface and have faced increasing scientific interest over the last few years. Apical-out organoids can thus be used for disease modelling to compare differing effects on the basolateral versus the apical cell surface. However, these 'inside-out' organoids die relatively quickly and cannot be propagated as long as their basal-out counterparts. Here, we show that apical-out organoids have drastically reduced proliferative potential, as evidenced by immunohistochemical staining and the incorporation of the thymidine analogue EdU. At the same time, cell death levels are increased. Nevertheless, these phenomena cannot be explained by an induction of differentiation, as the gene expression of key marker genes for various cell types does not change over time.

Mitochondrial polymorphism m.3017C>T of SHLP6 relates to heterothermy.

Heterothermic thermoregulation requires intricate regulation of metabolic rate and activation of pro-survival factors. Eliciting these responses and coordinating the necessary energy shifts likely involves retrograde signalling by mitochondrial-derived peptides (MDPs). Members of the group were suggested before to play a role in heterothermic physiology, a key component of hibernation and daily torpor. Here we studied the mitochondrial single-nucleotide polymorphism (SNP) m.3017C>T that resides in the evolutionarily conserved gene MT-SHLP6. The substitution occurring in several mammalian orders causes truncation of SHLP6 peptide size from twenty to nine amino acids. Public mass spectrometric (MS) data of human SHLP6 indicated a canonical size of 20 amino acids, but not the use of alternative translation initiation codons that would expand the peptide. The shorter isoform of SHLP6 was found in heterothermic rodents at higher frequency compared to homeothermic rodents (p < 0.001). In heterothermic mammals it was associated with lower minimal body temperature (T b, p < 0.001). In the thirteen-lined ground squirrel, brown adipose tissue-a key organ required for hibernation, showed dynamic changes of the steady-state transcript level of mt-Shlp6. The level was significantly higher before hibernation and during interbout arousal and lower during torpor and after hibernation. Our finding argues to further explore the mode of action of SHLP6 size isoforms with respect to mammalian thermoregulation and possibly mitochondrial retrograde signalling.

Corrigendum: Identification of extremely GC-rich micro RNAs for RT-qPCR data normalization in human plasma.

[This corrects the article DOI: 10.3389/fgene.2022.1058668.].

Development of a three-colour digital PCR for early and quantitative detection of benzimidazole resistance-associated single nucleotide polymorphisms in Haemonchus contortus.

Haemonchus contortus is the most pathogenic nematode in small ruminants and anthelmintic resistance (AR) hampers its efficient control. Early detection of AR status is required to reduce selection for AR and cannot be achieved using phenotypic tests. For benzimidazoles (BZs), the detection of AR-associated alleles characterised by single nucleotide polymorphisms (SNPs) in the isotype 1 ß-tubulin gene allows early AR detection in strongyles. The F200Y, F167Y, E198A and E198L polymorphisms have been described in BZ-resistant populations with a clear variation in frequencies between regions. A novel digital PCR (dPCR) enables the detection of all of the above-described polymorphisms in H. contortus. Assays were validated using synthetic DNA fragments containing these SNPs. Then, larvae obtained and pooled at farm level from 26 Austrian and 10 Italian sheep farms were analysed. For all assays a detection limit of 15 copies/µl of resistance alleles and a high level of accuracy were demonstrated, allowing to detect allele frequencies of 1% in most samples. In Austrian samples, elevated frequencies of F200Y resistance alleles were detected on all farms. Polymorphisms in codon 167 and codon 198 were identified in H. contortus from Austria for the first time. In Italian samples, the frequency of resistance alleles was still comparatively low, but F200Y resistance alleles were traceable. In conclusion we developed for the first time dPCR assays that target all SNPs of relevance associated with BZ-resistance in H. contortus. Future research on AR development could benefit from an early onset of SNP-based surveillance that would include the developed assays for all SNPs of relevance. Improved surveillance in the long term will include other important, though less pathogenic, nematode genera in the analyses.

Monoterpene synthases of three closely related sage species (Salvia officinalis, S. fruticosa and S. pomifera, Lamiaceae).

The diversity of plant monoterpenes is largely based on the catalytic activity of monoterpene synthases. Additionally, copy number variation of monoterpene synthase genes may contribute to the quantity of transcripts and hence to the essential oil profile. This study used whole-genome sequencing and digital PCR for the measurement of copy number variation and quantification of gene expression in three closely related Salvia species, namely Salvia officinalis, Salvia pomifera and Salvia fruticosa. Twelve, 13 and 15 monoterpene synthase-encoding open-reading frames were predicted for Salvia officinalis, Salvia pomifera and Salvia fruticosa, respectively. In Salvia officinalis, one of the open reading frames was disrupted indicating a pseudogene. Monoterpene synthase genes were generally single copy per haploid genome, only a few were double or triple copy genes. Expression levels of monoterpene synthases in leaves corresponded generally well with essential oil composition. In some cases, a higher expression level of a certain monoterpene synthase could be explained by its duplication or triplication. The very high content of thujones in Salvia pomifera, for example, was accompanied by gene duplication and increased gene expression of (+)-sabinene synthase responsible for the thujone precursor sabinene. In Salvia officinalis, three individuals different in their essential oil profile showed significant differences in their monoterpene synthase expression levels corresponding roughly to the profile of the essential oils. Transcript expression of monoterpene synthase genes were measured in leaf, calyx and corolla. The corolla differed significantly from leaves, while calyces usually showed a profile intermediary between leaf and corolla.

Identification of extremely GC-rich micro RNAs for RT-qPCR data normalization in human plasma.

We aimed at extending the repertoire of high-quality miRNA normalizers for reverse transcription-quantitative PCR (RT-qPCR) of human plasma with special emphasis on the extremely guanine-cytosine-rich portion of the miRNome. For high-throughput selection of stable candidates, microarray technology was preferred over small-RNA sequencing (sRNA-seq) since the latter underrepresented miRNAs with a guanine-cytosine (GC) content of at least 75% (p = 0.0002, n = 2). miRNA abundances measured on the microarray were ranked for consistency and uniformity using nine normalization approaches. The eleven most stable sequences included miRNAs of moderate, but also extreme GC content (45%-65%: miR-320d, miR-425-5p, miR-185-5p, miR-486-5p; 80%-95%: miR-1915-3p, miR-3656-5p, miR-3665-5p, miR-3960-5p, miR-4488-5p, miR-4497 and miR-4787-5p). In contrast, the seven extremely GC-rich miRNAs were not found in the two plasma miRNomes screened by sRNA-seq. Stem-loop RT-qPCR was employed for stability verification in 32 plasma samples of healthy male Caucasians (age range: 18-55 years). In general, inter-individual variance of miRNA abundance was low or very low as indicated by coefficient of variation (CV) values of 0.6%-8.2%. miR-3665 and miR-1915-3p outperformed in this analysis (CVs: 0.6 and 2.4%, respectively). The eight most stable sequences included four extremely GC-rich miRNAs (miR-1915-3p, miR-3665, miR-4787-5p and miR-4497). The best-performing duo normalization factor (NF) for the condition of human plasma, miR-320d and miR-4787-5p, also included a GC-extreme miRNA. In summary, the identification of extremely guanine-cytosine-rich plasma normalizers will help to increase accuracy of PCR-based miRNA quantification, thus raise the potential that miRNAs become markers for psychological stress reactions or early and precise diagnosis of clinical phenotypes. The novel miRNAs might also be useful for orthologous contexts considering their conservation in related animal genomes.

A Recurrent STAT5BN642H Driver Mutation in Feline Alimentary T Cell Lymphoma.

Alimentary lymphomas arising from T cells are rare and aggressive malignancies in humans. In comparison, they represent the most common anatomical form of lymphoma in cats. Due to the low prevalence in humans, the underlying pathomechanism for these diseases is poorly characterised, limiting experimental analysis and therapeutic exploration. To date, activating mutations of the JAK/STAT core cancer pathway and particularly the STAT5B oncoprotein have been identified in human enteropathy-associated T cell lymphoma. Here, we describe a high homology of human and feline STAT3 and STAT5B proteins and strong conservation at the genomic level. Analysis of 42 samples of feline T cell alimentary lymphoma reveals broad activation of STAT3 and STAT5B. Screening for known activating mutations in STAT3 or STAT5B identifies the presence of the STAT5BN642H driver mutation in feline enteropathy-associated T cell lymphoma in 7 out of 42 (16.67%) samples in total. Regarding lymphoma subtypes, the majority of mutations with 5 out of 17 (29.41%) cases were found in feline enteropathy-associated lymphoma type II (EATL II). This identification of an oncogenic STAT5B driver mutation in felines recapitulates the genetic situation in the corresponding human disease, thereby establishing the cat as a potential new model for a rare and incurable human T cell disease.

Equine Parvovirus-Hepatitis Screening in Horses and Donkeys with Histopathologic Liver Abnormalities.

There is strong evidence that equine parvovirus-hepatitis (EqPV-H) is associated with the onset of Theiler's disease, an acute hepatic necrosis, in horses. However, the impact of this virus on other hepatopathies remains unknown. The objective of this retrospective study was to evaluate the prevalence and quantify the viral loads of EqPV-H in formalin-fixed, paraffin-embedded equine and donkey livers with various histopathologic abnormalities. The pathologies included cirrhosis, circulatory disorders of the liver, toxic and metabolic hepatic diseases as well as neoplastic and inflammatory diseases (n = 84). Eight normal liver samples were included for comparison as controls. EqPV-H DNA was qualitatively and quantitatively measured by real-time PCR and digital PCR, respectively. The virus was detected in two livers originating from horses diagnosed with abdominal neoplasia and liver metastasis (loads of 5 × 103 and 9.5 × 103 genome equivalents per million cells). The amount of viral nucleic acids measured indicates chronic infection or persistence of EqPV-H, which might have been facilitated by the neoplastic disease. In summary, this study did not provide evidence for EqPV-H being involved in hepatopathies other than Theiler's disease.

Cytoplasmic Transfer Methods for Studying the Segregation of Mitochondrial DNA in Mice.

Heteroplasmic mice represent a valuable tool to study the segregation of different mtDNA haplotypes (mtDNAs with differing alleles) in vivo against a defined nuclear background. We describe two methods for the creation of such models, differing in the resulting initial heteroplasmy levels: (a) transfer of ooplasm and (b) fusion of two blastomeres. These methods result in typical heteroplasmy of 5% and 50% donor mtDNA , respectively. The choice of method depends on the aim of the study. By means of breeding even 100% donor mtDNA can be reached within a few generations.

Extension of Mitogenome Enrichment Based on Single Long-Range PCR: mtDNAs and Putative Mitochondrial-Derived Peptides of Five Rodent Hibernators.

Enriching mitochondrial DNA (mtDNA) for sequencing entire mitochondrial genomes (mitogenomes) can be achieved by single long-range PCR. This avoids interference from the omnipresent nuclear mtDNA sequences (NUMTs). The approach is currently restricted to the use of samples collected from humans and ray-finned fishes. Here, we extended the use of single long-range PCR by introducing back-to-back oligonucleotides that target a sequence of extraordinary homology across vertebrates. The assay was applied to five hibernating rodents, namely alpine marmot, Arctic and European ground squirrels, and common and garden dormice, four of which have not been fully sequenced before. Analysis of the novel mitogenomes focussed on the prediction of mitochondrial-derived peptides (MDPs) providing another level of information encoded by mtDNA. The comparison of MOTS-c, SHLP4 and SHLP6 sequences across vertebrate species identified segments of high homology that argue for future experimentation. In addition, we evaluated four candidate polymorphisms replacing an amino acid in mitochondrially encoded subunits of the oxidative phosphorylation (OXPHOS) system that were reported in relation to cold-adaptation. No obvious pattern was found for the diverse sets of mammalian species that either apply daily or multiday torpor or otherwise cope with cold. In summary, our single long-range PCR assay applying a pair of back-to-back primers that target a consensus sequence motif of Vertebrata has potential to amplify (intact) mitochondrial rings present in templates from a taxonomically diverse range of vertebrates. It could be promising for studying novel mitogenomes, mitotypes of a population and mitochondrial heteroplasmy in a sensitive, straightforward and flexible manner.

S100A4 mRNA-protein relationship uncovered by measurement noise reduction.

Intrinsic biological fluctuation and/or measurement error can obscure the association of gene expression patterns between RNA and protein levels. Appropriate normalization of reverse-transcription quantitative PCR (RT-qPCR) data can reduce technical noise in transcript measurement, thus uncovering such relationships. The accuracy of gene expression measurement is often challenged in the context of cancer due to the genetic instability and "splicing weakness" involved. Here, we sequenced the poly(A) cancer transcriptome of canine osteosarcoma using mRNA-Seq. Expressed sequences were resolved at the level of two consecutive exons to enable the design of exon-border spanning RT-qPCR assays and ranked for stability based on the coefficient of variation (CV). Using the same template type for RT-qPCR validation, i.e. poly(A) RNA, avoided skewing of stability assessment by circular RNAs (circRNAs) and/or rRNA deregulation. The strength of the relationship between mRNA expression of the tumour marker S100A4 and its proportion score of quantitative immunohistochemistry (qIHC) was introduced as an experimental readout to fine-tune the normalization choice. Together with the essential logit transformation of qIHC scores, this approach reduced the noise of measurement as demonstrated by uncovering a highly significant, strong association between mRNA and protein expressions of S100A4 (Spearman's coefficient ρ = 0.72 (p = 0.006)). KEY MESSAGES: • RNA-seq identifies stable pairs of consecutive exons in a heterogeneous tumour. • Poly(A) RNA templates for RT-qPCR avoid bias from circRNA and rRNA deregulation. • HNRNPL is stably expressed across various cancer tissues and osteosarcoma. • Logit transformed qIHC score better associates with mRNA amount. • Quantification of minor S100A4 mRNA species requires poly(A) RNA templates and dPCR.

Protein measurements in venous plasma, earlobe capillary plasma and in plasma stored on filter paper.

In this study, levels of inflammatory protein biomarkers in venous plasma, plasma derived from capillary blood from the earlobe, and capillary plasma stored as dried plasma spots (DPS) were compared. Samples from 12 male individuals were assessed with a panel of 92 inflammation-related proteins using multiplex proximity extension assay. Correlations between sample types varied greatly between analytes. A high correlation of ρ > 0.8 was observed between capillary plasma and DPS for 32 analytes. At this level of correlation, 13 analytes correlated between venous and capillary plasma and 5 analytes in the comparison of venous blood with DPS.

Oxygen Tension Strongly Influences Metabolic Parameters and the Release of Interleukin-6 of Human Amniotic Mesenchymal Stromal Cells In Vitro.

The human amniotic membrane (hAM) has been used for tissue regeneration for over a century. In vivo (in utero), cells of the hAM are exposed to low oxygen tension (1-4% oxygen), while the hAM is usually cultured in atmospheric, meaning high, oxygen tension (20% oxygen). We tested the influence of oxygen tensions on mitochondrial and inflammatory parameters of human amniotic mesenchymal stromal cells (hAMSCs). Freshly isolated hAMSCs were incubated for 4 days at 5% and 20% oxygen. We found 20% oxygen to strongly increase mitochondrial oxidative phosphorylation, especially in placental amniotic cells. Oxygen tension did not impact levels of reactive oxygen species (ROS); however, placental amniotic cells showed lower levels of ROS, independent of oxygen tension. In contrast, the release of nitric oxide was independent of the amniotic region but dependent on oxygen tension. Furthermore, IL-6 was significantly increased at 20% oxygen. To conclude, short-time cultivation at 20% oxygen of freshly isolated hAMSCs induced significant changes in mitochondrial function and release of IL-6. Depending on the therapeutic purpose, cultivation conditions of the cells should be chosen carefully for providing the best possible quality of cell therapy.

Large-scale genetic analysis reveals mammalian mtDNA heteroplasmy dynamics and variance increase through lifetimes and generations.

Vital mitochondrial DNA (mtDNA) populations exist in cells and may consist of heteroplasmic mixtures of mtDNA types. The evolution of these heteroplasmic populations through development, ageing, and generations is central to genetic diseases, but is poorly understood in mammals. Here we dissect these population dynamics using a dataset of unprecedented size and temporal span, comprising 1947 single-cell oocyte and 899 somatic measurements of heteroplasmy change throughout lifetimes and generations in two genetically distinct mouse models. We provide a novel and detailed quantitative characterisation of the linear increase in heteroplasmy variance throughout mammalian life courses in oocytes and pups. We find that differences in mean heteroplasmy are induced between generations, and the heteroplasmy of germline and somatic precursors diverge early in development, with a haplotype-specific direction of segregation. We develop stochastic theory predicting the implications of these dynamics for ageing and disease manifestation and discuss its application to human mtDNA dynamics.

Cellular and Site-Specific Mitochondrial Characterization of Vital Human Amniotic Membrane.

Over a century ago, clinicians started to use the human amniotic membrane for coverage of wounds and burn injuries. To date, literally thousands of different clinical applications exist for this biomaterial almost exclusively in a decellularized or denuded form. Recent reconsiderations for the use of vital human amniotic membrane for clinical applications would take advantage of the versatile cells of embryonic origin including the entirety of their cell organelles. Recently, more and more evidence was found, showing mitochondria to be involved in most fundamental cellular processes, such as differentiation and cell death. In this study, we focused on specific properties of mitochondria of vital human amniotic membrane and characterized bioenergetical parameters of 2 subregions of the human amniotic membrane, the placental and reflected amnion. We found significantly different levels of adenosine triphosphate (ATP) and extracellular reactive oxygen species, concentrations of succinate dehydrogenase, and lactate upon inhibition of ATP synthase in placental and reflected amnion. We also found significantly different rates of mitochondrial respiration in isolated human amniotic epithelial cells and human amniotic mesenchymal stromal cells, according to the subregions. Differences in metabolic activities were inversely related to mitochondrial DNA copy numbers in isolated cells of placental and reflected amnion. Based on significant differences of several key parameters of energy metabolism in 2 subregions of vital amnion, we propose that these metabolic differences of vital placental and reflected amnion could have critical impact on therapeutic applications. Inclusion of region-specific metabolic properties could optimize and fine-tune the clinical application of the human amniotic membrane and improve the outcome significantly.

Retinoic acid prevents immunogenicity of milk lipocalin Bos d 5 through binding to its immunodominant T-cell epitope.

The major cow's milk allergen Bos d 5 belongs to the lipocalin protein family, with an intramolecular pocket for hydrophobic ligands. We investigated whether Bos d 5 when loaded with the active vitamin A metabolite retinoic acid (RA), would elicit differential immune responses compared to the unloaded state. By in silico docking an affinity energy of -7.8 kcal/mol was calculated for RA into Bos d 5. Loading of RA to Bos d 5 could be achieved in vitro, as demonstrated by ANS displacement assay, but had no effect on serum IgE binding in tolerant or challenge-positive milk allergic children. Bioinformatic analysis revealed that RA binds to the immunodominant T-cell epitope region of Bos d 5. In accordance, Bos d 5 significantly suppressed the CD3+ CD4+ cell numbers, proliferative response and IL-10, IL-13 and IFN-γ secretion from stimulated human PBMCs only when complexed with RA. This phenomenon was neither associated with apoptosis of T-cells nor with the activation of Foxp3+ T-cells, but correlated likely with enhanced stability to lysosomal digestion due to a predicted overlap of Cathepsin S cleavage sites with the RA binding site. Taken together, proper loading of Bos d 5 with RA may suppress its immunogenicity and prevent its allergenicity.

Stem cell growth factor receptor in canine vs. feline osteosarcomas.

Osteosarcoma is considered the most common bone cancer in cats and dogs, with cats having a much better prognosis than dogs, since the great majority of dogs with osteosarcoma develop distant metastases. In search of a factor possibly contributing to this disparity, the stem cell growth factor receptor KIT was targeted, and the messenger (m)RNA and protein expression levels of KIT were compared in canine vs. feline osteosarcomas, as well as in normal bone. The mRNA expression of KIT was quantified by reverse transcription-quantitative polymerase chain reaction, and was observed to be significantly higher in canine (n=14) than in feline (n=5) osteosarcoma samples (P<0.001). KIT protein expression was evaluated by immunohistochemistry, which revealed that 21% of canine osteosarcoma samples did not exhibit KIT staining in their neoplastic cells, while in 14% of samples, a score of 1 (<10% positive tumour cells) was observed, and in 50% and 14% of samples, a score of 2 (10-50% positivity) and 3 (>50% positivity), respectively, was observed. By contrast, the cancer cells of all the feline bone tumour samples analysed were entirely negative for KIT. Notably, canine and feline osteocytes of healthy bone tissue lacked any KIT expression. These results could be the first evidence that KIT may be involved in the higher aggressiveness of canine osteosarcoma compared with feline osteosarcoma.

Allelic variation of the COMT gene in a despotic primate society: A haplotype is related to cortisol excretion in Macaca fuscata.

Sequence variations in genes of the monoamine neurotransmitter system and their common function in human and non-human primate species are an ongoing issue of investigation. However, the COMT gene, coding for the catechol-O-methyltransferase, has not yet attracted much scientific attention regarding its functional role in non-human primates. Considering that a polymorphism of the human COMT gene affects the enzyme activity and cortisol level in response to a social stressor, this study investigated the impact of COMT on endocrine stress and behavioural parameters in Japanese macaques (Macaca fuscata). The species exemplifies a despotic hierarchy in which males' social rank positions require an adaptation of behaviour strategies. During the mating period steroid secretion and the frequency of aggressive encounters between males increase. We addressed i) whether this species exhibits potential functional COMT variants, ii) whether these variants are associated with faecal cortisol excretion of males, iii) how they are distributed among different social rank positions and iv) whether they are associated with behavioural strategies during times of mate competition. By genotyping 26 males we identified three COMT haplotypes (HT), including a putative splice mutant (HT3). This variant was associated with increased cortisol excretion. Given the observed inverse correlation between cortisol and physical aggression, we assume that different COMT haplotypes may predispose individuals to pursue more or less aggressive strategies. How these gene-stress effects might favour a specific social role is discussed. Our study of non-invasive genotyping in combination with behavioural and endocrine parameters represents an important step towards the understanding of gene-stress effects in a hierarchically organised primate society.

The Orthology Clause in the Next Generation Sequencing Era: Novel Reference Genes Identified by RNA-seq in Humans Improve Normalization of Neonatal Equine Ovary RT-qPCR Data.

BACKGROUND: Vertebrate evolution is accompanied by a substantial conservation of transcriptional programs with more than a third of unique orthologous genes showing constrained levels of expression. Moreover, there are genes and exons exhibiting excellent expression stability according to RNA-seq data across a panel of eighteen tissues including the ovary (Human Body Map 2.0). RESULTS: We hypothesized that orthologs of these exons would also be highly uniformly expressed across neonatal ovaries of the horse, which would render them appropriate reference genes (RGs) for normalization of reverse transcription quantitative PCR (RT-qPCR) data in this context. The expression stability of eleven novel RGs (C1orf43, CHMP2A, EMC7, GPI, PSMB2, PSMB4, RAB7A, REEP5, SNRPD3, VCP and VPS29) was assessed by RT-qPCR in ovaries of seven neonatal fillies and compared to that of the expressed repetitive element ERE-B, two universal (OAZ1 and RPS29) and four traditional RGs (ACTB, GAPDH, UBB and B2M). Expression stability analyzed with the software tool RefFinder top ranked the normalization factor constituted of the genes SNRPD3 and VCP, a gene pair that is not co-expressed according to COEXPRESdb and GeneMANIA. The traditional RGs GAPDH, B2M, ACTB and UBB were only ranked 3rd and 12th to 14th, respectively. CONCLUSIONS: The functional diversity of the novel RGs likely facilitates expression studies over a wide range of physiological and pathological contexts related to the neonatal equine ovary. In addition, this study augments the potential for RT-qPCR-based profiling of human samples by introducing seven new human RG assays (C1orf43, CHMP2A, EMC7, GPI, RAB7A, VPS29 and UBB).