Gamma-rays induced genome wide stable mutations in cowpea deciphered through whole genome sequencing

PURPOSE

Gamma rays are the most widely exploited physical mutagen in plant mutation breeding. They are known to be involved in the development of more than 60% of global cowpea (Vigna unguiculata (L.) Walp.) mutant varieties. Nevertheless, the nature and type of genome-wide mutations induced by gamma rays have not been studied in cowpea and therefore, the present investigation was undertaken.

MATERIALS AND METHODS

Genomic DNAs from three stable gamma rays-induced mutants (large seed size, small seed size and disease resistant mutant) of cowpea cultivar 'CPD103' in M6 generation along with its progenitor were used for Illumina-based whole-genome resequencing.

RESULTS

Gamma rays induced a relatively higher frequency (88.9%) of single base substitutions (SBSs) with an average transition to transversion ratio (Ti/Tv) of 3.51 in M6 generation. A > G transitions, including its complementary T > C transitions, predominated the transition mutations, while all four types of transversion mutations were detected with frequencies over 6.5%. Indels (small insertions and deletions) constituted about 11% of the total induced variations, wherein small insertions (6.3%) were relatively more prominent than small deletions (4.8%). Among the indels, single-base indels and, in particular, those involving A/T bases showed a preponderance, albeit indels of up to three bases were detected in low proportions. Distributed across all 11 chromosomes, only a fraction of SBSs (19.45%) and indels (20.2%) potentially altered the encoded amino acids/peptides. The inherent mutation rate induced by gamma rays in cowpea was observed to be in the order of 1.4 × 10-7 per base pair in M6 generation.

CONCLUSION

Gamma-rays with a greater tendency to induce SBSs and, to a lesser extent, indels could be efficiently and effectively exploited in cowpea mutation breeding.

Phylogenomics and plastomics offer new evolutionary perspectives on Kalanchoideae (Crassulaceae)

BACKGROUND AND AIMS

Kalanchoideae is one of three subfamilies within Crassulaceae and contains four genera. Despite previous efforts, the phylogeny of Kalanchoideae remains inadequately resolved with persistent issues including low support, unstructured topologies and polytomies. This study aimed to address two central objectives: (1) resolving the pending phylogenetic questions within Kalanchoideae by using organelle-scale 'barcodes' (plastomes) and nuclear data; and (2) investigating interspecific diversity patterns among Kalanchoideae plastomes.

METHODS

To explore the plastome evolution in Kalanchoideae, we newly sequenced 38 plastomes representing all four constituent genera (Adromischus, Cotyledon, Kalanchoe and Tylecodon). We performed comparative analyses of plastomic features, including GC and gene contents, gene distributions at the IR (inverted repeat) boundaries, nucleotide divergence, plastomic tRNA (pttRNA) structures and codon aversions. Additionally, phylogenetic inferences were inferred using both the plastomic dataset (79 genes) and nuclear dataset (1054 genes).

KEY RESULTS

Significant heterogeneities were observed in plastome lengths among Kalanchoideae, strongly correlated with LSC (large single copy) lengths. Informative diversities existed in the gene content at SSC/IRa (small single copy/inverted repeat a), with unique patterns individually identified in Adromischus leucophyllus and one major Kalanchoe clade. The ycf1 gene was assessed as a shared hypervariable region among all four genera, containing nine lineage-specific indels. Three pttRNAs exhibited unique structures specific to Kalanchoideae and the genera Adromischus and Kalanchoe. Moreover, 24 coding sequences revealed a total of 41 lineage-specific unused codons across all four constituent genera. The phyloplastomic inferences clearly depicted internal branching patterns in Kalanchoideae. Most notably, by both plastid- and nuclear-based phylogenies, our research offers the first evidence that Kalanchoe section Eukalanchoe is not monophyletic.

CONCLUSIONS

This study conducted comprehensive analyses on 38 newly reported Kalanchoideae plastomes. Importantly, our results not only reconstructed well-resolved phylogenies within Kalanchoideae, but also identified highly informative unique markers at the subfamily, genus and species levels. These findings significantly enhance our understanding of the evolutionary history of Kalanchoideae.

New Genera and Species of Caulobacter and Brevundimonas Bacteriophages Provide Insights into Phage Genome Evolution

Previous studies have identified diverse bacteriophages that infect Caulobacter vibrioides strain CB15 ranging from small RNA phages to four genera of jumbo phages. In this study, we focus on 20 bacteriophages whose genomes range from 40 to 60 kb in length. Genome comparisons indicated that these diverse phages represent six Caulobacter phage genera and one additional genus that includes both Caulobacter and Brevundimonas phages. Within species, comparisons revealed that both single base changes and inserted or deleted genetic material cause the genomes of closely related phages to diverge. Among genera, the basic gene order and the orientation of key genes were retained with most of the observed variation occurring at ends of the genomes. We hypothesize that the nucleotide sequences of the ends of these phage genomes are less important than the need to maintain the size of the genome and the stability of the corresponding mRNAs.

Reversal and Transposition Distance on Unbalanced Genomes Using Intergenic Information

The most common way to calculate the rearrangement distance between two genomes is to use the size of a minimum length sequence of rearrangements that transforms one of the two given genomes into the other, where the genomes are represented as permutations using only their gene order, based on the assumption that genomes have the same gene content. With the advance of research in genome rearrangements, new works extended the classical models by either considering genomes with different gene content (unbalanced genomes) or including more genomic characteristics to the mathematical representation of the genomes, such as the distribution of intergenic regions sizes. In this study, we study the Reversal, Transposition, and Indel (Insertion and Deletion) Distance using intergenic information, which allows comparing unbalanced genomes, because indels are included in the rearrangement model (i.e., the set of possible rearrangements allowed when we compute the distance). For the particular case of transpositions and indels on unbalanced genomes, we present a 4-approximation algorithm, improving a previous 4.5 approximation. This algorithm is extended so as to deal with gene orientation and to maintain the 4-approximation factor for the Reversal, Transposition, and Indel Distance on unbalanced genomes. Furthermore, we evaluate the proposed algorithms using experiments on simulated data.

Experimental estimates of germline mutation rate in eukaryotes: a phylogenetic meta-analysis

Mutation is the ultimate source of all genetic variation, and over the last 10 years the ready availability of whole-genome sequencing has permitted direct estimation of mutation rate for many non-model species across the tree of life. In this meta-analysis, we make a comprehensive search of the literature for mutation rate estimates in eukaryotes, identifying 140 mutation accumulation (MA) and parent-offspring (PO) sequencing studies covering 134 species. Based on these data, we revisit differences in the single-nucleotide mutation (SNM) rate between different phylogenetic lineages and update the known relationships between mutation rate and generation time, genome size, and nucleotide diversity-while accounting for phylogenetic nonindependence. We do not find a significant difference between MA and PO in estimated mutation rates, but we confirm that mammal and plant lineages have higher mutation rates than arthropods and that unicellular eukaryotes have the lowest mutation rates. We find that mutation rates are higher in species with longer generation times and larger genome sizes, even when accounting for phylogenetic relationships. Moreover, although nucleotide diversity is positively correlated with mutation rate, the gradient of the relationship is significantly less than one (on a logarithmic scale), consistent with higher mutation rates in populations with smaller effective size. For the 29 species for which data are available, we find that indel mutation rates are positively correlated with nucleotide mutation rates and that short deletions are generally more common than short insertions. Nevertheless, despite recent progress, no estimates of either SNM or indel mutation rates are available for the majority of deeply branching eukaryotic lineages-or even for most animal phyla. Even among charismatic megafauna, experimental mutation rate estimates remain unknown for amphibia and scarce for reptiles and fish.

Recent Genome-Editing Approaches toward Post-Implanted Fetuses in Mice

