Complete mitochondrial genome of Episymploce splendens (Blattodea: Ectobiidae): A large intergenic spacer and lacking of two tRNA genes

The complete mitochondrial genome of Episymploce splendens, 15,802 bp in length, was determined and annotated in this study. The mito-genome included 13 PCGs, 20 tRNAs and 2 rRNAs. Unlike most typical mito-genomes with conservative gene arrangement and exceptional economic organization, E. splendens mito-genome has two tRNAs (tRNA-Gln and tRNA-Met) absence and a long intergenic spacer sequence (93 bp) between tRNA-Val and srRNA, showing the diversified features of insect mito-genomes. This is the first report of the tRNAs deletion in blattarian mito-genomes and we supported the duplication/random loss model as the origin mechanism of the long intergenic spacer. Two Numts, Numt-1 (557 bp) and Numt-2 (975 bp) transferred to the nucleus at about 14.15 Ma to 22.34 Ma, and 19.19 Ma to 24.06 Ma respectively, were found in E. splendens. They can be used as molecular fossils in insect phylogenetic relationship inference. Our study provided useful data for further studies on the evolution of insect mito-genome.

Complete mitochondrial genome of Ovalona pulchella (Branchiopoda, Cladocera) as the first representative in the family Chydoridae: Gene rearrangements and phylogenetic analysis of Cladocera

Chydoridae are phytophilic-benthic microcrustaceans that make up a significant proportion of species diversity and play an important role in the littoral zone of freshwater ecosystems worldwide. Here, we provide the complete mitochondrial genome of Ovalona pulchella (King, 1853), determined by next-generation sequencing. The entire mitochondrial genome is 15,362 bp in length; this is the first sequenced mitochondrial genome in the family Chydoridae. The base composition and codon usage were typical of Cladocera species. The mitochondrial gene arrangement (37 genes) was not consistent with that of other Branchiopoda. Both maximum likelihood and Bayesian analyses supported each suborder and family of Branchiopoda as monophyletic groups. The relationships among the families were as follows: [(Leptestheriidae + Limnadiidae) + (Sididae + (Bosminidae + (Chydoridae + Daphniidae)))] + Triopsidae. The newly sequenced O. pulchella was most closely related to the family Daphniidae. The complete mitochondrial genome of O. pulchella also provides valuable molecular information for further analysis of the phylogeny of the Chydoridae and the taxonomic status of the Branchiopoda.

Nonsense-mediated mRNA decay uses complementary mechanisms to suppress mRNA and protein accumulation

Nonsense-mediated mRNA decay (NMD) is an essential, highly conserved quality control pathway that detects and degrades mRNAs containing premature termination codons. Although the essentiality of NMD is frequently ascribed to its prevention of truncated protein accumulation, the extent to which NMD actually suppresses proteins encoded by NMD-sensitive transcripts is less well-understood than NMD-mediated suppression of mRNA. Here, we describe a reporter system that permits accurate quantification of both mRNA and protein levels via stable integration of paired reporters encoding NMD-sensitive and NMD-insensitive transcripts into the AAVS1 safe harbor loci in human cells. We use this system to demonstrate that NMD suppresses proteins encoded by NMD-sensitive transcripts by up to eightfold more than the mRNA itself. Our data indicate that NMD limits the accumulation of proteins encoded by NMD substrates by mechanisms beyond mRNA degradation, such that even when NMD-sensitive mRNAs escape destruction, their encoded proteins are still effectively suppressed.

IMGT®Homo sapiens IG and TR Loci, Gene Order, CNV and Haplotypes: New Concepts as a Paradigm for Jawed Vertebrates Genome Assemblies

IMGT^®, the international ImMunoGeneTics information system^®, created in 1989, by Marie-Paule Lefranc (Université de Montpellier and CNRS), marked the advent of immunoinformatics, a new science which emerged at the interface between immunogenetics and bioinformatics for the study of the adaptive immune responses. IMGT^® is based on a standardized nomenclature of the immunoglobulin (IG) and T cell receptor (TR) genes and alleles from fish to humans and on the IMGT unique numbering for the variable (V) and constant (C) domains of the immunoglobulin superfamily (IgSF) of vertebrates and invertebrates, and for the groove (G) domain of the major histocompatibility (MH) and MH superfamily (MhSF) proteins. IMGT^® comprises 7 databases, 17 tools and more than 25,000 pages of web resources for sequences, genes and structures, based on the IMGT Scientific chart rules generated from the IMGT-ONTOLOGY axioms and concepts. IMGT^® reference directories are used for the analysis of the NGS high-throughput expressed IG and TR repertoires (natural, synthetic and/or bioengineered) and for bridging sequences, two-dimensional (2D) and three-dimensional (3D) structures. This manuscript focuses on the IMGT^®Homo sapiens IG and TR loci, gene order, copy number variation (CNV) and haplotypes new concepts, as a paradigm for jawed vertebrates genome assemblies.

