Out with the old, in with the new: Meiotic driving of sex chromosome evolution.

Chromosomal regions with meiotic drivers exhibit biased transmission (> 50 %) over their competing homologous chromosomal region. These regions often have two prominent genetic features: suppressed meiotic crossing over and rapidly evolving multicopy gene families. Heteromorphic sex chromosomes (e.g., XY) often share these two genetic features with chromosomal regions exhibiting meiotic drive. Here, we discuss parallels between meiotic drive and sex chromosome evolution, how the divergence of heteromorphic sex chromosomes can be influenced by meiotic drive, experimental approaches to study meiotic drive on sex chromosomes, and meiotic drive in traditional and non-traditional model organisms with high-quality genome assemblies. The newly available diversity of high-quality sex chromosome sequences allows us to revisit conventional models of sex chromosome evolution through the lens of meiotic drive.

Expression of Multicopy Tandem Recombinant Ginseng Hexapeptide in Bacillus subtilis and the Evaluation of Antiaging Activity.

Ginseng oligopeptides are naturally occurring small-molecule peptides extracted from ginseng that exhibit positive effects on health and longevity. However, the current industrial production of ginseng oligopeptides primarily relies on plant extraction and chemical synthesis. In this study, we proposed a novel genetic engineering approach to produce active ginseng peptides through multicopy tandem insertion (5 and 15 times). The recombinant ginseng peptides were successfully produced from engineered Bacillus subtilis with an increasing yield from 356.55 to 2900 mg/L as the repeats multiple. Additionally, an oxidative stress-induced aging model caused by H2O2 was established to evaluate whether the recombinant ginseng peptides, without enzymatic hydrolysis into individual peptides, also have positive effects on antiaging. The results demonstrated that all two kinds of recombinant ginseng peptides could also delay cellular aging through various mechanisms, such as inhibiting cell cycle arrest, suppressing the expression of pro-inflammatory factors, and enhancing cellular antioxidant capacity.

High-Level De Novo Production of (2S)-Eriodictyol in Yarrowia Lipolytica by Metabolic Pathway and NADPH Regeneration Engineering.

(2S)-Eriodictyol, a polyphenolic flavonoid, has found widespread applications in health supplements and food additives. However, the limited availability of plant-derived (2S)-eriodictyol cannot meet the market demand. Microbial production of (2S)-eriodictyol faces challenges, including the low catalytic efficiency of flavone 3'-hydroxylase/cytochrome P450 reductase (F3'H/CPR), insufficient precursor supplementation, and inadequate NADPH regeneration. This study systematically engineered Yarrowia lipolytica for high-level (2S)-eriodictyol production. In doing this, the expression of F3'H/CPR was balanced, and the supply of precursors was enhanced by relieving feedback inhibition of the shikimate pathway, promoting fatty acid ß-oxidation, and increasing the copy number of synthetic pathway genes. These strategies, combined with NADPH regeneration, achieved an (2S)-eriodictyol titer of 423.6 mg/L. Finally, in fed-batch fermentation, a remarkable 6.8 g/L (2S)-eriodictyol was obtained, representing the highest de novo microbial titer reported to date and paving the way for industrial production.

YaliCMulti and YaliHMulti: Stable, efficient multi-copy integration tools for engineering Yarrowia lipolytica.

Yarrowia lipolytica is widely used in biotechnology to produce recombinant proteins, food ingredients and diverse natural products. However, unstable expression of plasmids, difficult and time-consuming integration of single and low-copy-number plasmids hampers the construction of efficient production pathways and application to industrial production. Here, by exploiting sequence diversity in the long terminal repeats (LTRs) of retrotransposons and ribosomal DNA (rDNA) sequences, a set of vectors and methods that can recycle multiple and high-copy-number plasmids was developed that can achieve stable integration of long-pathway genes in Y. lipolytica. By combining these sequences, amino acids and antibiotic tags with the Cre-LoxP system, a series of multi-copy site integration recyclable vectors were constructed and assessed using the green fluorescent protein (HrGFP) reporter system. Furthermore, by combining the consensus sequence with the vector backbone of a rapidly degrading selective marker and a weak promoter, multiple integrated high-copy-number vectors were obtained and high levels of stable HrGFP expression were achieved. To validate the universality of the tools, simple integration of essential biosynthesis modules was explored, and 7.3 g/L of L-ergothioneine and 8.3 g/L of (2S)-naringenin were achieved in a 5 L fermenter, the highest titres reported to date for Y. lipolytica. These novel multi-copy genome integration strategies provide convenient and effective tools for further metabolic engineering of Y. lipolytica.