Genome editing, as exemplified by the CRISPR/Cas9 system, has recently been employed to effectively generate genetically modified animals and cells for the purpose of gene function analysis and disease model creation. There are at least four ways to induce genome editing in individuals: the first is to perform genome editing at the early preimplantation stage, such as fertilized eggs (zygotes), for the creation of whole genetically modified animals; the second is at post-implanted stages, as exemplified by the mid-gestational stages (E9 to E15), for targeting specific cell populations through in utero injection of viral vectors carrying genome-editing components or that of nonviral vectors carrying genome-editing components and subsequent in utero electroporation; the third is at the mid-gestational stages, as exemplified by tail-vein injection of genome-editing components into the pregnant females through which the genome-editing components can be transmitted to fetal cells via a placenta-blood barrier; and the last is at the newborn or adult stage, as exemplified by facial or tail-vein injection of genome-editing components. Here, we focus on the second and third approaches and will review the latest techniques for various methods concerning gene editing in developing fetuses.

Quantitative Trait Loci (QTL) Analysis of Seed Protein and Oil Content in Wild Soybean (Glycine soja)

Soybean seeds consist of approximately 40% protein and 20% oil, making them one of the world's most important cultivated legumes. However, the levels of these compounds are negatively correlated with each other and regulated by quantitative trait loci (QTL) that are controlled by several genes. In this study, a total of 190 F₂ and 90 BC₁F₂ plants derived from a cross of Daepung (Glycine max) with GWS-1887 (G. soja, a source of high protein), were used for the QTL analysis of protein and oil content. In the F_2:3 populations, the average protein and oil content was 45.52% and 11.59%, respectively. A QTL associated with protein levels was detected at Gm20_29512680 on chr. 20 with a likelihood of odds (LOD) of 9.57 and an R² of 17.2%. A QTL associated with oil levels was also detected at Gm15_3621773 on chr. 15 (LOD: 5.80; R²: 12.2%). In the BC₁F_2:3 populations, the average protein and oil content was 44.25% and 12.14%, respectively. A QTL associated with both protein and oil content was detected at Gm20_27578013 on chr. 20 (LOD: 3.77 and 3.06; R² 15.8% and 10.7%, respectively). The crossover to the protein content of BC₁F_3:4 population was identified by SNP marker Gm20_32603292. Based on these results, two genes, Glyma.20g088000 (S-adenosyl-l-methionine-dependent methyltransferases) and Glyma.20g088400 (oxidoreductase, 2-oxoglutarate-Fe(II) oxygenase family protein), in which the amino acid sequence had changed and a stop codon was generated due to an InDel in the exon region, were identified.

Two sides of the same coin: the N-terminal and the receptor binding domains of SARS-CoV-2 Spike

The SARS-CoV-2 Spike receptor binding domain and N-terminal domain interact with each other in an intricate mechanism. Mutations modulate the interplay between the Spike and host molecules. This editorial comments on the intricacies of SARS-CoV-2 Spike interactions.

Gender-Difference in Hair Length as Revealed by Crispr-Based Production of Long-Haired Mice with Dysfunctional FGF5 Mutations

Fibroblast growth factor 5 (FGF5) is an important molecule required for the transition from anagen to catagen phase of the mammalian hair cycle. We previously reported that Syrian hamsters harboring a 1-bp deletion in the Fgf5 gene exhibit excessive hair growth in males. Herein, we generated Fgf5 mutant mice using genome editing via oviductal nucleic acid delivery (GONAD)/improved GONAD (i-GONAD), an in vivo genome editing system used to target early embryos present in the oviductal lumen, to study gender differences in hair length in mutant mice. The two lines (Fgf5go-malc), one with a 2-bp deletion (c.552_553del) and the other with a 1-bp insertion (c.552_553insA) in exon 3 of Fgf5, were successfully established. Each mutation was predicted to disrupt a part of the FGF domain through frameshift mutation (p.Glu184ValfsX128 or p.Glu184ArgfsX128). Fgf5go-malc1 mice had heterogeneously distributed longer hairs than wild-type mice (C57BL/6J). Notably, this change was more evident in males than in females (p < 0.0001). Immunohistochemical analysis revealed the presence of FGF5 protein in the dermal papilla and outer root sheath of the hair follicles from C57BL/6J and Fgf5go-malc1 mice. Histological analysis revealed that the prolonged anagen phase might be the cause of accelerated hair growth in Fgf5go-malc1 mice.

d-StructMAn: Containerized structural annotation on the scale from genetic variants to whole proteomes

Background

Structural annotation of genetic variants in the context of intermolecular interactions and protein stability can shed light onto mechanisms of disease-related phenotypes. Three-dimensional structures of related proteins in complexes with other proteins, nucleic acids, or ligands enrich such functional interpretation, since intermolecular interactions are well conserved in evolution.

Results

We present d-StructMAn, a novel computational method that enables structural annotation of local genetic variants, such as single-nucleotide variants and in-frame indels, and implements it in a highly efficient and user-friendly tool provided as a Docker container. Using d-StructMAn, we annotated several very large sets of human genetic variants, including all variants from ClinVar and all amino acid positions in the human proteome. We were able to provide annotation for more than 46% of positions in the human proteome representing over 60% proteins.

Conclusions

d-StructMAn is the first of its kind and a highly efficient tool for structural annotation of protein-coding genetic variation in the context of observed and potential intermolecular interactions. d-StructMAn is readily applicable to proteome-scale datasets and can be an instrumental building machine-learning tool for predicting genotype-to-phenotype relationships.

Refining the Camelus dromedarius Myostatin Gene Polymorphism through Worldwide Whole-Genome Sequencing

Myostatin (MSTN) is a highly conserved negative regulator of skeletal muscle in mammals. Inactivating mutations results in a hyper-muscularity phenotype known as "double muscling" in several livestock and model species. In Camelus dromedarius, the gene structure organization and the sequence polymorphisms have been previously investigated, using Sanger and Next-Generation Sequencing technologies on a limited number of animals. Here, we carried out a follow-up study with the aim to further expand our knowledge about the sequence polymorphisms at the myostatin locus, through the whole-genome sequencing data of 183 samples representative of the geographical distribution range for this species. We focused our polymorphism analysis on the ±5 kb upstream and downstream region of the MSTN gene. A total of 99 variants (77 Single Nucleotide Polymorphisms and 22 indels) were observed. These were mainly located in intergenic and intronic regions, with only six synonymous Single Nucleotide Polymorphisms in exons. A sequence comparative analysis among the three species within the Camelus genus confirmed the expected higher genetic distance of C. dromedarius from the wild and domestic two-humped camels compared to the genetic distance between C. bactrianus and C. ferus. In silico functional prediction highlighted: (i) 213 differential putative transcription factor-binding sites, out of which 41 relative to transcription factors, with known literature evidence supporting their involvement in muscle metabolism and/or muscle development; and (ii) a number of variants potentially disrupting the canonical MSTN splicing elements, out of which two are discussed here for their potential ability to generate a prematurely truncated (inactive) form of the protein. The distribution of the considered variants in the studied cohort is discussed in light of the peculiar evolutionary history of this species and the hypothesis that extremely high muscularity, associated with a homozygous condition for mutated (inactivating) alleles at the myostatin locus, may represent, in arid desert conditions, a clear metabolic disadvantage, emphasizing the thermoregulatory and water availability challenges typical of these habitats.

Independent acquisition of short insertions at the RIR1 site in the spike N-terminal domain of the SARS-CoV-2 BA.2 lineage

