A robust benchmark for evaluating and improving mosaic variant calling strategies

The rapid advances in sequencing and analysis technologies have enabled the accurate detection of diverse forms of genomic variants, including germline, somatic, and mosaic mutations. However, unlike for the former two mutations, the best practices for mosaic variant calling still remain chaotic due to the technical and conceptual difficulties faced in evaluation. Here, we present our benchmark of nine feasible strategies for mosaic variant detection based on a systematically designed reference standard that mimics mosaic samples, with 390,153 control positive and 35,208,888 negative single-nucleotide variants and insertion–deletion mutations. We identified the condition-dependent strengths and weaknesses of the current strategies, instead of a single winner, regarding variant allele frequencies, variant sharing, and the usage of control samples. Moreover, feature-level investigation directs the way for immediate to prolonged improvements in mosaic variant calling. Our results will guide researchers in selecting suitable calling algorithms and suggest future strategies for developers.

Characterization of Genomic Variation from Lotus (Nelumbo Adans.) Mutants with Wide and Narrow Tepals

Compared with rose, chrysanthemum, and water lily, the absence of short-wide and long-narrow tepals of ornamental lotus (Nelumbo Adans.) limits the commercial value of flowers. In this study, the genomes of two groups of lotus mutants with wide-short and narrow-long tepals were resequenced to uncover the genomic variation and candidate genes associated with tepal shape. In group NL (short for N. lutea, containing two mutants and one control of N. lutea), 716,656 single nucleotide polymorphisms (SNPs) and 221,688 insertion-deletion mutations (Indels) were obtained, while 639,953 SNPs and 134,6118 Indels were obtained in group WSH (short for ‘Weishan Hong’, containing one mutant and two controls of N. nucifera ‘Weishan Hong’). Only a small proportion of these SNPs and Indels was mapped to exonic regions of genome: 1.92% and 0.47%, respectively, in the NL group, and 1.66% and 0.48%, respectively, in the WSH group. Gene Ontology (GO) analysis showed that out of 4890 (NL group) and 1272 (WSH group) annotated variant genes, 125 and 62 genes were enriched (Q < 0.05), respectively. Additionally, in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, 104 genes (NL group) and 35 genes (WSH group) were selected (p < 0.05). Finally, there were 306 candidate genes that were sieved to determine the development of tepal shape in lotus plants. It will be an essential reference for future identification of tepal-shaped control genes in lotus plants. This is the first comprehensive report of genomic variation controlling tepal shape in lotus, and the mutants in this study are promising materials for breeding novel lotus cultivars with special tepals.

Phage Selective Pressure Reduces Virulence of Hypervirulent Klebsiella pneumoniae Through Mutation of the wzc Gene

Hypervirulent Klebsiella pneumoniae (hvKp), one of the major community-acquired pathogens, can cause invasive infections such as liver abscess. In recent years, bacteriophages have been used in the treatment of K. pneumoniae, but the characteristics of the phage-resistant bacteria produced in the process of phage therapy need to be evaluated. In this study, two Podoviridae phages, hvKpP1 and hvKpP2, were isolated and characterized. In vitro and in vivo experiments demonstrated that the virulence of the resistant bacteria was significantly reduced compared with that of the wild type. Comparative genomic analysis of monoclonal sequencing showed that nucleotide deletion mutations of wzc and wcaJ genes led to phage resistance, and the electron microscopy and mucoviscosity results showed that mutations led to the loss of the capsule. Meanwhile, animal assay indicated that loss of capsule reduced the virulence of hvKp. These findings contribute to a better understanding of bacteriophage therapy, which not only can kill bacteria directly but also can reduce the virulence of bacteria by phage screening.

Multiplex Fragment analysis detects all COVID-19 variants of concern