Mitogenome sequence comparison in the endangered congeneric Pinna nobilis and Pinna rudis bivalves

BACKGROUND

The pen shells Pinna nobilis and Pinna rudis are large wedge-shaped bivalve molluscs. Both species are threatened by different anthropogenic pressures. In the last few years, P. nobilis populations have significantly reduced due to massive mortality events. The complete mitochondrial DNA sequences of these congeneric species have been determined and compared for the first time.

RESULTS

The mitogenome sequences of P. nobilis and P. rudis were 18,919 bp and 18,264 bp in length, respectively. Each mitogenome is composed of 12 protein-coding genes, 2 ribosomal RNA, 22 transfer RNA (tRNAs) genes and non-coding regions. A putative Adenosine Triphosphate synthase subunit 8 gene could only be proposed for P. nobilis. Both newly sequenced mitogenomes present a conserved gene order between them, comparable to the closely related Atrina pectinata, but global arrangement greatly differs from other available bivalve mitochondrial sequences. Multiple copies of tRNA-Cys were identified, located in different positions probably due to mechanisms of mitochondrial genome rearrangements, and detected 2 and 3 times in P. rudis and in P. nobilis, respectively.

CONCLUSION

A close relationship was shown between Pinna species and Atrina pectinata and a consistent clustering showing a monophyletic origin of Pinnidae family sequences was evidenced. The mitochondrial genomes will provide a valuable genetic resource for further studies on population genetics and species identification.

Powassan Viruses Spread Cell to Cell during Direct Isolation from Ixodes Ticks and Persistently Infect Human Brain Endothelial Cells and Pericytes

Powassan viruses (POWVs) are neurovirulent tick-borne flaviviruses emerging in the northeastern United States, with a 2% prevalence in Long Island (LI) deer ticks (Ixodes scapularis). POWVs are transmitted within as little as 15 min of a tick bite and enter the central nervous system (CNS) to cause encephalitis (10% of cases are fatal) and long-term neuronal damage. POWV-LI9 and POWV-LI41 present in LI Ixodes ticks were isolated by directly inoculating VeroE6 cells with tick homogenates and detecting POWV-infected cells by immunoperoxidase staining. Inoculated POWV-LI9 and LI41 were exclusively present in infected cell foci, indicative of cell to cell spread, despite growth in liquid culture without an overlay. Cloning and sequencing establish POWV-LI9 as a phylogenetically distinct lineage II POWV strain circulating in LI deer ticks. Primary human brain microvascular endothelial cells (hBMECs) and pericytes form a neurovascular complex that restricts entry into the CNS. We found that POWV-LI9 and -LI41 and lineage I POWV-LB productively infect hBMECs and pericytes and that POWVs were basolaterally transmitted from hBMECs to lower-chamber pericytes without permeabilizing polarized hBMECs. Synchronous POWV-LI9 infection of hBMECs and pericytes induced proinflammatory chemokines, interferon-β (IFN-β) and proteins of the IFN-stimulated gene family (ISGs), with delayed IFN-β secretion by infected pericytes. IFN inhibited POWV infection, but despite IFN secretion, a subset of POWV-infected hBMECs and pericytes remained persistently infected. These findings suggest a potential mechanism for POWVs (LI9/LI41 and LB) to infect hBMECs, spread basolaterally to pericytes, and enter the CNS. hBMEC and pericyte responses to POWV infection suggest a role for immunopathology in POWV neurovirulence and potential therapeutic targets for preventing POWV spread to neuronal compartments. IMPORTANCE We isolated POWVs from LI deer ticks (I. scapularis) directly in VeroE6 cells, and sequencing revealed POWV-LI9 as a distinct lineage II POWV strain. Remarkably, inoculation of VeroE6 cells with POWV-containing tick homogenates resulted in infected cell foci in liquid culture, consistent with cell-to-cell spread. POWV-LI9 and -LI41 and lineage I POWV-LB strains infected hBMECs and pericytes that comprise neurovascular complexes. POWVs were nonlytically transmitted basolaterally from infected hBMECs to lower-chamber pericytes, suggesting a mechanism for POWV transmission across the blood-brain barrier (BBB). POWV-LI9 elicited inflammatory responses from infected hBMEC and pericytes that may contribute to immune cell recruitment and neuropathogenesis. This study reveals a potential mechanism for POWVs to enter the CNS by infecting hBMECs and spreading basolaterally to abluminal pericytes. Our findings reveal that POWV-LI9 persists in cells that form a neurovascular complex spanning the BBB and suggest potential therapeutic targets for preventing POWV spread to neuronal compartments.