Although the major SARS‐CoV‐2 omicron lineages share over 30 non‐synonymous substitutions in the spike glycoprotein, they show several unique mutations that were acquired after their ancestral split. One of the most intriguing mutations associated with BA.1 is the presence of the inserted tripeptide Glu‐Pro‐Glu within the N‐terminal domain, at a site that had previously independently acquired short insertions in several other SARS‐CoV‐2 lineages. Although the functional implications of the small nucleotide sequences found at this insertion hotspot, named RIR1, are still unclear, we have previously hypothesized that they may play a compensatory role in counterbalancing minor fitness deficits associated with other co‐occurring spike non‐synonymous mutations.

Here we show that similar insertion events have independently occurred at RIR1 at least 20 times in early 2022 within the BA.2 lineage, being occasionally associated with significant community transmission. One of these omicron sublineages, characterized by a Ser‐Gly‐Arg insertion in position 212, has been responsible of over 4,000 documented covid‐19 cases worldwide between January and July 2022, for the most part concentrated in Denmark, where it reached a national prevalence close to 4% (10% in the Nordjylland region) in mid‐May. Although the concurrent spread of the BA.2.12.1, BA.4 and BA.5 lineages led to the rapid decline of this BA.2 sublineage, the independent acquisition of several other RIR1 insertions on a BA.2 genomic background suggests that these events may provide a slight fitness advantage. Therefore, we they should be carefully monitored in the upcoming months in other emerging omicron‐related lineages, including BA.5.

This article is protected by copyright. All rights reserved

Contribution of rare mutational outcomes to broadly neutralizing antibodies

Antibodies are important immune molecules that are elicited by B cells to protect our bodies during viral infections or vaccinations. In humans, the antibody repertoire is diversified by programmed DNA lesion processes to ensure specific and high affinity binding to various antigens. Broadly neutralizing antibodies (bnAbs) are antibodies that have strong neutralizing activities against different variants of a virus. bnAbs such as anti-HIV bnAbs often have special characteristics including insertions and deletions, long complementarity determining region 3 (CDR3), and high frequencies of mutations, often at improbable sites of the variable regions. These unique features are rare mutational outcomes that are acquired during antibody diversification processes. In this review, we will discuss possible mechanisms that generate these rare antibody mutational outcomes. The understanding of the mechanisms that generate these rare mutational outcomes during antibody diversification will have implications in vaccine design strategies to elicit bnAbs.

A Comprehensive Review of Indel Detection Methods for Identification of Zebrafish Knockout Mutants Generated by Genome-Editing Nucleases

The use of zebrafish in functional genomics and disease modeling has become popular due to the ease of targeted mutagenesis with genome editing nucleases, i.e., zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and clustered regularly interspaced short palindromic repeats/Cas9 (CRISPR/Cas9). These nucleases, specifically CRISPR/Cas9, are routinely used to generate gene knockout mutants by causing a double stranded break at the desired site in the target gene and selecting for frameshift insertions or deletions (indels) caused by the errors during the repair process. Thus, a variety of methods have been developed to identify fish with indels during the process of mutant generation and phenotypic analysis. These methods range from PCR and gel-based low-throughput methods to high-throughput methods requiring specific reagents and/or equipment. Here, we provide a comprehensive review of currently used indel detection methods in zebrafish. By discussing the molecular basis for each method as well as their pros and cons, we hope that this review will serve as a comprehensive resource for zebrafish researchers, allowing them to choose the most appropriate method depending upon their budget, access to required equipment and the throughput needs of the projects.

Cas9-mediated maternal-effect and derived resistance alleles in a gene-drive strain of the African malaria vector mosquito, Anopheles gambiae

CRISPR/Cas9 technologies are important tools for the development of gene-drive systems to modify mosquito vector populations to control the transmission of pathogens that cause diseases such as malaria. However, one of the challenges for current Cas9-based drive systems is their ability to produce drive-resistant alleles resulting from insertions and deletions (indels) caused principally by nonhomologous end-joining following chromosome cleavage. Rapid increases in the frequency of such alleles may impair gene-drive dynamics. We explored the generation of indels in the germline and somatic cells in female gene-drive lineages using a series of selective crosses between a gene-drive line, AgNosCd-1, and wild-type mosquitoes. We find that potential drive-resistant mutant alleles are generated largely during embryonic development, most likely caused by deposition of the Cas9 endonuclease and guide RNAs in oocytes and resulting embryos by homozygous and hemizygous gene-drive mothers.

Allosteric Inter-Domain Contacts in Bacterial Hsp70 Are Located in Regions That Avoid Insertion and Deletion Events

Heat shock proteins 70 (Hsp70) are chaperones consisting of a nucleotide-binding domain (NBD) and a substrate-binding domain (SBD), the latter of which binds protein clients. After ATP binds to the NBD, the SBD α/β subdomains’ shared interface opens, and the open SBD docks to the NBD. Such allosteric effects are stabilized by the newly formed NBD-SBD interdomain contacts. In this paper, we examined how such an opening and formation of subdomain interfaces is affected during the evolution of Hsp70. In particular, insertion and deletion events (indels) can be highly disruptive for the mechanical events since such changes introduce a collective shift in the pairing interactions at communicating interfaces. Based on a multiple sequence alignment analysis of data collected from Swiss-Prot/UniProt database, we find several indel-free regions (IFR) in Hsp70. The two largest IFRs are located in interdomain regions that participate in allosteric structural changes. We speculate that the reason why the indels have a lower likelihood of occurrence in these regions is that indel events in these regions cause dysfunction in the protein due to perturbations of the mechanical balance. Thus, the development of functional allosteric machines requires including in the rational design a concept of the balance between structural elements.

High Polymorphism in the Dmrt2a Gene Is Incompletely Sex-Linked in Spotted Scat, Scatophagus argus

Unlike mammals and birds, many fishes have young sex chromosomes, providing excellent models to study sex chromosome differentiation at early stages. Previous studies showed that spotted scat possesses an XX-XY sex determination system. The X has a complete Dmrt3 copy (termed normal) and a truncated copy of Dmrt1 (called Dmrt1b), while the Y has the opposite (normal Dmrt1, which is male-specific, and a truncated Dmrt3 called Dmrt3△-Y). Dmrt1 is the candidate sex determination gene, while the differentiation of other sex-linked genes remains unknown. The spotted scat has proven to be a good model to study the evolution of sex chromosomes in vertebrates. Herein, we sequenced a neighbor gene of this family, Dmrt2, positioned farther from Dmrt1 and closer to Dmrt3 in the spotted scat, and analyzed its sequence variation and expression profiles. The physical locations of the three genes span across an estimated size of >40 kb. The open reading frames of Dmrt2a and its paralog Dmrt2b are 1578 bp and 1311 bp, encoding peptides of 525 and 436 amino acid residues, respectively. Dmrt2a is positioned close to Dmrt3 but farther from Dmrt1 on the same chromosome, while Dmrt2b is not. Sequence analysis revealed several mutations; insertions, and deletions (indels) on Dmrt2a non-coding regions and single-nucleotide polymorphisms (SNPs) on the Dmrt2a transcript. These indels and SNPs are sex-linked and showed high male heterogeneity but do not affect gene translation. The markers designed to span the mutation sites tested on four different populations showed varied concordance with the genetic sexes. Dmrt2a is transcribed solely in the gonads and gills, while Dmrt2b exists in the gonads, hypothalamus, gills, heart, and spleen. The Dmrt2a and Dmrt2b transcripts are profoundly expressed in the male gonads. Analyses of the transcriptome data from five other fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), Hong Kong catfish (Clarias fuscus), and spot-fin porcupine fish (Diodon hystrix)) revealed testes-biased expression of Dmrt1 in all, similar to spotted scat. Additionally, the expression of Dmrt2a is higher in the testes than the ovaries in spotted scat and Hainan medaka. The Dmrt2a transcript was not altered in the coding regions as found in Dmrt1 and Dmrt3 in spotted scat. This could be due to the functional importance of Dmrt2a in development. Another possibility is that because Dmrt2a is positioned farther from Dmrt1 and the chromosome is still young, meaning it is only a matter of time before it differentiates. This study undeniably will aid in understanding the functional divergence of the sex-linked genes in fish.