Introduction/Objective The majority of tracking methods have employed whole genome sequencing, which can be very expensive and time consuming. An alternative method has been to use genotyping of specific mutations to identify variants. However, tracking SARS-CoV-2 variants by targeted methods has been a moving target. Most methods only multiplex four targets per reaction, but we have multiplexed 8 targets in a single tube using fragment analysis. Methods/Case Report Fluorescently labeled primers targeted a combination of insertion/ deletion mutations and single nucleotide mutations. The PCR amplified products, amplicons, were separated by capillary electrophoresis. Primers were designed to detect changes in size indicative of insertion or deletion mutations including: ORF1A:Del3675_3677, S:Del69_70, S:Del144, S:Del157_158, S:Del242_244, ORF8:Del119_120, and ORF8:ins28269-28273. Allele-specific primers were designed to detect both the wild-type and mutated versions of S:N501Y, S:E484K, and S:L452R. Residual nasopharyngeal and nasal specimens testing positive for SARS-CoV-2 by RT-PCR or isothermal amplification (IDnow) methods were selected from May 1- June 24, 2021. Variant analysis was performed by multiplex targeted PCR and whole genome sequencing in parallel on the same specimens to determine positive percent agreement. Results (if a Case Study enter NA) Variant analysis was performed on 250 specimens detecting each of the major variants of concern Alpha (B.1.1.7, U.K. origin, n= 108), Beta (B.1.351, South Africa origin, n=3), Gamma (P.1, Brazil origin, n=12), Delta (B.1.617.2, Indian origin, n=17), and Iota (B.1.526, New York, n=5). Some specimens with low viral load were detected by only PCR (n=18), only WGS (n=41), or neither (n=20). Overall positive percent agreement was 95% (163/171). Conclusion This adjustable method robustly and accurately identifies COVID-19 VOCs utilizing a platform amenable to multiple targets (20-40 targets ranging from 100-500b.p. across four fluorescent channels) using equipment commonly found in routine molecular pathology laboratories. Future directions include adjusting targets to detect new variants.

A 104-bp Structural Variation of the ADPRHL1 Gene Is Associated With Growth Traits in Chickens

Analyzing marker-assisted breeding is an important method utilized in modern molecular breeding. Recent studies have determined that a large number of molecular markers appear to explain the impact of “lost heritability” on human height. Therefore, it is necessary to locate molecular marker sites in poultry and investigate the possible molecular mechanisms governing their effects. In this study, we found a 104-bp insertion/deletion polymorphism in the 5′UTR of the ADPRHL1 gene through resequencing. In cross-designed F2 resource groups, the indel was significantly associated with weight at 0, 2, 4, 6, and 10 weeks and a number of other traits [carcass weight (CW), semi-evisceration weight (SEW), evisceration weight (EW), claw weight (CLW), wings weight (DWW), gizzard weight (GW), pancreas weight (PW), chest muscle weight (CMW), leg weight (LW), leg muscle weight (LMW), shedding Weight (SW), liver rate (LR), and leg muscle rate (LMR)] (P &lt; 0.05). In brief, the insertion-insertion (II) genotype was significantly associated with the greatest growth traits and meat quality traits, whereas the values associated with the insertion-deletion (ID) genotype were the lowest in the F2 reciprocal cross chickens. The mutation sites were genotyped in 4,526 individuals from 12 different chicken breeds and cross-designed F2 resource groups. The II genotype is the most important genotype in commercial broilers, and the I allele frequency observed in these breeds is relatively high. Deletion mutations tend to be fixed in commercial broilers. However, there is still considerable great potential for breeding in dual-purpose chickens and commercial laying hens. A luciferase reporter assay showed that the II genotype of the ADPRHL1 gene possessed 2.49-fold higher promoter activity than the DD genotype (P &lt; 0.05). We hypothesized that this indel might affect the transcriptional activity of ADPRHL1, thereby affecting the growth traits of chickens. These findings may help to elucidate the function of the ADPRHL1 gene and facilitate enhanced reproduction in the chicken industry.

Mutability of mononucleotide repeats, not oxidative stress, explains the discrepancy between laboratory-accumulated mutations and the natural allele-frequency spectrum in C. elegans