Sequencing of the Complete Mitochondrial Genome of Pingus sinensis (Spirurina: Quimperiidae): Gene Arrangements and Phylogenetic Implications

Despite several decades of intensive research on spirurine nematodes, molecular data on some of the main lineages are still absent, which makes taxonomic classification insufficiently resolved. In the present study, we sequenced the first complete mitogenome for the family Quimperiidae, belonging to P. sinensis (Spirurina: Quimperiidae), a parasite living in the intestines of snakehead (Ophiocephalus argus). The circular mitogenome is 13,874 bp long, and it contains the standard nematode gene set: 22 transfer RNAs, 2 ribosomal RNAs and 12 protein-coding genes. There are also two long non-coding regions (NCR), in addition to only 8 other intergenic regions, ranging in size from 1 to 58 bp. To investigate its phylogenetic position and study the relationships among other available Spirurina, we performed the phylogenetic analysis using Bayesian inference and maximum likelihood approaches by concatenating the nucleotide sequences of all 36 genes on a dataset containing all available mitogenomes of the suborder Spirurina from NCBI and compared with gene order phylogenies using the MLGO program. Both supported the closer relationship of Ascaridoidea to Seuratoidea than to Spiruroidea. Pingus formed a sister-group with the Cucullanus genus. The results provide a new insights into the relationships within Spirurina.

The complete mitogenome of Lysmata vittata (Crustacea: Decapoda: Hippolytidae) with implication of phylogenomics and population genetics

In this study, the complete mitogenome of Lysmata vittata (Crustacea: Decapoda: Hippolytidae) has been determined. The genome sequence was 22003 base pairs (bp) and it included thirteen protein-coding genes (PCGs), twenty-two transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and three putative control regions (CRs). The nucleotide composition of AT was 71.50%, with a slightly negative AT skewness (-0.04). Usually the standard start codon of the PCGs was ATN, while cox1, nad4L and cox3 began with TTG, TTG and GTG. The canonical termination codon was TAA, while nad5 and nad4 ended with incomplete stop codon T, and cox1 ended with TAG. The mitochondrial gene arrangement of eight species of the Hippolytidae were compared with the order of genes of Decapoda ancestors, finding that the gene arrangement order of the Lebbeus groenlandicus had not changed, but the gene arrangement order of other species changed to varying degrees. The positions of the two tRNAs genes (trnA and trnR) of the L. vittata had translocations, which also showed that the Hippolytidae species were relatively unconserved in evolution. Phylogenetic analysis of 50 shrimp showed that L. vittata formed a monophyletic clade with Lysmata/Exhippolysmata species. This study should be helpful to better understand the evolutionary status, and population genetic diversity of L. vittata and related species.

Co-opting regulation bypass repair as a gene-correction strategy for monogenic diseases

With the development of CRISPR-Cas9-mediated gene-editing technologies, correction of disease-causing mutations has become possible. However, current gene-correction strategies preclude mutation repair in post-mitotic cells of human tissues, and a unique repair strategy must be designed and tested for each and every mutation that may occur in a gene. We have developed a novel gene-correction strategy, co-opting regulation bypass repair (CRBR), which can repair a spectrum of mutations in mitotic or post-mitotic cells and tissues. CRBR utilizes the non-homologous end joining (NHEJ) pathway to insert a coding sequence (CDS) and transcription/translation terminators targeted upstream of any CDS mutation and downstream of the transcriptional promoter. CRBR results in simultaneous co-option of the endogenous regulatory region and bypass of the genetic defect. We validated the CRBR strategy for human gene therapy by rescuing a mouse model of Wolcott-Rallison syndrome (WRS) with permanent neonatal diabetes caused by either a large deletion or a nonsense mutation in the PERK (EIF2AK3) gene. Additionally, we integrated a CRBR GFP-terminator cassette downstream of the human insulin promoter in cadaver pancreatic islets of Langerhans, which resulted in insulin promoter regulated expression of GFP, demonstrating the potential utility of CRBR in human tissue gene repair.