Extensively drug-resistant Enterobacter ludwigii co-harbouring MCR-9 and a multicopy of blaIMP-1 in South Korea.

In this study, we describe an Enterobacter ludwigii clinical isolate that is resistant to both carbapenems and colistin in South Korea. Antimicrobial susceptibility testing revealed that E. ludwigii CRE2104-31 was non-susceptible to all tested antibiotics except fosfomycin. Whole genome sequencing identified a 323-kbp IncHI2 plasmid, pCRE2104-31a, that was co-harbouring mobile colistin resistance (mcr)-9.1 and blaIMP-1. In comparison with other full plasmids, pCRE2104-31a exhibited the closest similarity to a plasmid from the Klebsiella pneumoniae strain CNR48 from France, with 19.9% query coverage and 99% identity. Notably, we observed five tandem repeats of blaIMP-1 and aac(6')-Il genes, accompanied by multiple attCs within a class I integron on the Tn402-like transposon. The unit of blaIMP-1-attC-aac(6')-Il-attC might have accumulated due to multiple convergent events. In addition to mcr-9.1 and blaIMP-1, various other antibiotic resistance-associated genes were identified in the plasmid, as follows: blaTEM-1B, aph(3')-I, aph(3')-Ia, aac(6')-Il, aac(6')-IIc, aac(6')-IIa, aph(6)-Id, aph(3'')-Ib, aadA2b, aac(6')-Ib3, sul, dfrA19, qnrB2, aac(6')-Ib-cr, ere(A), and qacE. A conjugation assay showed that the mcr-9.1/blaIMP-1-co-bearing plasmid was self-transmissible to E. coli J53. However, colistin and carbapenem resistance could not be transferred to E. coli due to high incompatibility. The convergence of mcr and carbapenemase genes is thought to be host-dependent among Enterobacteriaceae. The emergence of extensively drug-resistant E. ludwigii co-harbouring MCR-9.1 and a multicopy of blaIMP-1 would pose a significant threat within the compatible Enterobacteriaceae.

Specific detection of Phytophthora parasitica by recombinase polymerase amplification (RPA) assays based on a unique multi-copy genomic sequence.

Phytophthora parasitica is a highly destructive oomycete plant pathogen that is capable of infecting a wide range of hosts including many agricultural cash crops, fruit trees, and ornamental garden plants. One of the most important diseases caused by P. parasitica worldwide is black shank of tobacco. Rapid, sensitive, and specific pathogen detection is crucial for early rapid diagnosis which can facilitate effective disease management. In this study, we used a genomics approach to identify repeated sequences in the genome of P. parasitica by genome sequence alignment, and identified a 203 bp P. parasitica-specific sequence, PpM34, that is present in 31-60 copies in the genome. The P. parasitica genome-specificity of PpM34 was supported by PCR amplification of 24 genetically diverse strains of P. parasitica, 32 strains representing twelve other Phytophthora species, one Pythium specie, six fungal species and three bacterial species, all of which are plant pathogens. Our PCR and real-time PCR assays showed that the PpM34 sequence was highly sensitive in specifically detecting P. parasitica. Finally, we developed a PpM34-based high-efficiency Recombinase Polymerase Amplification (RPA) assay, which allowed us to specifically detect as little as 1 pg of P. parasitica total DNA from both pure cultures and infected Nicotiana benthamiana at 39°C using a fluorometric thermal cycler. The sensitivity, specificity, convenience and rapidity of this assay represents a major improvement for early diagnosis of P. parasitica infection.

Novel genetic tools improve Penicillium expansum patulin synthase production in Aspergillus niger.

Since the first CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) system was developed for creating double-stranded DNA breaks, it has been adapted and improved for different biotechnological applications. In this issue of The FEBS Journal, Arentshorst et al. developed a novel approach to enhance transgene expression of a specific protein, patulin synthase (PatE) from Penicillium expansum, in the important industrial filamentous fungus Aspergillus niger. Their technique involved the disruption of selected genes with counter-effects on targeted protein production and simultaneous integration of glucoamylase landing sites into the disrupted gene locus such as protease regulator (prtT) in an ATP-dependent DNA helicase II subunit 1 (kusA or ku70)-deletion strain. Multiple copies of the PatE transgene expression cassette were introduced by CRISPR-Cas9-mediated insertion. The purified PatE was further used for structural and functional studies, and the technique laid the foundation for elevating the overall production of various proteins or chemicals in those industrially important fungi.