Forensic DNA phenotyping: Canis familiaris breed classification and skeletal phenotype prediction using functionally significant skeletal SNPs and indels

This review highlights a novel application of breed identification and prediction of skeletal traits in forensic investigations using canine DNA evidence. Currently, genotyping methods used for canine breed classification involve the application of highly polymorphic short tandem repeats in addition to larger commercially available SNP arrays. Both applications face technical challenges. An additional approach to breed identification could be through genotyping SNPs and indels that characterise the array of skeletal differences displayed across domestic dog populations. Research has shown that a small number of genetic variants of large effect drive differences in skeletal phenotypes among domestic dog breeds. This feature makes functionally significant canine skeletal variants a cost-effective target for forensic investigators to classify individuals according to their breed. Further analysis of these skeletal variants would enable the prediction of external appearance. To date, functionally significant genes with genetic variants associated with differences in size, bulk, skull shape, ear shape, limb length, digit type, and tail morphology have been uncovered. Recommendations of a cost-effective genotyping method that can be readily designed and applied by forensic investigators have been given. Further advances to improve the field of canine skeletal forensic DNA phenotyping include the refinement of phenotyping methods, further biological validation of the skeletal genetic variants and establishing a publicly available database for storage of allele frequencies of the skeletal genetic variants in the wider domestic dog population.

Evolutionary Dynamics of Indels in SARS-CoV-2 Spike Glycoprotein

SARS-CoV-2, responsible for the current COVID-19 pandemic that claimed over 5.0 million lives, belongs to a class of enveloped viruses that undergo quick evolutionary adjustments under selection pressure. Numerous variants have emerged in SARS-CoV-2, posing a serious challenge to the global vaccination effort and COVID-19 management. The evolutionary dynamics of this virus are only beginning to be explored. In this work, we have analysed 1.79 million spike glycoprotein sequences of SARS-CoV-2 and found that the virus is fine-tuning the spike with numerous amino acid insertions and deletions (indels). Indels seem to have a selective advantage as the proportions of sequences with indels steadily increased over time, currently at over 89%, with similar trends across countries/variants. There were as many as 420 unique indel positions and 447 unique combinations of indels. Despite their high frequency, indels resulted in only minimal alteration of N-glycosylation sites, including both gain and loss. As indels and point mutations are positively correlated and sequences with indels have significantly more point mutations, they have implications in the evolutionary dynamics of the SARS-CoV-2 spike glycoprotein.

The potential of family-free rearrangements towards gene orthology inference

Recently, we proposed an efficient ILP formulation [Rubert DP, Martinez FV, Braga MDV, Natural family-free genomic distance, Algorithms Mol Biol 16:4, 2021] for exactly computing the rearrangement distance of two genomes in a family-free setting. In such a setting, neither prior classification of genes into families, nor further restrictions on the genomes are imposed. Given two genomes, the mentioned ILP computes an optimal matching of the genes taking into account simultaneously local mutations, given by gene similarities, and large-scale genome rearrangements. Here, we explore the potential of using this ILP for inferring groups of orthologs across several species. More precisely, given a set of genomes, our method first computes all pairwise optimal gene matchings, which are then integrated into gene families in the second step. Our approach is implemented into a pipeline incorporating the pre-computation of gene similarities. It can be downloaded from gitlab.ub.uni-bielefeld.de/gi/FFGC. We obtained promising results with experiments on both simulated and real data.

Incorporating intergenic regions into reversal and transposition distances with indels

Problems in the genome rearrangement field are often formulated in terms of pairwise genome comparison: given two genomes [Formula: see text] and [Formula: see text], find the minimum number of genome rearrangements that may have occurred during the evolutionary process. This broad definition lacks at least two important considerations: the first being which features are extracted from genomes to create a useful mathematical model, and the second being which types of genome rearrangement events should be represented. Regarding the first consideration, seminal works in the genome rearrangement field solely used gene order to represent genomes as permutations of integer numbers, neglecting many important aspects like gene duplication, intergenic regions, and complex interactions between genes. Regarding the second consideration, some rearrangement events are widely studied such as reversals and transpositions. In this paper, we shed light on the first consideration and created a model that takes into account gene order and the number of nucleotides in intergenic regions. In addition, we consider events of reversals, transpositions, and indels (insertions and deletions) of genomic material. We present a 4-approximation algorithm for reversals and indels, a [Formula: see text]-approximation algorithm for transpositions and indels, and a 6-approximation for reversals, transpositions, and indels.

Labeled Cycle Graph for Transposition and Indel Distance

In the comparative genomics field, one way to infer the evolutionary distance between two organisms of related species is by finding the minimum number of large-scale mutations, called genome rearrangements, that transform one genome into the other. This number is referred to as the rearrangement distance. Since problems in this area emerged in the mid-1990s, several genome rearrangements have been proposed. Rearrangements that do not alter the genome content are called conservative, and in this group we have the following: the reversal, which inverts a segment of the genome; the transposition, which exchanges two consecutive segments; and the double cut and join, which cuts two different pairs of adjacent blocks and joins them differently. Seminal works compared genomes sharing the same set of conserved blocks, but nowadays, researchers started looking at genomes with unequal gene content, by allowing the use of nonconservative rearrangements such as insertion and deletion (jointly called indel). The transposition distance and the transposition and indel distance are both NP-hard. We investigate the transposition and indel distance and present a structure called labeled cycle graph, representing an instance of rearrangement distance problems for genomes with unequal gene content. This structure is used to devise a lower bound and a 2-approximation algorithm for the transposition and indel distance.

DNA Replication-Transcription Conflicts Do Not Significantly Contribute to Spontaneous Mutations Due to Replication Errors in Escherichia coli

Encounters between DNA replication and transcription can cause genomic disruption, particularly when the two meet head-on. Whether these conflicts produce point mutations is debated. This paper presents detailed analyses of a large collection of mutations generated during mutation accumulation experiments with mismatch repair (MMR)-defective Escherichia coli. With MMR absent, mutations are primarily due to DNA replication errors. Overall, there were no differences in the frequencies of base pair substitutions or small indels (i.e., insertion and deletions of ≤4 bp) in the coding sequences or promoters of genes oriented codirectionally versus head-on to replication. Among a subset of highly expressed genes, there was a 2- to 3-fold bias for indels in genes oriented head-on to replication, but this difference was almost entirely due to the asymmetrical genomic locations of tRNA genes containing mononucleotide runs, which are hot spots for indels. No additional orientation bias in mutation frequencies occurred when MMR⁻ strains were also defective for transcription-coupled repair (TCR). However, in contrast to other reports, loss of TCR slightly increased the overall mutation rate, meaning that TCR is antimutagenic. There was no orientation bias in mutation frequencies among the stress response genes that are regulated by RpoS or induced by DNA damage. Thus, biases in the locations of mutational targets can account for most, if not all, apparent biases in mutation frequencies between genes oriented head-on versus codirectional to replication. In addition, the data revealed a strong correlation of the frequency of base pair substitutions with gene length but no correlation with gene expression levels.

A probabilistic model for indel evolution: differentiating insertions from deletions