Full-length dystrophin restoration via targeted exon integration by AAV-CRISPR in a humanized mouse model of Duchenne muscular dystrophy

Targeted gene-editing strategies have emerged as promising therapeutic approaches for the permanent treatment of inherited genetic diseases. However, precise gene correction and insertion approaches using homology-directed repair are still limited by low efficiencies. Consequently, many gene-editing strategies have focused on removal or disruption, rather than repair, of genomic DNA. In contrast, homology-independent targeted integration (HITI) has been reported to effectively insert DNA sequences at targeted genomic loci. This approach could be particularly useful for restoring full-length sequences of genes affected by a spectrum of mutations that are also too large to deliver by conventional adeno-associated virus (AAV) vectors. Here, we utilize an AAV-based, HITI-mediated approach for correction of full-length dystrophin expression in a humanized mouse model of Duchenne muscular dystrophy (DMD). We co-deliver CRISPR-Cas9 and a donor DNA sequence to insert the missing human exon 52 into its corresponding position within the DMD gene and achieve full-length dystrophin correction in skeletal and cardiac muscle. Additionally, as a proof-of-concept strategy to correct genetic mutations characterized by diverse patient mutations, we deliver a superexon donor encoding the last 28 exons of the DMD gene as a therapeutic strategy to restore full-length dystrophin in >20% of the DMD patient population. This work highlights the potential of HITI-mediated gene correction for diverse DMD mutations and advances genome editing toward realizing the promise of full-length gene restoration to treat genetic disease.

The dynamic history of gymnosperm plastomes: Insights from structural characterization, comparative analysis, phylogenomics, and time divergence

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution.

Phylogenetic relationships in the Sorghum genus based on sequencing of the chloroplast and nuclear genes

Sorghum [Sorghum bicolor (L.) Moench] is an important food crop with a diverse gene pool residing in its wild relatives. A total of 15 sorghum accessions from the unexploited wild gene pool of the Sorghum genus, representing the five subgenera, were sequenced, and the complete chloroplast genomes and 99 common single-copy concatenated nuclear genes were assembled. Annotation of the chloroplast genomes identified a total of 81 protein-coding genes, 38 tRNA, and four rRNA genes. The gene content and gene order among the species was identical. A total of 153 nonsynonymous amino acid changes in 40 genes were identified across the species. Phylogenetic analysis of both the whole chloroplast genome and nuclear genes revealed a similar topology with two distinct clades within the genus. The species within the subgenera Eusorghum, Chaetosorghum, and Heterosorghum clustered in one clade, whereas the species within the subgenera Parasorghum and Stiposorghum clustered in a second clade. However, the subgenera Parasorghum and Stiposorghum were not monophyletic, suggesting the need for further research to resolve the relationships within this group. The close relationship between the two monotypic subgenera Chaetosorghum and Heterosorghum suggests that species within these subgenera could be considered as one group. This analysis provides an improved understanding of the genetic relationships within the Sorghum genus and defines diversity in wild sorghum species that may be useful for crop improvement.

A chromosome-scale genome assembly and karyotype of the ctenophore Hormiphora californensis

Here, we present a karyotype, a chromosome-scale genome assembly, and a genome annotation from the ctenophore Hormiphora californensis (Ctenophora: Cydippida: Pleurobrachiidae). The assembly spans 110 Mb in 44 scaffolds and 99.47% of the bases are contained in 13 scaffolds. Chromosome micrographs and Hi-C heatmaps support a karyotype of 13 diploid chromosomes. Hi-C data reveal three large heterozygous inversions on chromosome 1, and one heterozygous inversion shares the same gene order found in the genome of the ctenophore Pleurobrachia bachei. We find evidence that H. californensis and P. bachei share thirteen homologous chromosomes, and the same karyotype of 1n = 13. The manually curated PacBio Iso-Seq-based genome annotation reveals complex gene structures, including nested genes and trans-spliced leader sequences. This chromosome-scale assembly is a useful resource for ctenophore biology and will aid future studies of metazoan evolution and phylogenetics.