Alternative PCR-Based Approaches for Generation of Komagataella phaffii Strains.

Komagataella phaffii (Pichia pastoris) is a widely known microbial host for recombinant protein production and an emerging model organism in fundamental research. The development of new materials and techniques on this yeast improves heterologous protein synthesis. One of the most prominent ways to enhance protein production efficiency is to select K. phaffii strains with multiple expression cassettes and generate MutS strains using various vectors. In this study, we demonstrate approaches to expand the applications of pPICZ series vectors. Procedures based on PCR amplification and in vivo cloning allow rapid exchange of selectable markers. The combination of PCR amplification with split-marker-mediated transformation allows the development of K. phaffii MutS strains with two expression cassettes using pPICZ vectors. Both PCR-based approaches were applied to efficiently produce interleukin-2 mimetic Neo-2/15 in K. phaffii. The described techniques provide alternative ways to generate and improve K. phaffii strains without the need for obtaining new specific vectors or additional cloning of expression cassettes.

Hijacking a small plasmid to confer high-level resistance to aztreonam-avibactam and ceftazidime-avibactam.

Acquired ß-lactamase-encoding genes are typically carried by large plasmids in Gram-negative bacteria, which also commonly carry multi-copy small plasmids. This study found that mobile genetic elements carrying antimicrobial resistance genes are capable of hijacking small plasmids. This study focused on aztreonam-avibactam (ATM-AVI) as this combination can be used to effectively counter almost all ß-lactamases produced by bacteria, and has been recommended against carbapenem-resistant Enterobacterales. A clinical strain (085003) of carbapenem-resistant Escherichia coli was investigated, and mutants (085003R32 and 085003R512) able to grow under 32/4 and 512/4 mg/L of ATM-AVI were obtained as representatives of low- and high-level resistance, respectively, by induction. Comparative genomics showed that 085003R32 and 085003R512 had a single nucleotide mutation of ß-lactamase gene blaCMY-2, encoding a novel CMY with a Thr319Ile substitution, assigned 'CMY-2R'. Cloning and enzyme kinetics were used to verify that CMY-2R conferred ATM-AVI resistance by compromising binding of AVI and subsequent protection of ATM. Mechanisms for the discrepant resistance between 085003R32 and 085003R512 were investigated. Three tandem copies of blaCMY-2R were identified on a self-transmissible IncP1 plasmid of 085003R32 due to IS1294 misrecognizing its end terIS and rolling-circle replication. 085003R512 had only a single copy of blaCMY-2R on the IncP1 plasmid, but possessed anther blaCMY-2R on an already present 4-kb small plasmid. IS1294-mediated mobilization on to this multi-copy small plasmid increased the copy number of blaCMY-2R significantly, rendering higher resistance. This study shows that bacteria can employ multiple approaches to accommodate selection pressures imposed by exposure to varied concentrations of antimicrobial agents.