Insertions and deletions (indels) are common molecular evolutionary events. However, probabilistic models for indel evolution are under-developed due to their computational complexity. Here, we introduce several improvements to indel modeling: 1) While previous models for indel evolution assumed that the rates and length distributions of insertions and deletions are equal, here we propose a richer model that explicitly distinguishes between the two; 2) we introduce numerous summary statistics that allow approximate Bayesian computation-based parameter estimation; 3) we develop a method to correct for biases introduced by alignment programs, when inferring indel parameters from empirical data sets; and 4) using a model-selection scheme, we test whether the richer model better fits biological data compared with the simpler model. Our analyses suggest that both our inference scheme and the model-selection procedure achieve high accuracy on simulated data. We further demonstrate that our proposed richer model better fits a large number of empirical data sets and that, for the majority of these data sets, the deletion rate is higher than the insertion rate.

MSI-WES: a simple approach for microsatellite instability testing using whole exome sequencing

Aim: To demonstrate that MSI-WES is an accurate testing method for microsatellite instability (MSI). Materials & methods: Microsatellite-based indels were counted in the variant call-formatted whole exome sequencing (WES) data of 441 gastric cancer cases using Unix-based algorithms, and the counts expressed as a fraction of the genome sequenced to obtain next-generation sequencing-based MSI indices. Results: The next-generation sequencing-based MSI indices showed a near-perfect concordance with PCR-based MSI status, and moderate to good correlations with the molecular targets of MSI index, MLH1 expression and MLH1 methylation status, at a level comparable to the strengths of correlation between PCR-based MSI status and molecular targets of MSI index/MLH1 expression and methylation. Conclusion: MSI-WES is a valid, adequate and sensitive approach for testing MSI in cancer.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

Molecular Characterization of the Coat Protein Gene of Greek Apple Stem Pitting Virus Isolates: Evolution through Deletions, Insertions, and Recombination Events

A RT–PCR assay developed to amplify the full coat protein (CP) gene of apple stem pitting virus (ASPV) was evaluated using 180 Greek apple and pear samples and showed a broad detection range. This method was used to investigate the presence of ASPV in quince in Greece and showed a high incidence of 52%. The sequences of 14 isolates from various hosts with a distinct RFLP profile were determined. ASPV population genetics and the factors driving ASPV evolution were analyzed using the Greek ASPV sequences, novel sequences from Brazilian apple trees and Chinese botanical Pyrus species, and homologous sequences retrieved from GenBank. Fourteen variant types of Greek, Brazilian and botanical isolates, which differ in CP gene length and presence of indels, were identified. In addition, these analyses showed high intra- and inter-group variation among isolates from different countries and hosts, indicating the significant variability present in ASPV. Recombination events were detected in four isolates originating from Greek pear and quince and two from Brazilian apples. In a phylogenetic analysis, there was a tendency for isolates to cluster together based on CP gene length, the isolation host, and the detection method applied. Although there was no strict clustering based on geographical origin, most isolates from a given country tended to regroup in specific clusters. Interestingly, it was found that the phylogeny was correlated to the type, position, and pattern of indels, which represent hallmarks of specific lineages and indicate their possible role in virus diversification, rather than the CP size itself. Evidence of recombination between isolates from botanical and cultivated species and the clustering of isolates from botanical species and isolates from cultivated species suggest the existence of a possible undetermined transmission mechanism allowing the exchange of ASPV isolates between the cultivated and wild/ornamental hosts.

The tree of life describes a tripartite cellular world

The canonical view of a 3-domain (3D) tree of life was recently challenged by the discovery of Asgardarchaeota encoding eukaryote signature proteins (ESPs), which were treated as missing links of a 2-domain (2D) tree. Here we revisit the debate. We discuss methodological limitations of building trees with alignment-dependent approaches, which often fail to satisfactorily address the problem of ''gaps.'' In addition, most phylogenies are reconstructed unrooted, neglecting the power of direct rooting methods. Alignment-free methodologies lift most difficulties but require employing realistic evolutionary models. We argue that the discoveries of Asgards and ESPs, by themselves, do not rule out the 3D tree, which is strongly supported by comparative and evolutionary genomic analyses and vast genomic and biochemical superkingdom distinctions. Given uncertainties of retrodiction and interpretation difficulties, we conclude that the 3D view has not been falsified but instead has been strengthened by genomic analyses. In turn, the objections to the 2D model have not been lifted. The debate remains open. Also see the video abstract here: https://youtu.be/-6TBN0bubI8.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

Genome diversity in Ukraine

Background

The main goal of this collaborative effort is to provide genome-wide data for the previously underrepresented population in Eastern Europe, and to provide cross-validation of the data from genome sequences and genotypes of the same individuals acquired by different technologies. We collected 97 genome-grade DNA samples from consented individuals representing major regions of Ukraine that were consented for public data release. BGISEQ-500 sequence data and genotypes by an Illumina GWAS chip were cross-validated on multiple samples and additionally referenced to 1 sample that has been resequenced by Illumina NovaSeq6000 S4 at high coverage.

Results

The genome data have been searched for genomic variation represented in this population, and a number of variants have been reported: large structural variants, indels, copy number variations, single-nucletide polymorphisms, and microsatellites. To our knowledge, this study provides the largest to-date survey of genetic variation in Ukraine, creating a public reference resource aiming to provide data for medical research in a large understudied population.

Conclusions

Our results indicate that the genetic diversity of the Ukrainian population is uniquely shaped by evolutionary and demographic forces and cannot be ignored in future genetic and biomedical studies. These data will contribute a wealth of new information bringing forth a wealth of novel, endemic and medically related alleles.

A Model of Indel Evolution by Finite-State, Continuous-Time Machines

We introduce a systematic method of approximating finite-time transition probabilities for continuous-time insertion-deletion models on sequences. The method uses automata theory to describe the action of an infinitesimal evolutionary generator on a probability distribution over alignments, where both the generator and the alignment distribution can be represented by pair hidden Markov models (HMMs). In general, combining HMMs in this way induces a multiplication of their state spaces; to control this, we introduce a coarse-graining operation to keep the state space at a constant size. This leads naturally to ordinary differential equations for the evolution of the transition probabilities of the approximating pair HMM. The TKF91 model emerges as an exact solution to these equations for the special case of single-residue indels. For the more general case of multiple-residue indels, the equations can be solved by numerical integration. Using simulated data, we show that the resulting distribution over alignments, when compared to previous approximations, is a better fit over a broader range of parameters. We also propose a related approach to develop differential equations for sufficient statistics to estimate the underlying instantaneous indel rates by expectation maximization. Our code and data are available at https://github.com/ihh/trajectory-likelihood.

Effect of indel variants within the sorting nexin 29 (SNX29) gene on growth traits of goats

The sorting nexin 29 gene (SNX29) is a well-known regulator of myocyte differentiation and proliferation. In this work, two indels (17-bp and 21-bp) were identified in the goat SNX29 gene, and their effects on the growth traits of 1,759 Shaanbei white cashmere (SBWC) goats were analyzed. Both indels had three genotypes [homozygote wild type (II), heterozygote (ID), and homozygote mutation (DD)] and displayed medium genetic diversity (0.25 < polymorphism information content (PIC) < 0.50) in the population. The 17-bp indel was significantly associated with chest width (p = 0.009), body weight (p = 0.021), and chest depth (p = 0.032), with the II genotype dominant. The 21-bp indel was significantly associated with chest width (p = 0.001), chest depth (p = 4.8E-5), heart girth (p = 0.007), and hip width (p = 0.002). Because the two indels were in the upstream (17-bp) and intron (21-bp) regions of the SNX29 gene, transcription factor binding sites were predicted. The IRF5 and MYC could bind with the 17-bp indel and 21-bp indel sequences, respectively. This study indicates that SNX29 is a promising candidate gene that can be used to improve meat production in goat breeding.

Genome Rearrangement Distance with Reversals, Transpositions, and Indels