Comparative mitogenomics of Drosophilidae and the evolution of the Zygothrica genus group (Diptera, Drosophilidae)

The Zygothrica genus group of Drosophilidae encompasses more than 437 species and five genera. Although knowledge regarding its diversity has increased, uncertainties about its monophyly and position within Drosophilidae remain. Genomic approaches have been widely used to address different phylogenetic questions and analyses involving the mitogenome have revealed a cost-efficient tool to these studies. Thus, this work aims to characterize mitogenomes of three species of the Zygothrica genus group (from the Hirtodrosophila, Paraliodrosophila and Zygothrica genera), while comparing them with orthologous sequences from other 23 Drosophilidae species and addressing their phylogenetic position. General content concerning gene order and overlap, nucleotide composition, start and stop codon, codon usage and tRNA structures were compared, and phylogenetic trees were constructed under different datasets. The complete mitogenomes characterized for H. subflavohalterata affinis H002 and P. antennta present the PanCrustacea gene order with 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, 13 protein coding genes and an A+T rich region with two T-stretched elements. Some peculiarities such as the almost complete overlap of genes tRNAH/ND4, tRNAF/ND5 and tRNAS2/ND1 are reported for different Drosophilidae species. Non-canonical secondary structures were encountered for tRNAS1 and tRNAY, revealing patterns that apply at different phylogenetic scales. According to the best depiction of the mitogenomes evolutionary history, the three Neotropical species of the Zygothrica genus group encompass a monophyletic lineage sister to Zaprionus, composing with this genus a clade that is sister to the Drosophila subgenus.

Chimeric antigen receptor-modified human regulatory T cells that constitutively express IL-10 maintain their phenotype and are potently suppressive

Clinical trials of Treg therapy in transplantation are currently entering phases IIa and IIb, with the majority of these employing polyclonal Treg populations that harbor a broad specificity. Enhancing Treg specificity is possible with the use of chimeric antigen receptors (CARs), which can be customized to respond to a specific human leukocyte antigen (HLA). In this study, we build on our previous work in the development of HLA-A2 CAR-Tregs by further equipping cells with the constitutive expression of interleukin 10 (IL-10) and an imaging reporter as additional payloads. Cells were engineered to express combinations of these domains and assessed for phenotype and function. Cells expressing the full construct maintained a stable phenotype after transduction, were specifically activated by HLA-A2, and suppressed alloresponses potently. The addition of IL-10 provided an additional advantage to suppressive capacity. This study therefore provides an important proof-of-principle for this cell engineering approach for next-generation Treg therapy in transplantation.

Efficient biallelic knock-in in mouse embryonic stem cells by in vivo-linearization of donor and transient inhibition of DNA polymerase θ/DNA-PK

CRISPR/Cas9-mediated homology-directed repair (HDR) is used for error-free targeted knock-in of foreign donor DNA. However, the low efficiency of HDR-mediated knock-in hinders establishment of knock-in clones. Double-strand breaks (DSBs) induced by CRISPR/Cas9 are preferentially repaired by non-homologous end joining (NHEJ) or microhomology-mediated end joining (MMEJ) before HDR can occur, thereby preventing HDR-mediated knock-in. NHEJ/MMEJ also cause random integrations, which give rise to false-positive knock-in events, or silently disrupt the genome. In this study, we optimized an HDR-mediated knock-in method for mouse embryonic stem cells (mESCs). We succeeded in improving efficiency of HDR-mediated knock-in of a plasmid donor while almost completely suppressing NHEJ/MMEJ-based integration by combining in vivo-linearization of the donor plasmid, transient knockdown of DNA polymerase θ, and chemical inhibition of DNA-dependent protein kinase (DNA-PK) by M3814. This method also dramatically improved the efficiency of biallelic knock-in; at the Rosa26a locus, 95% of HDR-mediated knock-in clones were biallelic. We designate this method BiPoD (Biallelic knock-in assisted by Pol θ and DNA-PK inhibition). BiPoD achieved simultaneous efficient biallelic knock-in into two loci. BiPoD, therefore, enables rapid and easy establishment of biallelic knock-in mESC lines.