Nanopore sequencing of PCR products enables multicopy gene family reconstruction.

The importance of gene amplifications in evolution is more and more recognized. Yet, tools to study multi-copy gene families are still scarce, and many such families are overlooked using common sequencing methods. Haplotype reconstruction is even harder for polymorphic multi-copy gene families. Here, we show that all variants (or haplotypes) of a multi-copy gene family present in a single genome, can be obtained using Oxford Nanopore Technologies sequencing of PCR products, followed by steps of mapping, SNP calling and haplotyping. As a proof of concept, we acquired the sequences of highly similar variants of the cidA and cidB genes present in the genome of the Wolbachia wPip, a bacterium infecting Culex pipiens mosquitoes. Our method relies on a wide database of cid genes, previously acquired by cloning and Sanger sequencing. We addressed problems commonly faced when using mapping approaches for multi-copy gene families with highly similar variants. In addition, we confirmed that PCR amplification causes frequent chimeras which have to be carefully considered when working on families of recombinant genes. We tested the robustness of the method using a combination of bioinformatics (read simulations) and molecular biology approaches (sequence acquisitions through cloning and Sanger sequencing, specific PCRs and digital droplet PCR). When different haplotypes present within a single genome cannot be reconstructed from short reads sequencing, this pipeline confers a high throughput acquisition, gives reliable results as well as insights of the relative copy numbers of the different variants.

A landing pad system for multicopy gene integration in Issatchenkia orientalis.

The robust nature of the non-conventional yeast Issatchenkia orientalis allows it to grow under highly acidic conditions and therefore, has gained increasing interest in producing organic acids using a variety of carbon sources. Recently, the development of a genetic toolbox for I. orientalis, including an episomal plasmid, characterization of multiple promoters and terminators, and CRISPR-Cas9 tools, has eased the metabolic engineering efforts in I. orientalis. However, multiplex engineering is still hampered by the lack of efficient multicopy integration tools. To facilitate the construction of large, complex metabolic pathways by multiplex CRISPR-Cas9-mediated genome editing, we developed a bioinformatics pipeline to identify and prioritize genome-wide intergenic loci and characterized 47 gRNAs located in 21 intergenic regions. These loci are screened for guide RNA cutting efficiency, integration efficiency of a gene cassette, the resulting cellular fitness, and GFP expression level. We further developed a landing pad system using components from these well-characterized loci, which can aid in the integration of multiple genes using single guide RNA and multiple repair templates of the user's choice. We have demonstrated the use of the landing pad for simultaneous integrations of 2, 3, 4, or 5 genes to the target loci with efficiencies greater than 80%. As a proof of concept, we showed how the production of 5-aminolevulinic acid can be improved by integrating five copies of genes at multiple sites in one step. We have further demonstrated the efficiency of this tool by constructing a metabolic pathway for succinic acid production by integrating five gene expression cassettes using a single guide RNA along with five different repair templates, leading to the production of 9 g/L of succinic acid in batch fermentations. This study demonstrates the effectiveness of a single gRNA-mediated CRISPR platform to build complex metabolic pathways in a non-conventional yeast. This landing pad system will be a valuable tool for the metabolic engineering of I. orientalis.

Improving solubility and copy number of taxadiene synthase to enhance the titer of taxadiene in Yarrowia lipolytica.

Taxadiene is an important precursor for the biosynthesis of highly effective anticancer drug paclitaxel, but its microbial biosynthesis yield is very low. In this study, we employed Yarrowia lipolytica as a microbial host to produce taxadiene. First, a "push-pull" strategy was adopted to increase taxadiene production by 234%. Then taxadiene synthase was fused with five solubilizing tags respectively, leading a maximum increase of 62.3% in taxadiene production when fused with SUMO. Subsequently, a multi-copy iterative integration method was used to further increase taxadiene titer, achieving the maximum titer of 23.7 mg/L in shake flask culture after three rounds of integration. Finally, the taxadiene titer was increased to 101.4 mg/L by optimization of the fed-batch fermentation conditions. This is the first report of taxadiene biosynthesis accomplished in Y. lipolytica, serving as a good example for the sustainable production of taxadiene and other terpenoids in this oleaginous yeast.