The rearrangement distance is a well-known problem in the field of comparative genomics. Given two genomes, the rearrangement distance is the minimum number of rearrangements in a set of allowed rearrangements (rearrangement model), which transforms one genome into the other. In rearrangement distance problems, a genome is modeled as a string, where each element represents a conserved region within the two genomes. When the orientation of the genes is known, it is represented by (plus or minus) signs assigned to the elements of the string. Two of the most studied rearrangements are reversals, which invert a segment of the genome, and transpositions, which exchange the relative positions of two adjacent segments of the genome. The first works in genome rearrangements considered that the genomes being compared had the same genetic material and that rearrangement events were restricted to reversals, transpositions, or both. El-Mabrouk extended the reversal model on signed strings to include the operations of insertion and deletion of segments in the genome, which allowed the comparison of genomes with different genetic material. Other studies also addressed this problem and, recently, this problem was proved to be solvable in polynomial time by Willing et al. For unsigned strings, we still observe a lack of results. That said, in this study we prove that computing the rearrangement distance for the following models is NP-Hard: reversals and indels on unsigned strings; transpositions and indels on unsigned strings; and reversals, transpositions, and indels on signed and unsigned strings. Along with the NP-hardness proofs, we present a 2-approximation algorithm for reversals on unsigned strings and 3-approximation algorithms for the other models.

Comparative Plastomics of Ashwagandha (Withania, Solanaceae) and Identification of Mutational Hotspots for Barcoding Medicinal Plants

Within the family Solanaceae, Withania is a small genus belonging to the Solanoideae subfamily. Here, we report the de novo assembled chloroplast genome sequences of W. coagulans, W. adpressa, and W. riebeckii. The length of these genomes ranged from 154,162 to 154,364 base pairs (bp). These genomes contained a pair of inverted repeats (IRa and IRb) ranging from 25,029 to 25,071 bp that were separated by a large single-copy (LSC) region of 85,635-85,765 bp and a small single-copy (SSC) region of 18,457-18,469 bp. We analyzed the structural organization, gene content and order, guanine-cytosine content, codon usage, RNA-editing sites, microsatellites, oligonucleotide and tandem repeats, and substitutions of Withania plastomes, which revealed high similarities among the species. Comparative analysis among the Withania species also highlighted 10 divergent hotspots that could potentially be used for molecular marker development, phylogenetic analysis, and species identification. Furthermore, our analyses showed that even three mutational hotspots (rps4-trnT, trnM-atpE, and rps15) were sufficient to discriminate the Withania species included in current study.

First isolation and characterization of Brucella suis from yak

Brucella spp., facultative intracellular pathogens that can persistently colonize animal host cells and cause zoonosis, affect public health and safety. A Brucella strain was isolated from yak in Qinghai Province. To detect whether this isolate could cause an outbreak of brucellosis and to reveal its genetic characteristics, several typing and whole-genome sequencing methods were applied to identify its species and genetic characteristics. Phylogenetic analysis based on MLVA and whole-genome sequencing revealed the genetic characteristics of the isolated strain. The results showed that the isolated strain is a B. suis biovar 1 smooth strain, and this isolate was named B. suis QH05. The results of comparative genomics and MLVA showed that B. suis QH05 is not a vaccine strain. Comparison with other B. suis strains isolated from humans and animals indicated that B. suis QH05 may be linked to specific animal and human sources. In conclusion, B. suis QH05 does not belong to the Brucella epidemic species in China, and as the first isolation of B. suis from yak, this strain expands the host range of B. suis.

Identification of Structural Variants in Two Novel Genomes of Maize Inbred Lines Possibly Related to Glyphosate Tolerance

To study genetic variations between genomes of plants that are naturally tolerant and sensitive to glyphosate, we used two Zea mays L. lines traditionally bred in Poland. To overcome the complexity of the maize genome, two sequencing technologies were employed: Illumina and Single Molecule Real-Time (SMRT) PacBio. Eleven thousand structural variants, 4 million SNPs and approximately 800 thousand indels differentiating the two genomes were identified. Detailed analyses allowed to identify 20 variations within the EPSPS gene, but all of them were predicted to have moderate or unknown effects on gene expression. Other genes of the shikimate pathway encoding bifunctional 3-dehydroquinate dehydratase/shikimate dehydrogenase and chorismate synthase were altered by variants predicted to have a high impact on gene expression. Additionally, high-impact variants located within the genes involved in the active transport of glyphosate through the cell membrane encoding phosphate transporters as well as multidrug and toxic compound extrusion have been identified.

Sequential i-GONAD: An Improved In Vivo Technique for CRISPR/Cas9-Based Genetic Manipulations in Mice

Improved genome-editing via oviductal nucleic acid delivery (i-GONAD) is a technique capable of inducing genomic changes in preimplantation embryos (zygotes) present within the oviduct of a pregnant female. i-GONAD involves intraoviductal injection of a solution containing genome-editing components via a glass micropipette under a dissecting microscope, followed by in vivo electroporation using tweezer-type electrodes. i-GONAD does not involve ex vivo handling of embryos (isolation of zygotes, microinjection or electroporation of zygotes, and egg transfer of the treated embryos to the oviducts of a recipient female), which is required for in vitro genome-editing of zygotes. i-GONAD enables the generation of indels, knock-in (KI) of ~ 1 kb sequence of interest, and large deletion at a target locus. i-GONAD is usually performed on Day 0.7 of pregnancy, which corresponds to the late zygote stage. During the initial development of this technique, we performed i-GONAD on Days 1.4–1.5 (corresponding to the 2-cell stage). Theoretically, this means that at least two GONAD steps (on Day 0.7 and Day 1.4–1.5) must be performed. If this is practically demonstrated, it provides additional options for various clustered regularly interspaced palindrome repeats (CRISPR)/Caspase 9 (Cas9)-based genetic manipulations. For example, it is usually difficult to induce two independent indels at the target sites, which are located very close to each other, by simultaneous transfection of two guide RNAs and Cas9 protein. However, the sequential induction of indels at a target site may be possible when repeated i-GONAD is performed on different days. Furthermore, simultaneous introduction of two mutated lox sites (to which Cre recombinase bind) for making a floxed allele is reported to be difficult, as it often causes deletion of a sequence between the two gRNA target sites. However, differential KI of lox sites may be possible when repeated i-GONAD is performed on different days. In this study, we performed proof-of-principle experiments to demonstrate the feasibility of the proposed approach called “sequential i-GONAD (si-GONAD).”

Metagenomics of Imported Multidrug-Resistant Mycobacterium leprae, Saudi Arabia, 2017

Using shotgun metagenomics, we identified an imported case of multidrug-resistant Mycobacterium leprae in a Filipino resident of Saudi Arabia in 2017. We determined the phylogenomic lineage (3K1) and identified mutations in rpoB and rrs corresponding to the multidrug-resistance phenotype clinically observed. Metagenomics sequencing can be used to identify multidrug-resistant M. leprae.

Complete Chloroplast Genomes of Three Medicinal Alpinia Species: Genome Organization, Comparative Analyses and Phylogenetic Relationships in Family Zingiberaceae