Characterization of the complete mitochondrial genomes of two Ixodes ticks, I. nipponensis and Ixodes (Pholeoixodes) sp

In this study, the authors sequenced and characterized the complete mitochondrial (mt) genomes of two hard ticks of the genus Ixodes, I. nipponensis and Ixodes (Pholeoixodes) sp., which were 14 505 and 14 543 bp in length, respectively. Their mt genomes encoded 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes and two ribosomal RNA genes, and have only one non-coding region. The gene order in their mt genomes was the same as that of other Ixodes spp. mt genomes. The average sequence identity, combined nucleotide diversity, non-synonymous/synonymous substitutions ratio analyses consistently demonstrated that cox1, rrnS, cox2, cox3 and cytb were the most conserved and atp8, nad6 and nad2 were the most variable genes across Ixodes mitogenomes. Phylogeny of the present Ixodes spp., and other selected hard tick species, based on concatenated amino acid sequences of PCGs, confirmed their position within the genus Ixodes and sub-family Ixodinae. The novel mt markers described herein will be useful for further studies of the population genetics, molecular epidemiology and systematics of hard ticks.

Mitogenomics of Cladocera (Branchiopoda): Marked gene order rearrangements and independent predation roots

Cladocera (Crustacea: Branchiopoda) is a key group of invertebrates. Despite a long history of phylogenetic research, relationships within this group remain disputed. We here provide new insights based on 15 new mitochondrial genomes obtained from high-throughput sequencing (HTS) and 40 mitogenomes extracted from published HTS datasets. Together with 25 mitogenomes from GenBank, we generated a matrix of 80 mitogenomes, 44 of them belonging to Cladocera. We also obtained a matrix with 168 nuclear orthologous genes to further assess the phylogenetic result from mitogenomes based on published data and one new HTS data ofLeptodora. Maximum likelihood and Bayesian phylogenetic analyses recovered all Branchiopoda orders as monophyletic and supported a sister-group relationship between Anomopoda and Onychopoda, making the taxon Gymnomera paraphyletic and supporting an independent origin of predatory Haplopoda and Onychopoda. The nuclear phylogeny and topological tests also support Gymnomera as paraphyletic, and the nuclear phylogeny strongly supports a sister-group relationship between Ctenopoda and Haplopoda. We provide a fossil-calibrated time tree, congruent with a Carboniferous origin for Cladocera and a subsequent diversification of the crown group of Anomopoda, Onychopoda, and Ctenopoda, at least in the Triassic. Despite their long evolutionary history, non-Cladoceran Branchiopoda exhibited high mitogenome structural stability. On the other hand, 21 out of 24 gene rearrangements occurred within the relatively younger Cladocera. We found the differential base compositional skewness patterns between Daphnia s.s. and Ctenodaphnia, which might be related to the divergence between these taxa. We also provide evidence to support the recent finding that Spinicaudata possesses mitogenomes with inversed compositional skewness without gene rearrangement. Such a pattern has only been reported in Spinicaudata.

The complete mitochondrial genome of Cycas debaoensis revealed unexpected static evolution in gymnosperm species

Mitochondrial genomes of vascular plants are well known for their liability in architecture evolution. However, the evolutionary features of mitogenomes at intra-generic level are seldom studied in vascular plants, especially among gymnosperms. Here we present the complete mitogenome of Cycas debaoensis, an endemic cycad species to the Guangxi region in southern China. In addition to assemblage of draft mitochondrial genome, we test the conservation of gene content and mitogenomic stability by comparing it to the previously published mitogenome of Cycas taitungensis. Furthermore, we explored the factors such as structural rearrangements and nuclear surveillance of double-strand break repair (DSBR) proteins in Cycas in comparison to other vascular plant groups. The C. debaoensis mitogenome is 413,715 bp in size and encodes 69 unique genes, including 40 protein coding genes, 26 tRNAs, and 3 rRNA genes, similar to that of C. taitungensis. Cycas mitogenomes maintained the ancestral intron content of seed plants (26 introns), which is reduced in other lineages of gymnosperms, such as Ginkgo biloba, Taxus cuspidata and Welwitschia mirabilis due to selective pressure or retroprocessing events. C. debaoensis mitogenome holds 1,569 repeated sequences (> 50 bp), which partially account for fairly large intron size (1200 bp in average) of Cycas mitogenome. The comparison of RNA-editing sites revealed 267 shared non-silent editing site among predicted vs. empirically observed editing events. Another 33 silent editing sites from empirical data increase the total number of editing sites in Cycas debaoensis mitochondrial protein coding genes to 300. Our study revealed unexpected conserved evolution between the two Cycas species. Furthermore, we found strict collinearity of the gene order along with the identical set of genomic content in Cycas mt genomes. The stability of Cycas mt genomes is surprising despite the existence of large number of repeats. This structural stability may be related to the relative expansion of three DSBR protein families (i.e., RecA, OSB, and RecG) in Cycas nuclear genome, which inhibit the homologous recombinations, by monitoring the accuracy of mitochondrial chromosome repair.