Efficient traceless modification of the P1 bacteriophage genome through homologous recombination with enrichment in double recombinants: A new perspective on the functional annotation of uncharacterized phage genes.

The advent of high-throughput omic technologies has caused unprecedented progress in research on bacteriophages, the most abundant and still the least explored entities on earth. Despite the growing number of phage genomes sequenced and the rejuvenation of interest in phage therapy, the progress in the functional analysis of phage genes is slow. Simple and efficient techniques of phage genome targeted mutagenesis that would allow one to knock out particular genes precisely without polar effects in order to study the effect of these knock-outs on phage functions are lacking. Even in the case of model phages, the functions of approximately half of their genes are unknown. P1 is an enterobacterial temperate myophage of clinical significance, which lysogenizes cells as a plasmid. It has a long history of studies, serves as a model in basic research, is a gene transfer vector, and is a source of genetic tools. Its gene products have structural homologs in several other phages. In this perspective article, we describe a simple and efficient procedure of traceless P1 genome modification that could also serve to acquire targeted mutations in the genomes of certain other temperate phages and speed up functional annotations of phage genes.

Enhanced Expression of Alcohol Dehydrogenase I in Pichia pastoris Reduces the Content of Acetaldehyde in Wines.

Acetaldehyde is an important carbonyl compound commonly detected in wines. A high concentration of acetaldehyde can affect the flavor of wines and result in adverse effects on human health. Alcohol dehydrogenase I (ADH1) in Saccharomyces cerevisiae catalyzes the reduction reaction of acetaldehyde into ethanol in the presence of cofactors, showing the potential to reduce the content of acetaldehyde in wines. In this study, ADH1 was successfully expressed in Pichia pastoris GS115 based on codon optimization. Then, the expression level of ADH1 was enhanced by replacing its promoter with optimized promoters and increasing the copy number of the expression cassette, with ADH1 being purified using nickel column affinity chromatography. The enzymatic activity of purified ADH1 reached 605.44 ± 44.30 U/mg. The results of the effect of ADH1 on the content of acetaldehyde in wine revealed that the acetaldehyde content of wine samples was reduced from 168.05 ± 0.55 to 113.17 ± 6.08 mg/L with the addition of 5 mM NADH and the catalysis of ADH1, and from 135.53 ± 4.08 to 52.89 ± 2.20 mg/L through cofactor regeneration. Our study provides a novel approach to reducing the content of acetaldehyde in wines through enzymatic catalysis.

Heterologous Expression of Extracellular Proteinase pAsPs of Aspergillus pseudotamarii in Komagataella phaffii.

Neutral protease pAsPs gene was obtained by sequence optimization of NpI protease from Aspergillus pseudotamarii. pAsPs was for the first time integrated in the genome of yeast strain Komagataella phaffii T07, and then produced in a 5 L bioreactor with an enzyme yield of 150,800 U/mL of culture liquid towards casein. The specific activity of the pAsPs was 7,657,000 U/mg toward casein, 2320 U/mg toward hemoglobin, and 25,344 U/mg toward azocasein per 1 mg of the protein. The enzyme was found to be inhibited by Cu2+. Optimal activity pH was shown in the range of pH 6.5-8.0, and optimal temperature-50-60 °C. The molecular mass of the recombinant protease pAsPs was shown to be 67.5 kDa. Mass-spectrometric analysis confirmed the identity of the amino acid sequence of the obtained pAsPs preparation with the predicted sequence, with 17% coverage and protein score 288. Thus, the novel neutral protease pAsPs is a promising candidate for large-scale use in manufacturing, including the food industry.

Copy number variation of horse Y chromosome genes in normal equine populations and in horses with abnormal sex development and subfertility: relationship of copy number variations with Y haplogroups.