Alpinia katsumadai (A. katsumadai), Alpinia oxyphylla (A. oxyphylla) and Alpinia pumila (A. pumila), which belong to the family Zingiberaceae, exhibit multiple medicinal properties. The chloroplast genome of a non-model plant provides valuable information for species identification and phylogenetic analysis. Here, we sequenced three complete chloroplast genomes of A. katsumadai, A. oxyphylla sampled from Guangdong and A. pumila, and analyzed the published chloroplast genomes of Alpinia zerumbet (A. zerumbet) and A. oxyphylla sampled from Hainan to retrieve useful chloroplast molecular resources for Alpinia. The five Alpinia chloroplast genomes possessed typical quadripartite structures comprising of a large single copy (LSC, 87,248-87,667 bp), a small single copy (SSC, 15,306-18,295 bp) and a pair of inverted repeats (IR, 26,917-29,707 bp). They had similar gene contents, gene orders and GC contents, but were slightly different in the numbers of small sequence repeats (SSRs) and long repeats. Interestingly, fifteen highly divergent regions (rpl36, ycf1, rps15, rpl22, infA, psbT-psbN, accD-psaI, petD-rpoA, psaC-ndhE, ccsA-ndhD, ndhF-rpl32, rps11-rpl36, infA-rps8, psbC-psbZ, and rpl32-ccsA), which could be suitable for species identification and phylogenetic studies, were detected in the Alpinia chloroplast genomes. Comparative analyses among the five chloroplast genomes indicated that 1891 mutational events, including 304 single nucleotide polymorphisms (SNPs) and 118 insertion/deletions (indels) between A. pumila and A. katsumadai, 367 SNPs and 122 indels between A. pumila and A. oxyphylla sampled from Guangdong, 331 SNPs and 115 indels between A. pumila and A. zerumbet, 371 SNPs and 120 indels between A. pumila and A. oxyphylla sampled from Hainan, and 20 SNPs and 23 indels between the two accessions of A. oxyphylla, were accurately located. Additionally, phylogenetic relationships based on SNP matrix among 28 whole chloroplast genomes showed that Alpinia was a sister branch to Amomum in the family Zingiberaceae, and that the five Alpinia accessions were divided into three groups, one including A. pumila, another including A. zerumbet and A. katsumadai, and the other including two accessions of A. oxyphylla. In conclusion, the complete chloroplast genomes of the three medicinal Alpinia species in this study provided valuable genomic resources for further phylogeny and species identification in the family Zingiberaceae.

A Single "All-in-One" Helper-Dependent Adenovirus to Deliver Donor DNA and CRISPR/Cas9 for Efficient Homology-Directed Repair

In this study, we developed a single helper-dependent adenovirus (HDAd) to deliver all of the components (donor DNA, CRISPR-associated protein 9 [Cas9], and guide RNA [gRNA]) needed to achieve high-efficiency gene targeting and homology-directed repair in transduced cells. We show that these "all-in-one" HDAds are up to 117-fold more efficient at gene targeting than donor HDAds that do not express CRISPR/Cas9 in human induced pluripotent stem cells (iPSCs). The vast majority (>90%) of targeted recombinants had only one allele targeted, and this was accompanied by high-frequency indel formation in the non-targeted allele at the site of Cas9 cleavage. These indels varied in size and nature, and included large deletions of ∼8 kb. The remaining minority of recombinants had both alleles targeted (so-called bi-allelic targeting). These all-in-one HDAds represent an important platform for accomplishing and expanding the utility of homology-directed repair, especially for difficult-to-transfect cells and for in vivo applications.

Elucidation of de novo small insertion/deletion biology with parent-of-origin phasing

The mechanisms underlying de novo insertion/deletion (indel) genesis, such as polymerase slippage, have been hypothesized but not well characterized in the human genome. We implemented two methodological improvements, which were leveraged to dissect indel mutagenesis. We assigned de novo variants to parent-of-origin (i.e., phasing) with low-coverage long-read whole-genome sequencing, achieving better phasing compared to short-read sequencing (medians of 84% and 23%, respectively). We then wrote an application programming interface to classify indels into three subtypes according to sequence context. Across three cohorts with different phasing methods (N_trios  = 540, all cohorts), we observed that one de novo indel subtype, change in copy count (CCC), was significantly correlated with father's (p = 7.1 × 10^-4 ) but not mother's (p = .45) age at conception. We replicated this effect in three cohorts without de novo phasing (p_paternal  = 1.9 × 10^-9 , p_maternal  = .61; N_trios  = 3,391, all cohorts). Although this is consistent with polymerase slippage during spermatogenesis, the percentage of variance explained by paternal age was low, and we did not observe an association with replication timing. These results suggest that spermatogenesis-specific events have a minor role in CCC indel mutagenesis, one not observed for other indel subtypes nor for maternal age in general. These results have implications for indel modeling in evolution and disease.

Transposons played a major role in the diversification between the closely related almond and peach genomes: results from the almond genome sequence

We sequenced the genome of the highly heterozygous almond Prunus dulcis cv. Texas combining short- and long-read sequencing. We obtained a genome assembly totaling 227.6 Mb of the estimated almond genome size of 238 Mb, of which 91% is anchored to eight pseudomolecules corresponding to its haploid chromosome complement, and annotated 27 969 protein-coding genes and 6747 non-coding transcripts. By phylogenomic comparison with the genomes of 16 additional close and distant species we estimated that almond and peach (Prunus persica) diverged around 5.88 million years ago. These two genomes are highly syntenic and show a high degree of sequence conservation (20 nucleotide substitutions per kb). However, they also exhibit a high number of presence/absence variants, many attributable to the movement of transposable elements (TEs). Transposable elements have generated an important number of presence/absence variants between almond and peach, and we show that the recent history of TE movement seems markedly different between them. Transposable elements may also be at the origin of important phenotypic differences between both species, and in particular for the sweet kernel phenotype, a key agronomic and domestication character for almond. Here we show that in sweet almond cultivars, highly methylated TE insertions surround a gene involved in the biosynthesis of amygdalin, whose reduced expression has been correlated with the sweet almond phenotype. Altogether, our results suggest a key role of TEs in the recent history and diversification of almond and its close relative peach.

The genomic pool of standing structural variation outnumbers single nucleotide polymorphism by threefold in the marine teleost Chrysophrys auratus