Gene Position Index Mutation Detection Algorithm Based on Feedback Fast Learning Neural Network

In the detection of genome variation, the research on the internal correlation of reference genome is deepening; the detection of variation in genome sequence has become the focus of research, and it has also become an effective path to find new genes and new functional proteins. The targeted sequencing sequence is used to sequence the exon region of a specific gene in cancer gene detection, and the sequencing depth is relatively large. Traditional alignment algorithms will lose some sequences, which will lead to inaccurate mutation detection. This paper proposes a mutation detection algorithm based on feedback fast learning neural network position index. By establishing a position index relationship for ACGT in the DNA sequence, the subsequence is decomposed into the position relationship of different subsequences corresponding to the main sequence. The positional relationship of the subsequence in the main sequence is determined by the positional relationship. Analyzing SNP and InDel mutations, even structural mutations, through the position correlation of sequences has the advantages of high precision and easy implementation by personal computers. The feedback fast learning neural network is used to verify whether there is a linear relationship between two or more positions. Experimental results show that the mutation points detected by position index are more than those detected by Bcftools, Freebye, Vanscan2, and Gatk.

Characterization and Phylogenetic Analysis of the Mitochondrial Genome Sequence of Nisia fuliginosa (Hemiptera: Fulgoroidea: Meenoplidae)

We explored characterization of the mitochondrial genome (mitogenome or mtGenome) and phylogenetic analysis between 32 Fulgoroid species by sequencing and analyzing the mitogenome of Nisia fuliginosa Yang and Hu, 1985 (Hemiptera: Fulgoroidea: Meenoplidae), thereby making it the first determined mitogenome from the family Meenoplidae. The mitogenome was found to be 15,754 bp in length and contained 13 protein-coding genes (PCGs), 22 tRNA genes, two ribosomal RNA genes (rRNAs), and a control region. All PCGs started with typical ATN codons, except for nad1, which used GTG as the start codon. Canonical TAA termination codons were found in 10 PCGs and the remaining three genes (cox2, nad6, and nad1) had incomplete stop codons T. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnC, trnV, and trnS1. Additionally, we compared the AT and GC skews of 13 PCGs of 32 Fulgoroidea mitogenomes, on the L-strand, the AT and GC skews were negative and positive, respectively. However, on the H-strand, the AT skew could be positive or negative and the GC skew was always negative. Phylogenetic results showed that the eight families of Fulgoroidea were divided into two large groups. Delphacidae formed a monophyletic group sister to a clade comprising Meenoplidae and other six families (Fulgoridae, Ricaniidae, Flatidae, Issidae, Caliscelidae, and Achilidae). Meenoplidae was located near the clade of Delphacidae, and Fulgoridae was located near the clade of Meenoplidae. Furthermore, Caliscelidae, Issidae, Ricaniidae, and Flatidae are closely related and they collectively formed a sister group to Achilidae.

Chromosomal-scale de novo genome assemblies of Cynomolgus Macaque and Common Marmoset

Cynomolgus macaque (Macaca fascicularis) and common marmoset (Callithrix jacchus) have been widely used in human biomedical research. Long-standing primate genome assemblies used the human genome as a reference for ordering and orienting the assembled fragments into chromosomes. Here we performed de novo genome assembly of these two species without any human genome-based bias observed in the genome assemblies released earlier. We assembled PacBio long reads, and the resultant contigs were scaffolded with Hi-C data, which were further refined based on Hi-C contact maps and alternate de novo assemblies. The assemblies achieved scaffold N50 lengths of 149 Mb and 137 Mb for cynomolgus macaque and common marmoset, respectively. The high fidelity of our assembly is also ascertained by BAC-end concordance in common marmoset. Our assembly of cynomolgus macaque outperformed all the available assemblies of this species in terms of contiguity. The chromosome-scale genome assemblies produced in this study are valuable resources for non-human primate models and provide an important baseline in human biomedical research.