Structural rearrangements like copy number variations in the male-specific Y chromosome have been associated with male fertility phenotypes in human and mouse but have been sparsely studied in other mammalian species. Here, we designed digital droplet PCR assays for 7 horse male-specific Y chromosome multicopy genes and SRY and evaluated their absolute copy numbers in 209 normal male horses of 22 breeds, 73 XY horses with disorders of sex development and/or infertility, 5 Przewalski's horses and 2 kulans. This established baseline copy number for these genes in horses. The TSPY gene showed the highest copy number and was the most copy number variable between individuals and breeds. SRY was a single-copy gene in most horses but had 2-3 copies in some indigenous breeds. Since SRY is flanked by 2 copies of RBMY, their copy number variations were interrelated and may lead to SRY-negative XY disorders of sex development. The Przewalski's horse and kulan had 1 copy of SRY and RBMY. TSPY and ETSTY2 showed significant copy number variations between cryptorchid and normal males (P < 0.05). No significant copy number variations were observed in subfertile/infertile males. Notably, copy number of TSPY and ETSTY5 differed between successive male generations and between cloned horses, indicating germline and somatic mechanisms for copy number variations. We observed no correlation between male-specific Y chromosome gene copy number variations and male-specific Y chromosome haplotypes. We conclude that the ampliconic male-specific Y chromosome reference assembly has deficiencies and further studies with an improved male-specific Y chromosome assembly are needed to determine selective constraints over horse male-specific Y chromosome gene copy number and their relation to stallion reproduction and male biology.

Accessing the Variability of Multicopy Genes in Complex Genomes using Unassembled Next-Generation Sequencing Reads: The Case of Trypanosoma cruzi Multigene Families.

Repetitive elements cause assembly fragmentation in complex eukaryotic genomes, limiting the study of their variability. The genome of Trypanosoma cruzi, the parasite that causes Chagas disease, has a high repetitive content, including multigene families. Although many T. cruzi multigene families encode surface proteins that play pivotal roles in host-parasite interactions, their variability is currently underestimated, as their high repetitive content results in collapsed gene variants. To estimate sequence variability and copy number variation of multigene families, we developed a read-based approach that is independent of gene-specific read mapping and de novo assembly. This methodology was used to estimate the copy number and variability of MASP, TcMUC, and Trans-Sialidase (TS), the three largest T. cruzi multigene families, in 36 strains, including members of all six parasite discrete typing units (DTUs). We found that these three families present a specific pattern of variability and copy number among the distinct parasite DTUs. Inter-DTU hybrid strains presented a higher variability of these families, suggesting that maintaining a larger content of their members could be advantageous. In addition, in a chronic murine model and chronic Chagasic human patients, the immune response was focused on TS antigens, suggesting that targeting TS conserved sequences could be a potential avenue to improve diagnosis and vaccine design against Chagas disease. Finally, the proposed approach can be applied to study multicopy genes in any organism, opening new avenues to access sequence variability in complex genomes. IMPORTANCE Sequences that have several copies in a genome, such as multicopy-gene families, mobile elements, and microsatellites, are among the most challenging genomic segments to study. They are frequently underestimated in genome assemblies, hampering the correct assessment of these important players in genome evolution and adaptation. Here, we developed a new methodology to estimate variability and copy numbers of repetitive genomic regions and employed it to characterize the T. cruzi multigene families MASP, TcMUC, and transsialidase (TS), which are important virulence factors in this parasite. We showed that multigene families vary in sequence and content among the parasite's lineages, whereas hybrid strains have a higher sequence variability that could be advantageous to the parasite's survivability. By identifying conserved sequences within multigene families, we showed that the mammalian host immune response toward these multigene families is usually focused on the TS multigene family. These TS conserved and immunogenic peptides can be explored in future works as diagnostic targets or vaccine candidates for Chagas disease. Finally, this methodology can be easily applied to any organism of interest, which will aid in our understanding of complex genomic regions.

Remodelling metabolism for high-level resveratrol production in Yarrowia lipolytica.