Important clues about natural selection can be gleaned from discrepancies between the properties of segregating genetic variants and of mutations accumulated experimentally under minimal selection, provided the mutational process is the same in the laboratory as in nature. The base-substitution spectrum differs between C. elegans laboratory mutation accumulation (MA) experiments and the standing site-frequency spectrum, which has been argued to be in part owing to increased oxidative stress in the laboratory environment. Using genome sequence data from C. elegans MA lines carrying a mutation (mev-1) that increases the cellular titer of reactive oxygen species (ROS), leading to increased oxidative stress, we find the base-substitution spectrum is similar between mev-1, its wild-type progenitor (N2), and another set of MA lines derived from a different wild strain (PB306). Conversely, the rate of short insertions is greater in mev-1, consistent with studies in other organisms in which environmental stress increased the rate of insertion–deletion mutations. Further, the mutational properties of mononucleotide repeats in all strains are different from those of nonmononucleotide sequence, both for indels and base-substitutions, and whereas the nonmononucleotide spectra are fairly similar between MA lines and wild isolates, the mononucleotide spectra are very different, with a greater frequency of A:T → T:A transversions and an increased proportion of ±1-bp indels. The discrepancy in mutational spectra between laboratory MA experiments and natural variation is likely owing to a consistent (but unknown) effect of the laboratory environment that manifests itself via different modes of mutability and/or repair at mononucleotide loci.

Molecular Dynamics of SARS-CoV-2 Nsp1 Structural Changes Associated with Variant Mutations

SARS-CoV-2 non-structural protein 1 (Nsp1) is a virulence factor that inhibits the translation of host mRNAs and interact with viral RNA. Despite the relevance of Nsp1, few studies have been conducted to understand the effect of mutations on Nsp1 structure and function. Here, we provide a molecular dynamics simulation of SARS-CoV-2 Nsp1, wild type and variants. We found that SARS-CoV-2 Nsp1 has a more Rg value than SARS-CoV-1 Nsp1, with indicate an effect on the folding protein. This result suggest that SARS-CoV-2 Nsp1 can more easily approach the active site of the ribosome compared to SARS-CoV-1 Nsp1. In addition, we found that the C-terminal of the SARS-CoV-2 Nsp1, in particular residues 164 to 170, are more flexible than other regions of SARS-CoV-2 Nsp1 and SARS-CoV-1 Nsp1, confirming the role of this region in the interaction with the 40S subunit. Moreover, multiple deletion mutations have been found in the N/C-terminal of the SARS-CoV-2 Nsp1, which seems the effect of SARS-CoV-2 Nsp1 multiple deletions is greater than that of substitutions. Among all deletions, D156-158 and D80-90 may destabilize the protein structure and possibly increase the virulence of the SARS-CoV-2. Overall, our findings reinforce the importance of studying Nsp1 conformational changes in new variants and its effect on virulence of SARS-CoV-2.

Profiling a single-stranded DNA region within an rDNA segment that is a loading site for bacterial condensin

Bacterial condensin preferentially loads to single-stranded DNA (ssDNA) in vitro and loads onto rDNA in vivo to support proper chromosome compaction. Thus, the actively transcribing rDNA would provide the ssDNA region for the topological loading of bacterial condensin. We attempted to detect the ssDNA region in the rrnI gene in situ. Non-denaturing sodium bisulfite treatment catalyzed the conversion of cytosines to thymines via uracils (CT-conversion) at locally melted DNA of a bacterial genome. Using next-generation sequencing, we generated an average of 11,000 reads covering each cytosine on the PCR-amplified rDNA segment to obtain the actual CT-conversion rate. In principle, the CT-conversion rate is an accurate guide to detect the formation of the ssDNA segment. We expected that an increment of the CT-conversion rate would reflect a trend toward ssDNA accumulation at a given site within the rDNA. We detected multiple ssDNA segments throughout the rDNA. The deletion mutations of the rDNA that affect the bacterial-condensin loading hindered the ssDNA formation only at the 100-500 bp segment downstream of the promoter. These data support the idea that the ssDNA segment plays a crucial role as the bacterial condensin-loading site and suggest the mechanism of condensin loading onto rDNA.