Recent studies have highlighted an important role of structural variation (SV) in ecological and evolutionary processes, but few have studied nonmodel species in the wild. As part of our long-term research programme on the nonmodel teleost fish Australasian snapper (Chrysophrys auratus), we aim to build one of the first catalogues of genomic variants (SNPs and indels, and deletions, duplications and inversions) in fishes and evaluate overlap of genomic variants with regions under putative selection (Tajima's D and π), and coding sequences (genes). For this, we analysed six males and six females from three locations in New Zealand and generated a high-resolution genomic variation catalogue. We characterized 20,385 SVs and found they intersected with almost a third of all annotated genes. Together with small indels, SVs account for three times more variation in the genome in terms of bases affected compared to SNPs. We found that a sizeable portion of detected SVs was in the upper and lower genomic regions of Tajima's D and π, indicating that some of these have an effect on the phenotype. Together, these results shed light on the often neglected genomic variation that is produced by SVs and highlights the need to go beyond the mere measure of SNPs when investigating evolutionary processes, such as species diversification and adaptation.

Accuracy of de novo assembly of DNA sequences from double-digest libraries varies substantially among software

Advances in DNA sequencing have made it feasible to gather genomic data for non-model organisms and large sets of individuals, often using methods for sequencing subsets of the genome. Several of these methods sequence DNA associated with endonuclease restriction sites (various RAD and GBS methods). For use in taxa without a reference genome, these methods rely on de novo assembly of fragments in the sequencing library. Many of the software options available for this application were originally developed for other assembly types and we do not know their accuracy for reduced representation libraries. To address this important knowledge gap, we simulated data from the Arabidopsis thaliana and Homo sapiens genomes and compared de novo assemblies by six software programs that are commonly used or promising for this purpose (ABySS, CD-HIT, Stacks, Stacks2, Velvet and VSEARCH). We simulated different mutation rates and types of mutations, and then applied the six assemblers to the simulated data sets, varying assembly parameters. We found substantial variation in software performance across simulations and parameter settings. ABySS failed to recover any true genome fragments, and Velvet and VSEARCH performed poorly for most simulations. Stacks and Stacks2 produced accurate assemblies of simulations containing SNPs, but the addition of insertion and deletion mutations decreased their performance. CD-HIT was the only assembler that consistently recovered a high proportion of true genome fragments. Here, we demonstrate the substantial difference in the accuracy of assemblies from different software programs and the importance of comparing assemblies that result from different parameter settings.

NHEJ-Mediated Repair of CRISPR-Cas9-Induced DNA Breaks Efficiently Corrects Mutations in HSPCs from Patients with Fanconi Anemia

Non-homologous end-joining (NHEJ) is the preferred mechanism used by hematopoietic stem cells (HSCs) to repair double-stranded DNA breaks and is particularly increased in cells deficient in the Fanconi anemia (FA) pathway. Here, we show feasible correction of compromised functional phenotypes in hematopoietic cells from multiple FA complementation groups, including FA-A, FA-C, FA-D1, and FA-D2. NHEJ-mediated repair of targeted CRISPR-Cas9-induced DNA breaks generated compensatory insertions and deletions that restore the coding frame of the mutated gene. NHEJ-mediated editing efficacy was initially verified in FA lymphoblastic cell lines and then in primary FA patient-derived CD34⁺ cells, which showed marked proliferative advantage and phenotypic correction both in vitro and after transplantation. Importantly, and in contrast to homologous directed repair, NHEJ efficiently targeted primitive human HSCs, indicating that NHEJ editing approaches may constitute a sound alternative for editing self-renewing human HSCs and consequently for treatment of FA and other monogenic diseases affecting the hematopoietic system.

Transplacental delivery of genome editing components causes mutations in embryonic cardiomyocytes of mid-gestational murine fetuses

Genome editing, as exemplified by CRISPR/Cas9, is now recognized as a powerful tool for the engineering of endogenous target genes. It employs only two components, namely, Cas9 in the form of DNA, mRNA, or protein; and guide RNA (gRNA), which is specific to a target gene. When these components are transferred to cells, they create insertion/deletion mutations (indels) within a target gene. Therefore, when fetuses within the uteri of pregnant murine females are exposed to these reagents, fetal cells incorporating them should show mutations in the target gene. To examine a possible genome editing of fetal cells in vivo, we intravenously administered a solution containing plasmid DNA-FuGENE complex to pregnant wild-type female mice [which had been successfully mated with enhanced green fluorescent protein (EGFP)-expressing male transgenic mice] on day 12.5 of gestation. The plasmid DNA induces the expression of gRNA, which was targeted at the EGFP cDNA, and that of the Cas9 gene. All fetuses in the pregnant females should express EGFP systemically, since they are heterozygous (Tg/+) for the transgene. Thus, the delivery of CRISPR system targeted at EGFP in the fetuses will cause a reduced expression of EGFP as a result of the genome editing of EGFP genomic sequence. Of the 24 fetuses isolated from three pregnant females 2 days after gene delivery, 3 were found to have reduced fluorescence in their hearts. Genotyping of the dissected hearts revealed the presence of the transgene construct (Cas9 gene) in all the samples. Furthermore, all the three samples exhibited mutations at the target loci, although normal cells were also present. Thus, transplacental delivery of gene editing components may be a useful tool for developing animal models with heart disorder for heart-related disease research, and gene therapy in congenital heart defects such as hypertrophic cardiomyopathy (HCM). © 2019 IUBMB Life, 9999(9999):1-10, 2019.

Linearization of Median Genomes Under the Double-Cut-and-Join-Indel Model

Reconstruction of the median genome consisting of linear chromosomes from three given genomes is known to be intractable. There exist efficient methods for solving a relaxed version of this problem, where the median genome is allowed to have circular chromosomes. We propose a method for construction of an approximate solution to the original problem from a solution to the relaxed problem and prove a bound on its approximation error. Our method also provides insights into the combinatorial structure of genome transformations with respect to appearance of circular chromosomes.

Never the Two Shall Mix: Robust Indel Markers to Ensure the Fidelity of Two Pivotal and Closely-Related Accessions of Brachypodium distachyon

Brachypodium distachyon is an established model for monocotyledonous plants. Numerous markers intended for gene discovery and population genetics have been designed. However to date, very few indel markers with larger and easily scored length polymorphism differences, that distinguish between the two morphologically similar and highly utilized B. distachyon accessions, Bd21, the reference genome accession, and Bd21-3, the transformation-optimal accession, are publically available. In this study, 22 indel markers were designed and utilized to produce length polymorphism differences of 150 bp or more, for easy discrimination between Bd21 and Bd21-3. When tested on four other B. distachyon accessions, one case of multiallelism was observed. It was also shown that the markers could be used to determine homozygosity and heterozygosity at specific loci in a Bd21 x Bd3-1 F2 population. The work done in this study allows researchers to maintain the fidelity of Bd21 and Bd21-3 stocks for both transgenic and nontransgenic studies. It also provides markers that can be utilized in conjunction with others already available for further research on population genetics, gene discovery and gene characterization, all of which are necessary for the relevance of B. distachyon as a model species.

Prognostic implications of a molecular classifier derived from whole-exome sequencing in nasopharyngeal carcinoma

The aim of this study was to use whole-exome sequencing to derive a molecular classifier for nasopharyngeal carcinoma (NPC) and evaluate its clinical performance. We performed whole-exome sequencing on 82 primary NPC tumors from Sun Yat-sen University Cancer Center (Guangzhou cohort) to obtain somatic single-nucleotide variants, indels, and copy number variants. A novel molecular classifier was then developed and validated in another NPC cohort (Hong Kong cohort, n = 99). Survival analysis was estimated by the Kaplan-Meier method and compared using the log-rank test. Cox proportional hazards model was adopted for univariate and multivariate analyses. We identified three prominent NPC genetic subtypes: RAS/PI3K/AKT (based on RAS, AKT1, and PIK3CA mutations), cell-cycle (based on CDKN2A/CDKN2B deletions, and CDKN1B and CCND1 amplifications), and unclassified (based on dominant mutations in epigenetic regulators, such as KMT2C/2D, or the Notch signaling pathway, such as NOTCH1/2). These subtypes differed in survival analysis, with good, intermediate, and poor progression-free survival in the unclassified, cell-cycle, and RAS/PI3K/AKT subgroups, respectively, among the Guangzhou, Hong Kong, and combined cohorts (n = 82, P = 0.0342; n = 99, P = 0.0372; and n = 181, P = 0.0023; log-rank test). We have uncovered genetic subtypes of NPC with distinct mutations and/or copy number changes, reflecting discrete paths of NPC tumorigenesis and providing a roadmap for developing new prognostic biomarkers and targeted therapies.

Production of CRISPR/Cas9-Mediated Self-Cleaving Helper-Dependent Adenoviruses

Prolonged expression of CRISPR/Cas9 raises concerns about off-target cleavage, cytotoxicity, and immune responses. To address these issues, we have developed a system to produce helper-dependent adenoviruses that express CRISPR/Cas9 to direct cleavage of the vectors’ own genome after transduction of target cells. To prevent self-cleavage during vector production, it was necessary to downregulate Cas9 mRNA as well as inhibit Cas9 protein activity. Cas9 mRNA downregulation was achieved by inserting the target sequences for the helper-virus-encoded miRNA, mivaRNAI, and producer-cell-encoded miRNAs, hsa-miR183-5p, and hsa-miR218-5p, into the 3′ UTR of the HDAd-encoded Cas9 expression cassette. Cas9 protein activity was inhibited by expressing anti-CRISPR proteins AcrIIA2 and AcrAII4 from both the producer cells and the helper virus. After purification, these helper-dependent adenoviruses will perform CRISPR/Cas9-mediated self-cleavage in the transduced target cells, thereby limiting the duration of Cas9 expression and thus represent an important platform for improving the safety of gene editing by CRISPR/Cas9.