Resveratrol is a polyphenol with numerous applications in food, pharma, and cosmetics. Lack of precursors and low titer are the main problems hindering industrial scale resveratrol production. Based on previous prescreening, expressing the combination of FjTAL, Pc4CL1 and VvSTS achieved the best resveratrol titer. This was further improved to 235.1 mg/L through engineering the shikimic acid pathway, applying a modular enzyme assembly of Pc4CL1 and VvSTS, enhancing p-coumaric acid supply and diverting glycolytic flux toward erythrose-4-phosphate. The titer was increased to 819.1 mg/L following two rounds of multicopy integration of resveratrol biosynthesis and malonyl-CoA supply, respectively. The titer reached 22.5 g/L with a yield on glucose of 65.5 mg/g using an optimum fed-batch strategy in a 5 L bioreactor with morphology control. This research is the highest report on the de novo production of resveratrol in Yarrowia lipolytica and the findings lay a solid foundation for other producing polyphenols.

Cistanche Species Mitogenomes Suggest Diversity and Complexity in Lamiales-Order Mitogenomes.

The extreme diversity and complexity of angiosperms is well known. Despite the fact that parasitic plants are angiosperms, little is known about parasitic plant mitogenomic diversity, complexity, and evolution. In this study, we obtained and characterized the mitogenomes of three Cistanche species (holoparasitic plants) from China to compare the repeats, segment duplication and multi-copy protein-coding genes (PCGs), to clarify the phylogenetic and evolution relationship within the Lamiales order, and to identify the mitochondrial plastid insertions (MTPT) in Cistanche mitogenomes. The results showed that the mitogenome sizes of the three Cistanche species ranged from 1,708,661 to 3,978,341 bp. The Cistanche species genome encodes 75-126 genes, including 37-65 PCGs, 31-58 tRNA genes and 3-5 rRNA genes. Compared with other Lamiales and parasitic species, the Cistanche species showed extremely high rates of multi-copy PCGs, ranging from 0.13 to 0.58 percent of the total number of PCGs. In addition, 37-133 Simple Sequence Repeat (SSRs) were found in these three mitogenomes, the majority of which were the mononucleotides Adenine/Thymine. The interspersed repeats contained forward and palindromic repeats. Furthermore, the segment-duplication sequence size ranged from 199,584 to 2,142,551 bp, accounting for 24.9%, 11.7% and 53.9% of the Cistanche deserticola, Cistanche salsa and Cistanche tubulosa mitogenome, respectively. Furthermore, the Ka/Ks analysis suggested that the atp4, ccmB, ccmFc and matR were probably positively selected during Lamiales evolution. The Cistanche plastome suggested the presence of MTPT. Moreover, 6-12 tRNA, 9-15 PCGs fragments and 3 rRNA gene fragments in the Cistanche mitogenomes were found in the MTPT regions. This work reports the Cistanche species mitogenome for the first time, which will be invaluable for study the mitogenome evolution of Orobanchaceae family.

High-level secretory production of leghemoglobin in Pichia pastoris through enhanced globin expression and heme biosynthesis.

Soy leghemoglobin is a key food additive that imparts meaty flavor and color to meat analogs. Here, a Pichia pastoris strain capable of high-yield secretory production of functional leghemoglobin was developed through gene dosage optimization and heme pathway consolidation. First, multi-copy integration of LegH expression cassette was achieved via both post-transformational vector amplification and CRISPR/Cas9 mediated genome editing methods. A combination of inducible expression and constitutive expression resulted in the highest production of leghemoglobin. Then, heme biosynthetic pathway was engineered to address challenges in heme depletion and leghemoglobin secretion. Finally, the disruption of ku70 was complemented in engineered P. pastoris strain to enable high-density fermentation in a 10-L bioreactor. These engineering strategies increased the secretion of leghemoglobin by more than 83-fold, whose maximal leghemoglobin titer and heme binding ratio reached as high as 3.5 g/L and 93 %, respectively. This represents the highest secretory production of heme-containing proteins ever reported.