A transgenic female killing system for the genetic control of Drosophila suzukii

The spotted wing Drosophila (Drosophila suzukii) is an invasive pest of soft-skinned fruit crops. It is rapidly transmitted in Europe and North America, causing widespread agricultural losses. Genetic control strategies such as the sterile insect technique (SIT) have been proposed as environment-friendly and species-restricted approaches for this pest. However, females are inefficient agents in SIT programs. Here we report a conditional female-killing (FK) strategy based on the tetracycline-off system. We assembled sixteen genetic constructs for testing in vitro and in vivo. Twenty-four independent transgenic strains of D. suzukii were generated and tested for female-specific lethality. The strongest FK effect in the absence of tetracycline was achieved by the construct containing D. suzukii nullo promoter for early gene expression, D. suzukii pro-apoptotic gene hidAla4 for lethality, and the transformer gene intron from the Mediterranean fruit fly Ceratitis capitata for female-specific splicing. One strain carrying this construct eliminated 100% of the female offspring during embryogenesis and produced only males. However, homozygous females from these FK strains were not viable on a tetracycline-supplemented diet, possibly due to the basal expression of hidAla4. Potential improvements to the gene constructs and the use of such FK strains in an SIT program are discussed.

Phylogeny and Comparative Analysis for the Plastid Genomes of Five Tulipa (Liliaceae)

Species of Tulipa (Liliaceae) are of great horticultural importance and are distributed across Europe, North Africa, and Asia. The Tien Shan Mountain is one of the primary diversity centres of Tulipa, but the molecular studies of Tulipa species from this location are lacking. In our study, we assembled four Tulipa plastid genomes from the Tien Shan Mountains, T. altaica, T. iliensis, T. patens, and T. thianschanica, combined with the plastid genome of T. sylvestris to compare against other Liliaceae plastid genomes. We focussed on the species diversity and evolution of their plastid genomes. The five Tulipa plastid genomes proved highly similar in overall size (151,691-152,088 bp), structure, gene order, and content. With comparative analysis, we chose 7 mononucleotide SSRs from the Tulipa species that could be used in further population studies. Phylogenetic analyses based on 24 plastid genomes robustly supported the monophyly of Tulipa and the sister relationship between Tulipa and Amana, Erythronium. T. iliensis, T. thianschanica, and T. altaica were clustered together, and T. patens was clustered with T. sylvestris, with our results clearly demonstrating the relationships between these five Tulipa species. Our results provide a more comprehensive understanding of the phylogenomics and comparative genomics of Tulipa.

A Novel Plasmid DNA-Based Foot and Mouth Disease Virus Minigenome for Intracytoplasmic mRNA Production

Picornaviruses are non-enveloped, single-stranded RNA viruses that cause highly contagious diseases, such as polio and hand, foot-and-mouth disease (HFMD) in human, and foot-and-mouth disease (FMD) in animals. Reverse genetics and minigenome of picornaviruses mainly depend on in vitro transcription and RNA transfection; however, this approach is inefficient due to the rapid degradation of RNA template. Although DNA-based reverse genetics systems driven by mammalian RNA polymerase I and/or II promoters display the advantage of rescuing the engineered FMDV, the enzymatic functions are restricted in the nuclear compartment. To overcome these limitations, we successfully established a novel DNA-based vector, namely pKLS3, an FMDV minigenome containing the minimum cis-acting elements of FMDV essential for intracytoplasmic transcription and translation of a foreign gene. A combination of pKLS3 minigenome and the helper plasmids yielded the efficient production of uncapped-green florescent protein (GFP) mRNA visualized in the transfected cells. We have demonstrated the application of the pKLS3 for cell-based antiviral drug screening. Not only is the DNA-based FMDV minigenome system useful for the FMDV research and development but it could be implemented for generating other picornavirus minigenomes. Additionally, the prospective applications of this viral minigenome system as a vector for DNA and mRNA vaccines are also discussed.