Over the past thirty years, the identification of species-specific molecular markers has significantly advanced our understanding of genetic diversity in both plants and animals. Among these, short InDels have emerged as vital genomic features, contributing more to sequence divergence than single nucleotide polymorphisms do in closely related species. This study aimed to identify specific InDels for Bos taurus, Bubalus bubalis, Capra hircus, and Ovis aries via an in silico approach and validated them in 400 individuals (100 for each species).

We identified and characterized short, specific InDels in the sequences of the CSN1S1, CSN1S2, MSTN, and PRLR genes, which can be used for species identification of Capra hircus, Ovis aries, Bos taurus, and Bubalus bubalis, respectively. We developed a Tetraplex Specific PCR assay to enable efficient discrimination among these species.

This study highlights the utility of InDels as biallelic, codominant markers that are cost-effective and easy to analyse, providing valuable tools for genetic diversity analysis and species identification.

The effect of pressurized carbon dioxide microbubbles (CO2MB) on DNA in Saccharomyces pastorianus cells was investigated. The DNA concentration extracted from S. pastorianus cells treated with CO2MB by phenol extraction increased and decreased with the heating coil at 35 °C and at 40 and 45 °C, respectively. However, the density of DNA bands on agarose gel electrophoresis increased when using the heating coil at 40 and 45 °C. In addition, DNA fragmentation measured by quantitative PCR decreased with increasing temperature and prolonged exposure time in a heating coil. Furthermore, the generation of apurinic/apyrimidinic sites, an indicator of oxidative DNA damage, increased with increasing temperature or prolonged exposure time in a heating coil. In the phylogenetic tree constructed based on single-nucleotide polymorphisms obtained through DNA sequence analysis, S. pastorianus cells treated with CO2MB exhibited greater divergence from untreated with increasing temperature or prolonged exposure time in a heating coil. Similarly, insertions/deletions, showed genetic variations, of S. pastorianus cells treated with CO2MB increased with increasing temperature or prolonged exposure time in the heating coil. Therefore, these findings indicated that CO2MB induced partial oxidative damage and mutations of DNA in S. pastorianus cells, and the phenomena were suggested to be due to the oxidative stress generated within the cells.

Antiviral proteins often evolve rapidly at virus-binding interfaces to defend against new viruses. We investigated whether antiviral adaptation via missense mutations might face limits, which insertion or deletion mutations (indels) could overcome. Using high-throughput saturation missense mutagenesis, we identify one such case of a nearly insurmountable evolutionary challenge: the human anti-retroviral protein TRIM5α requires more than five missense mutations in its specificity-determining v1 loop to restrict a divergent simian immunodeficiency virus (SIV). However, through a novel saturating indel scanning methodology, we find that duplicating just one amino acid in v1 enables human TRIM5α to potently restrict SIV in a single evolutionary step. Moreover, natural primate TRIM5α v1 loops have evolved indels that confer novel antiviral specificities. Thus, indels enable antiviral proteins to overcome viral challenges otherwise inaccessible by missense mutations. Our findings reveal the potential of often-overlooked indel mutations in driving protein innovation.

Amino acid insertions and deletions (indels) are an abundant class of genetic variants. However, compared to substitutions, the effects of indels on protein stability are not well understood. To better understand indels here we analyse new and existing large-scale deep indel mutagenesis (DIM) of structurally diverse proteins. The effects of indels on protein stability vary extensively among and within proteins and are not well predicted by existing computational methods. To address this shortcoming we present INDELi, a series of models that combine experimental or predicted substitution effects and secondary structure information to provide good prediction of the effects of indels on both protein stability and pathogenicity. Moreover, quantifying the effects of indels on protein-protein interactions suggests that insertions can be an important class of gain-of-function variants. Our results provide an overview of the impact of indels on proteins and a method to predict their effects genome-wide.

Cancer neoantigens are peptides that originate from alterations in the genome, transcriptome, or proteome. These peptides can elicit cancer-specific T-cell recognition, making them potential candidates for cancer vaccines. The rapid advancement of proteomics technology holds tremendous potential for identifying these neoantigens. Here, we provided an up-to-date survey about database-based search methods and de novo peptide sequencing approaches in proteomics, and we also compared these methods to recommend reliable analytical tools for neoantigen identification. Unlike previous surveys on mass spectrometry-based neoantigen discovery, this survey summarizes the key advancements in de novo peptide sequencing approaches that utilize artificial intelligence. From a comparative study on a dataset of the HepG2 cell line and nine mixed hepatocellular carcinoma proteomics samples, we demonstrated the potential of proteomics for the identification of cancer neoantigens and conducted comparisons of the existing methods to illustrate their limits. Understanding these limits, we suggested a novel workflow for neoantigen discovery as perspectives.

Polymorphic short insertions and deletions (INDELs ≤ 50 bp) are abundant, although less common than single nucleotide polymorphisms (SNPs). Evidence from model organisms shows INDELs to be more strongly influenced by purifying selection than SNPs. Partly for this reason, INDELs are rarely used as markers for demographic processes or to detect divergent selection. Here, we compared INDELs and SNPs in the intertidal snail Littorina saxatilis, focussing on hybrid zones between ecotypes, in order to test the utility of INDELs in the detection of divergent selection. We computed INDEL and SNP site frequency spectra using capture sequencing data. We assessed the impact of divergent selection by analyzing allele frequency clines across habitat boundaries. We also examined the influence of GC-biased gene conversion because it may be confounded with signatures of selection. We show evidence that short INDELs are affected more by purifying selection than SNPs, but part of the observed site frequency spectra difference can be attributed to GC-biased gene conversion. We did not find a difference in the impact of divergent selection between short INDELs and SNPs. Short INDELs and SNPs were similarly distributed across the genome and so are likely to respond to indirect selection in the same way. A few regions likely affected by divergent selection were revealed by INDELs and not by SNPs. Short INDELs can be useful (additional) genetic markers helping to identify genomic regions important for adaptation and population divergence.

Engineering protein dynamics is a challenging and unsolved problem in protein design. Loop transplantation or loop grafting has been previously employed to transfer dynamic properties between proteins. We recently released a LoopGrafter Web server to execute the loop grafting task, employing eight computational tools and one database. The LoopGrafter method relies on the prediction of the local dynamic behavior of the elements to be transplanted and has successfully reconstructed previously engineered sequences. However, it was unclear whether catalytically competitive previously uncharacterized designs could be obtained by this method. Here, we address this question, showing how LoopGrafter generates viable loop-grafted chimeras of luciferases, how these chimeras encompass the activity of interest and unique kinetic properties, and how all this process is done fully automatically and agnostic of any previous knowledge. All constructed designs proved to be catalytically active, and the most active one improved the activity of the template enzyme by 4 orders of magnitude. The computational details and parameter optimization of the sequence pairing step of the LoopGrafter workflow are revealed. The optimized and experimentally validated loop grafting workflow available as a fully automated Web server represents a powerful approach for engineering catalytically efficient enzymes by modification of protein dynamics.

Gene variants resulting in insertions or deletions of amino acid residues (indels) have important consequences for evolution and are often linked to disease, yet, compared to missense variants, the effects of indels are poorly understood and predicted. We developed a sensitive protein folding sensor based on the complementation of uracil auxotrophy in yeast by circular permutated orotate phosphoribosyltransferase (CPOP). The sensor reports on the folding of disease-linked missense variants and de-novo-designed proteins. Applying the folding sensor to a saturated library of single-residue indels in human dihydrofolate reductase (DHFR) revealed that most regions that tolerate indels are confined to internal loops, the termini, and a central α helix. Several indels are temperature sensitive, and folding is rescued upon binding to methotrexate. Rosetta and AlphaFold2 predictions correlate with the observed effects, suggesting that most indels destabilize the native fold and that these computational tools are useful for the classification of indels observed in population sequencing.

Most of the risk factors associated with chronic and complex diseases, such as cancer, stem from exogenous and endogenous exposures experienced throughout an individual's life, collectively known as the exposome. These exposures can modify DNA, which can subsequently lead to the somatic mutations found in all normal and tumor tissues. Understanding the precise origins of specific somatic mutations has been challenging due to multitude of DNA adducts (i.e. the DNA adductome) and their diverse positions within the genome. Thus far, this limitation has prevented researchers from precisely linking exposures to DNA adducts and DNA adducts to subsequent mutational outcomes. Indeed, many common mutations observed in human cancers appear to originate from error-prone endogenous processes. Consequently, it remains unclear whether these mutations result from exposure-induced DNA adducts, or arise indirectly from endogenous processes or are a combination of both. In this review, we summarize approaches that aim to bridge our understanding of the mechanism by which exposure leads to DNA damage and then to mutation and highlight some of the remaining challenges and shortcomings to fully supporting this paradigm. We emphasize the need to integrate cellular DNA adductomics, long read-based mapping, single-molecule duplex sequencing of native DNA molecules and advanced computational analysis. This proposed holistic approach aims to unveil the causal connections between key DNA modifications and the mutational landscape, whether they originate from external exposures, internal processes or a combination of both, thereby addressing key questions in cancer biology.

Genome-wide association study (GWAS) have identified a large number of SNPs associated with milk production traits in dairy cattle. Behind SNPs, INDELs are the second most abundant genetic polymorphisms in the genome, which may exhibit an independent association with complex traits in humans and other species. However, there are no reports on GWASs of INDELs for milk production traits in dairy cattle. In this study, using imputed sequence data, we performed INDEL-based and SNP-based GWASs for milk production traits in a Holstein cattle population. We identified 58 unique significant INDELs for one or multiple traits. The majority of these INDELs are in considerable LD with nearby significant SNPs. However, through conditional association analysis, we identified nine INDELs which showed independent associations. Genomic annotations of these INDELs indicated some novel associated genes, i.e., TRNAG-CCC, EPPK1, PPM1K, PTDSS1, and mir-10163, which were not reported in previous SNP-based GWASs. Our findings suggest that INDEL-based GWASs could be valuable complement to SNP-based GWASs for milk production traits.

To investigate the potential regulatory role of gene insertion or deletion (in/del) polymorphism in the occurrence of acute T cell-mediated rejection (aTCMR) after kidney transplantation.

We retrospectively analyzed the 5-year follow-up data of 133 recipients who underwent renal transplantation at the First Affiliated Hospital of Nanjing Medical University between February 1, 2010, and December 1, 2015. With target sequencing based on next-generation sequencing (NGS), tagger in/dels selection involved calculating the Hardy-Weinberg equilibrium (HWE), Minor Allele Frequency (MAF), and the linkage disequilibrium (LD) blocks. Significant in/dels associated with aTCMR were identified by intersecting the results obtained through analysis of covariance (ANCOVA) of clinical cofounders and model analysis in Rstudio using the "SNPassoc" package. Additionally, logistic models were employed to assess the associations between genotypes and the aTCMR occurrence in 5 years after surgery.

NFATc1 rs55741427 insertion was identified to be significantly associated with the post-surgery aTCMR(OR = 2.66, P < 0.001). We constructed a conclusive model containing the occurrence of delayed graft function (DGF) and the insertion polymorphism of rs55741427, showing a favorable predictive ability (AUC = 0.766) for aTCMR after surgery. Based on the receiver operating characteristic (ROC) curve, all cases were stratified into aTCMR high-risk and low-risk groups. Kaplan-Meier curves for two groups revealed that the aTCMR high-risk group exhibited a more unfavorable graft survival outcome (P = 0.0048).

Insertion mutation of rs55741427 was found to be statistically correlated with the post-surgery aTCMR during 5 years of follow-up. Our model identified DGF and insertion of rs55741427 as two crucial aTCMR-related hazards, and aTCMR high-risk group showed a worse graft prognosis.

In India, 20 breeds of buffalo have been identified and registered, yet limited studies have been conducted to explore the performance potential of these breeds, especially in the Indian native breeds. This study is a maiden attempt to delineate the important variants and unique genes through exome sequencing for milk yield, milk composition, fertility, and adaptation traits in Indian local breeds of buffalo. In the present study, whole exome sequencing was performed on Chhattisgarhi (n = 3), Chilika (n = 4), Gojri (n = 3), and Murrah (n = 4) buffalo breeds and after stringent quality control, 4333, 6829, 4130, and 4854 InDels were revealed, respectively. Exome-wide FST along 100-kb sliding windows detected 27, 98, 38, and 35 outlier windows in Chhattisgarhi, Chilika, Gojri, and Murrah, respectively. The comparative exome analysis of InDels and subsequent gene ontology revealed unique breed specific genes for milk yield (CAMSAP3), milk composition (CLCN1, NUDT3), fertility (PTGER3) and adaptation (KCNA3, TH) traits. Study provides insight into mechanism of how these breeds have evolved under natural selection, the impact of these events on their respective genomes, and their importance in maintaining purity of these breeds for the traits under study. Additionally, this result will underwrite to the genetic acquaintance of these breeds for breeding application, and in understanding of evolution of these Indian local breeds.

Genetic mutations, both favorable and unfavorable, are the raw material for improvement in livestock populations. The random inheritance of these mutations is essential for generating progenies with genetic potential greater than their parents. These mutations can act either in a simple manner, such that a single alteration disrupts phenotype, or in a complex manner where hundreds or thousands of mutations of small effect create a continuous distribution of phenotypes. Selection tools leverage phenotypic records, pedigrees, and genomics to estimate the genetic potential of individual animals. This more accurate accounting of genetic potential has generated enormous gains in livestock populations.

The most recent version of ClinVar was utilized to filter variants of the MC4R gene based on location, condition, and clinical significance with the goal of obtaining benign and disease-associated variants of the MC4R gene. MC4R gene variants can lead to dysregulation of energy expenditure and appetite control, which prompted this study to delineate the distinctive features of MC4R gene variants submitted to the ClinVar repository regarding their association with obesity and related phenotypes.

A thorough search was conducted in the ClinVar repository for clinically significant MC4R variants through the utilization of the gene name MC4R[gene] and MeSH terms "MC4R[gene]" and "single gene"[properties]" in the search box. Leading to the identification of clinically significant genetic variants associated with obesity.

Utilizing the ClinVar clinical significance ranking system, the MC4R variants were categorized into six groups based on ClinVar/ClinGen's ranking system: pathogenic (P), likely pathogenic (LP), variant of uncertain significance (VUS), benign (B), likely benign (LB), and conflicting classifications (CC). A total of 103 pathogenic variants were observed. These variants have different clinical significance that are associated with monogenic obesity, monogenic diabetes, and body mass index quantitative traits. It was observed that over 80% of the mutations were single nucleotide variants, with nearly half being missense mutations spread throughout the topological and transmembrane domains. Furthermore, TM7 had the highest number of single nucleotide missense mutations.

Further analysis of the relationships between monogenic obesity and diabetes requires additional investigation to discover the underlying causes of these conditions. The study findings imply that mutations in MC4R's topological and transmembrane regions may significantly influence receptor activation and signaling. As more MC4R variants are discovered and their correlation with obesity is established, there is potential to definitively establish a strong connection between MC4R pathogenic variants and the development of obesity.

The search for genetic variants that act as causative factors in human diseases by disrupting the normal splicing process has primarily focused on single nucleotide variants (SNVs). It is worth noting that insertions or deletions (indels) have also been sporadically reported as causative disease variants through their potential impact on the splicing process. In this study, to perform identification of indels inducing exon extension/shrinkage events, we used individual-specific genomes and RNA sequencing (RNA-seq) data pertaining to the corresponding individuals and identified 12 exon extension/shrinkage events that were potentially induced by indels that disrupted authentic splice sites or created novel splice sites in 235 normal individuals. By evaluating the impact of these abnormal splicing events on the resulting transcripts, we found that five events led to the generation of premature termination codons (PTCs), including those occurring within genes associated with genetic disorders. Our analysis revealed that the potential functions of indels have been underexamined, and it is worth considering the possibility that indels may affect splice site usage, using RNA-seq data to discover novel potentially disease-associated mutations.

Insertions and deletions constitute the second most important source of natural genomic variation. Insertions and deletions make up to 25% of genomic variants in humans and are involved in complex evolutionary processes including genomic rearrangements, adaptation, and speciation. Recent advances in long-read sequencing technologies allow detailed inference of insertions and deletion variation in species and populations. Yet, despite their importance, evolutionary studies have traditionally ignored or mishandled insertions and deletions due to a lack of comprehensive methodologies and statistical models of insertions and deletion dynamics. Here, we discuss methods for describing insertions and deletion variation and modeling insertions and deletions over evolutionary time. We provide practical advice for tackling insertions and deletions in genomic sequences and illustrate our discussion with examples of insertions and deletion-induced effects in human and other natural populations and their contribution to evolutionary processes. We outline promising directions for future developments in statistical methodologies that would allow researchers to analyze insertions and deletion variation and their effects in large genomic data sets and to incorporate insertions and deletions in evolutionary inference.

Therapeutic cancer vaccines have been considered in recent decades as important immunotherapeutic strategies capable of leading to tumor regression. In the development of these vaccines, the identification of neoepitopes plays a critical role, and different computational methods have been proposed and employed to direct and accelerate this process. In this context, this review identified and systematically analyzed the most recent studies published in the literature on the computational prediction of epitopes for the development of therapeutic vaccines, outlining critical steps, along with the associated program's strengths and limitations. A scoping review was conducted following the PRISMA extension (PRISMA-ScR). Searches were performed in databases (Scopus, PubMed, Web of Science, Science Direct) using the keywords: neoepitope, epitope, vaccine, prediction, algorithm, cancer, and tumor. Forty-nine articles published from 2012 to 2024 were synthesized and analyzed. Most of the identified studies focus on the prediction of epitopes with an affinity for MHC I molecules in solid tumors, such as lung carcinoma. Predicting epitopes with class II MHC affinity has been relatively underexplored. Besides neoepitope prediction from high-throughput sequencing data, additional steps were identified, such as the prioritization of neoepitopes and validation. Mutect2 is the most used tool for variant calling, while NetMHCpan is favored for neoepitope prediction. Artificial/convolutional neural networks are the preferred methods for neoepitope prediction. For prioritizing immunogenic epitopes, the random forest algorithm is the most used for classification. The performance values related to the computational models for the prediction and prioritization of neoepitopes are high; however, a large part of the studies still use microbiome databases for training. The in vitro/in vivo validations of the predicted neoepitopes were verified in 55% of the analyzed studies. Clinical trials that led to successful tumor remission were identified, highlighting that this immunotherapeutic approach can benefit these patients. Integrating high-throughput sequencing, sophisticated bioinformatics tools, and rigorous validation methods through in vitro/in vivo assays as well as clinical trials, the tumor neoepitope-based vaccine approach holds promise for developing personalized therapeutic vaccines that target specific tumor cancers.

Short-Patch Double Illegitimate Recombination (SPDIR) has been recently identified as a rare mutation mechanism. During SPDIR, ectopic DNA single-strands anneal with genomic DNA at microhomologies and get integrated during DNA replication, presumably acting as primers for Okazaki fragments. The resulting microindel mutations are highly variable in size and sequence. In the soil bacterium Acinetobacter baylyi, SPDIR is tightly controlled by genome maintenance functions including RecA. It is thought that RecA scavenges DNA single-strands and renders them unable to anneal. To further elucidate the role of RecA in this process, we investigate the roles of the upstream functions DprA, RecFOR, and RecBCD, all of which load DNA single-strands with RecA. Here we show that all three functions suppress SPDIR mutations in the wildtype to levels below the detection limit. While SPDIR mutations are slightly elevated in the absence of DprA, they are strongly increased in the absence of both DprA and RecA. This SPDIR-avoiding function of DprA is not related to its role in natural transformation. These results suggest a function for DprA in combination with RecA to avoid potentially harmful microindel mutations, and offer an explanation for the ubiquity of dprA in the genomes of naturally non-transformable bacteria.

CD59 deficiency due to rare germline variants in the CD59 gene causes disabilities, ischemic strokes, neuropathy, and hemolysis. CD59 deficiency due to common somatic variants in the PIG-A gene in hematopoietic stem cells causes paroxysmal nocturnal hemoglobinuria. The ISBT database lists one nonsense and three missense germline variants that are associated with the CD59-null phenotype. To analyze the genetic diversity of the CD59 gene, we determined long-range CD59 haplotypes among individuals from different ethnicities.

We determined a 22.7 kb genomic fragment of the CD59 gene in 113 individuals using next-generation sequencing (NGS), which covered the whole NM_203330.2 mRNA transcript of 7796 base pairs. Samples came from an FDA reference repository and our Ethiopia study cohorts. The raw genotype data were computationally phased into individual haplotype sequences.

Nucleotide sequencing of the CD59 gene of 226 chromosomes identified 216 positions with single nucleotide variants. Only three haplotypes were observed in homozygous form, which allowed us to assign them unambiguously as experimentally verified CD59 haplotypes. They were also the most frequent haplotypes among both cohorts. An additional 140 haplotypes were imputed computationally.

We provided a large set of haplotypes and proposed three verified long-range CD59 reference sequences, based on a population approach, using a generalizable rationale for our choice. Correct long-range haplotypes are useful as template sequences for allele calling in high-throughput NGS and precision medicine approaches, thus enhancing the reliability of clinical diagnostics. Long-range haplotypes can also be used to evaluate the influence of genetic variation on the risk of transfusion reactions or diseases.

The use of genetic markers, specifically Short Tandem Repeats (STRs), has been a valuable tool for identifying persons of interest. However, the ability to analyze additional markers including Single Nucleotide Polymorphisms (SNPs) and Insertion/Deletion (INDELs) polymorphisms allows laboratories to explore other investigative leads. INDELs were chosen in this study because large panels can be differentiated by size, allowing them to be genotyped by capillary electrophoresis. Moreover, these markers do not produce stutter and are smaller in size than STRs, facilitating the recovery of genetic information from degraded samples. The INDEL Ancestry Informative Markers (AIMs) in this study were selected from the 1000 Genomes Project based on a fixation index (FST) greater than 0.50, high allele frequency divergence, and genetic distance. A total of 25 INDEL-AIMs were optimized and validated according to SWGDAM guidelines in a five-dye multiplex. To validate the panel, genotyping was performed on 155 unrelated individuals from four ancestral groups (Caucasian, African, Hispanic, and East Asian). Bayesian clustering and principal component analysis (PCA) were performed revealing clear separation among three groups, with some observed overlap within the Hispanic group. Additionally, the PCA results were compared against a training set of 793 samples from the 1000 Genomes Project, demonstrating consistent results. Validation studies showed the assay to be reproducible, tolerant to common inhibitors, robust with challenging casework type samples, and sensitive down to 125 pg. In conclusion, our results demonstrated the robustness and effectiveness of a 25 loci INDEL system for ancestry inference of four ancestries commonly found in the United States.

Antiviral proteins often evolve rapidly at virus-binding interfaces to defend against new viruses. We investigated whether antiviral adaptation via missense mutations might face limits, which insertion or deletion mutations (indels) could overcome. We report one such case of a nearly insurmountable evolutionary challenge: the human anti-retroviral protein TRIM5α requires more than five missense mutations in its specificity-determining v1 loop to restrict a divergent simian immunodeficiency virus (SIV). However, duplicating just one amino acid in v1 enables human TRIM5α to potently restrict SIV in a single evolutionary step. Moreover, natural primate TRIM5α v1 loops have evolved indels that confer novel antiviral specificities. Thus, indels enable antiviral proteins to overcome viral challenges inaccessible by missense mutations, revealing the potential of these often-overlooked mutations in driving protein innovation.

DNA analysis of forensic case samples relies on short tandem repeats (STRs), a key component of the combined DNA index system (CODIS) used to identify individuals. However, limitations arise when dealing with challenging samples, prompting the exploration of alternative markers such as single nucleotide polymorphisms (SNPs) and insertion/deletion (INDELs) polymorphisms. Unlike SNPs, INDELs can be differentiated easily by size, making them compatible with electrophoresis methods. It is possible to design small INDEL amplicons (<200 bp) to enhance recovery from degraded samples. To this end, a set of INDEL Human Identification Markers (HID) was curated from the 1000 Genomes Project, employing criteria including a fixation index (FST) ≤ 0.06, minor allele frequency (MAF) >0.2, and high allele frequency divergence. A panel of 33 INDEL-HIDs was optimized and validated following the Scientific Working Group on DNA Analysis Methods (SWGDAM) guidelines, utilizing a five-dye multiplex electrophoresis system. A small sample set (n = 79 unrelated individuals) was genotyped to assess the assay's performance. The validation studies exhibited reproducibility, inhibition tolerance, ability to detect a two-person mixture from a 4:1 to 1:6 ratio, robustness with challenging samples, and sensitivity down to 125 pg of DNA. In summary, the 33-loci INDEL-HID panel exhibited robust recovery with low-template and degraded samples and proved effective for individualization within a small sample set.

Despite high abundance of small indels in human genomes, their precise roles and underlying mechanisms of mutagenesis in Mendelian disorders require further investigation.

To profile the distribution, functional implications, and mechanisms of small indels in the androgen receptor (AR) gene in individuals with androgen insensitivity syndrome (AIS).

We conducted a systematic review of previously reported indels within the coding region of the AR gene, including 3 novel indels. Distribution throughout the AR coding region was examined and compared with genomic population data. Additionally, we assessed their impact on the AIS phenotype and investigated potential mechanisms driving their occurrence.

A total of 82 indels in AIS were included. Notably, all frameshift indels exhibited complete AIS. The distribution of indels across the AR gene showed a predominance in the N-terminal domain, most leading to frameshift mutations. Small deletions accounted for 59.7%. Most indels occurred in nonrepetitive sequences, with 15.8% situated within triplet regions. Gene burden analysis demonstrated significant enrichment of frameshift indels in AIS compared with controls (P < .00001), and deletions were overrepresented in AIS (P < .00001).

Our findings underscore a robust genotype-phenotype relationship regarding small indels in the AR gene in AIS, with a vast majority presenting complete AIS. Triplet regions and homopolymeric runs emerged as prone loci for small indels within the AR. Most were frameshift indels, with polymerase slippage potentially explaining half of AR indel occurrences. Complex frameshift indels exhibited association with palindromic runs. These discoveries advance understanding of the genetic basis of AIS and shed light on potential mechanisms underlying pathogenic small indel events.

Microhaplotypes (MHs), small sets of linked single nucleotide polymorphisms (SNPs), are becoming a valuable tool for paternity testing, personal identification and other different forensic purposes due to their advantages of both short tandem repeats (STRs) and SNPs. However, only a small part of MHs with small segments have been developed and reported so far. And the current population genetic data of MHs are still insufficient. MHs with small segments possess unique advantages in mixture deconvolution, degradation material identification, noninvasive prenatal paternity testing and even medical tumor diagnostic applications. In the present study, a set of 90 autosomal MHs whose PCR amplicon lengths are from 90-150 bp, of which 58 MHs are less than or equal to 100 bp are selected, and assembled into an amplification multiplex system optimized for Ion S5™ System for forensic application. Genetic diversity study of 90 MHs in the populations from different intercontinental regions shows that the polymorphism information content (PIC) values of 83 MHs are greater than 0.4 in populations from East Asia (EAS), and the average PIC value of 90 MHs is greater than 0.5. A total of EAS populations shows the highest cumulative match probability (CMP) and cumulative probability of exclusion (CPE) values in five intercontinental populations. The CMP and CPE values of 90 MHs in EAS are 1.1688 × 10-54 and 0.999999999998954. The informativeness for assignment (In) values of the 90 MHs are calculated based on data from five intercontinental populations, and the In values of 20 MHs have greater than 0.1, indicating that the 20 MHs are high effectiveness in distinguishing different intercontinental populations, which can be used as candidate ancestry informative markers. Further, we have studied the polymorphisms of the 90 MHs based on 224 unrelated individuals of Henan Han population, China, and obtained the frequency data of the 90 MHs. In the Henan Han population, the effective number of alleles (Ae) of the 90 MHs ranges from 1.7649 (MH45) to 3.9792 (MH50), and the Ae values of 10 MHs reach to 3.0; the Ae values of 80 MHs are greater than 2, and the average Ae value for these MHs is 2.422. The average expected heterozygosity, observed heterozygosity, PIC, matching probability, discrimination power and probability of exclusion values of 90 MHs in the Henan Han population are 0.5788, 0.5851, 0.5039, 0.2608, 0.7392 and 0.2806, respectively. The CMP value of 90 MHs in the study population is less than 10-54, and their CPE value reaches 0.999999999999999923. Moreover, the results of the depth of coverage, allele coverage ratio and noise level indicate that the 90 MH amplification system has well sequencing performance, and the sequencing results are reliable. The results indicate the 90 MHs show higher polymorphisms in the study population. The present panel can be well used in paternity testing and individual identification in the study population and even the populations from EAS.

Within anatomical willed body programs and skeletal collections, whole bodies and their disassociated limbs and organs are identified and tracked. However, if these tracking mechanisms fail, DNA recovered from the formalin-fixed tissues/organs could provide an additional layer of quality assurance. Embalming fluids preserve biological tissues; however, they also damage, fragment, and cross-link DNA and protein molecules. This project investigated the success of STR-typing from various soft tissue and bone samples that were fixed with embalming solutions with a range of formaldehyde concentrations. Formalin-fixed samples dissected from five cadavers, including skin, muscle, bone, heart, and kidney were used in Phase 1 of this study. In Phase 2, an additional 57 tissue samples from various embalmed organs and body parts were collected to demonstrate long-term fixation and direct applicability within a body donor program. DNA was extracted from the samples using the QIAamp® FFPE Tissue Kit (QIAGEN), quantified with the Investigator® QuantiPlex® Pro RGQ qPCR Kit (QIAGEN), and amplified using the Investigator® 24plex and 26plex QS Kits and the Investigator® DIPplex Kit (QIAGEN). The results show the DNA was severely damaged, degraded, and often in low amounts (after one year post-embalming). Sampling from skin and muscle tissues embalmed with ~2.5%-5% formaldehyde solutions appears to be the best strategy for identification, while also maintaining the preservation of the tissues. The results of this project can provide informative data when determining which genotyping strategy may be best suited for the identification, re-association, and establishment of a database for the provenance of formalin-fixed human remains.

Despite clinically demonstrated accuracy in next generation sequencing (NGS) data, many clinical laboratories continue to confirm variants with Sanger sequencing, which increases cost of testing and turnaround time. Several studies have assessed the accuracy of NGS in detecting single nucleotide variants; however, less has been reported about insertion, deletion, and deletion-insertion variants (indels).

We performed a retrospective analysis from 2015-2022 of indel results from a subset of NGS targeted gene panel tests offered through the Mayo Clinic Genomics Laboratories. We compared results from NGS and Sanger sequencing of indels observed in clinical runs and during the intra-assay validation of the tests.

Results demonstrated 100% concordance between NGS and Sanger sequencing for over 490 indels (217 unique), ranging in size from 1 to 68 basepairs (bp). The majority of indels were deletions (77%) and 1 to 5 bp in length (90%). Variant frequencies ranged from 11.4% to 67.4% and 85.1% to 100% for heterozygous and homozygous variants, respectively, with a median depth of coverage of 2562×. A subset of indels (7%) were located in complex regions of the genome, and these were accurately detected by NGS. We also demonstrated 100% reproducibility of indel detection (n = 179) during intra-assay validation.

Together this data demonstrates that reportable indel variants up to 68 bp can be accurately assessed using NGS, even when they occur in complex regions. Depending on the complexity of the region or variant, Sanger sequence confirmation of indels is usually not necessary if the variants meet appropriate coverage and allele frequency thresholds.

Renal cell carcinoma (RCC) comprises a group of malignancies arising from the kidney with unique tumour-specific antigen (TSA) signatures that can trigger cytotoxic immunity. Two classes of TSAs are now considered potential drivers of immunogenicity in RCC: small-scale insertions and deletions (INDELs) that result in coding frameshift mutations, and activation of human endogenous retroviruses. The presence of neoantigen-specific T cells is a hallmark of solid tumours with a high mutagenic burden, which typically have abundant TSAs owing to non-synonymous single nucleotide variations within the genome. However, RCC exhibits high cytotoxic T cell reactivity despite only having an intermediate non-synonymous single nucleotide variation mutational burden. Instead, RCC tumours have a high pan-cancer proportion of INDEL frameshift mutations, and coding frameshift INDELs are associated with high immunogenicity. Moreover, cytotoxic T cells in RCC subtypes seem to recognize tumour-specific endogenous retrovirus epitopes, whose presence is associated with clinical responses to immune checkpoint blockade therapy. Here, we review the distinct molecular landscapes in RCC that promote immunogenic responses, discuss clinical opportunities for discovery of biomarkers that can inform therapeutic immune checkpoint blockade strategies, and identify gaps in knowledge for future investigations.

The insertion or deletion (indel) of amino acids has a variety of effects on protein function, ranging from disease-forming changes to gaining new functions. Despite their importance, indels have not been systematically characterized towards protein engineering or modification goals. In the present work, we focus on deletions composed of multiple contiguous amino acids (mAA-dels) and their effects on the protein (mutant) folding ability. Our analysis reveals that the mutant retains the native fold when the mAA-del obeys well-defined structural dynamics properties: localization in intrinsically flexible regions, showing low resistance to mechanical stress, and separation from allosteric signaling paths. Motivated by the possibility of distinguishing the features that underlie the adaptability of proteins to mAA-dels, and by the rapid evaluation of these features using elastic network models, we developed a positive-unlabeled learning-based classifier that can be adopted for protein design purposes. Trained on a consolidated set of features, including those reflecting the intrinsic dynamics of the regions where the mAA-dels occur, the new classifier yields a high recall of 84.3% for identifying mAA-dels that are stably tolerated by the protein. The comparative examination of the relative contribution of different features to the prediction reveals the dominant role of structural dynamics in enabling the adaptation of the mutant to mAA-del without disrupting the native fold.

Yunnan is one of the main residences of the Zhuang group which is one of the 55 ethnic minorities in China. At present, there are relatively few researches on population genetics and forensic science of the Yunnan Zhuang group. Therefore, this study used a self-constructed panel containing 41 multi-InDel markers to analyze the genetic polymorphisms of 173 individuals from Yunnan Zhuang group. The results indicated that these 41 multi-InDels in Yunnan Zhuang group were highly polymorphic markers expect for three markers. The cumulative match probability and combined exclusion probability values of the 40 multi-InDels (MI38 marker was excluded) were 8.0671E-26 and 0.9999995959, respectively. In addition, population genetic analyses were performed on genotyping data of 41 multi-InDel markers among the Yunnan Zhuang and 26 reference populations, revealing that the Yunnan Zhuang group was genetically close to the five populations in East Asia. According to the STRUCTURE analysis, the Yunnan Zhuang group presented similar ancestral compositions to the five populations from East Asia, and when the K value was three, the five intercontinental populations showed their different genetic structures. In conclusion, the 41 multi-InDel markers could be used as an effective tool for individual identification and paternity testing of the Zhuang group in Yunnan province, as well as for their ancestry information inference studies.

The use of in silico pathogenicity predictions as evidence when interpreting genetic variants is widely accepted as part of standard variant classification guidelines. Although numerous algorithms have been developed and evaluated for classifying missense variants, in-frame insertions/deletions (indels) have been much less well studied.

We created a dataset of 3964 small (< 100 bp) indels predicted to result in in-frame amino acid insertions or deletions using data from gnomAD v3.1 (minor allele frequency of 1-5%), ClinVar and the Deciphering Developmental Disorders (DDD) study. We used this dataset to evaluate the performance of nine pathogenicity predictor tools: CADD, CAPICE, FATHMM-indel, MutPred-Indel, MutationTaster2021, PROVEAN, SIFT-indel, VEST-indel and VVP.

Our dataset consisted of 2224 benign/likely benign and 1740 pathogenic/likely pathogenic variants from gnomAD (n = 809), ClinVar (n = 2882) and, DDD (n = 273). We were able to generate scores across all tools for 91% of the variants, with areas under the ROC curve (AUC) of 0.81-0.96 based on the published recommended thresholds. To avoid biases caused by inclusion of our dataset in the tools' training data, we also evaluated just DDD variants not present in either gnomAD or ClinVar (70 pathogenic and 81 benign). Using this subset, the AUC of all tools decreased substantially to 0.64-0.87. Several of the tools performed similarly however, VEST-indel had the highest AUCs of 0.93 (full dataset) and 0.87 (DDD subset).

Algorithms designed for predicting the pathogenicity of in-frame indels perform well enough to aid clinical variant classification in a similar manner to missense prediction tools.

Structural variation (SV), which ranges from 50 bp to [Formula: see text] 3 Mb in size, is an important type of genetic variations. Deletion is a type of SV in which a part of a chromosome or a sequence of DNA is lost during DNA replication. Three types of signals, including discordant read-pairs, reads depth and split reads, are commonly used for SV detection from high-throughput sequence data. Many tools have been developed for detecting SVs by using one or multiple of these signals.

In this paper, we develop a new method called EigenDel for detecting the germline submicroscopic genomic deletions. EigenDel first takes advantage of discordant read-pairs and clipped reads to get initial deletion candidates, and then it clusters similar candidates by using unsupervised learning methods. After that, EigenDel uses a carefully designed approach for calling true deletions from each cluster. We conduct various experiments to evaluate the performance of EigenDel on low coverage sequence data.

Our results show that EigenDel outperforms other major methods in terms of improving capability of balancing accuracy and sensitivity as well as reducing bias. EigenDel can be downloaded from https://github.com/lxwgcool/EigenDel .

Over the years, protein engineers have studied nature and borrowed its tricks to accelerate protein evolution in the test tube. While there have been considerable advances, our ability to generate new proteins in the laboratory is seemingly limited. One explanation for these shortcomings may be that insertions and deletions (indels), which frequently arise in nature, are largely overlooked during protein engineering campaigns. The profound effect of indels on protein structures, by way of drastic backbone alterations, could be perceived as "saltation" events that bring about significant phenotypic changes in a single mutational step. Should we leverage these effects to accelerate protein engineering and gain access to unexplored regions of adaptive landscapes? In this Perspective, we describe the role played by indels in the functional diversification of proteins in nature and discuss their untapped potential for protein engineering, despite their often-destabilizing nature. We hope to spark a renewed interest in indels, emphasizing that their wider study and use may prove insightful and shape the future of protein engineering by unlocking unique functional changes that substitutions alone could never achieve.

BMPR1B (Bone morphogenetic protein receptor type-1B) is a receptor in the bone morphogenetic protein (BMP) family and has been identified as a candidate gene for reproductive traits in pigs. Our previous study in Taihu pigs found a specific estrogen response element (ERE) in the first intron of the BMPR1B gene that is associated with the number born alive trait. However, little is known about the mechanism by which the ERE regulates the expression of BMPR1B in the endometrium.

Here, a 15-bp InDel (insertion/deletion) (AGCCAGAAAGGAGGA) was identified as a unique variation in Taihu pigs, and was shown to be responsible for the binding of the type I receptor of estrogen (ESR1) to the ERE using dual-luciferase assays. Four BMPR1B transcripts (T1, T2, T3, and T4) were identified by 5' RACE in endometrial tissue. Expression of T3 and T4 in the endometrium of Meishan pigs was significantly higher than in Duroc pigs during pregnancy. Luciferase assays showed that three distinct BMPR1B promoters may drive expression of T1, T3, and T4. Interestingly, ERE-mediated enhancement of T4 promoter activity significantly increased expression of Transcript T4 in the endometrium of Taihu pigs (P < 0.05). In contrast, the ERE inhibited activity of the T3 promoter and decreased expression of the T3 transcript in the Duroc background (P < 0.05). In summary, we identified a 15-bp InDel in the Taihu ERE that can be used as a molecular marker for the number born alive trait, characterized the 5' untranslated regions (UTRs) of BMPR1B transcripts in the endometrium, and determined how the transcripts are processed by alternative splicing events.

Our results provide a foundation for understanding the transcriptional regulation of BMPR1B and its contributions to the unique breeding prolificacy characteristics of Taihu pigs.

Small in-frame insertion-deletion (indel) variants are a common form of genomic variation whose impact on rare disease phenotypes has been understudied. The prediction of the pathogenicity of such variants remains challenging. X-linked incomplete congenital stationary night blindness type 2 (CSNB2) is a nonprogressive, inherited retinal disorder caused by variants in CACNA1F, encoding the Ca_v1.4α1 channel protein. Here, structural analysis was used through homology modeling to interpret 10 disease-correlated and 10 putatively benign CACNA1F in-frame indel variants. CSNB2-correlated changes were found to be more highly conserved compared with putative benign variants. Notably, all 10 disease-correlated variants but none of the benign changes were within modeled regions of the protein. Structural analysis revealed that disease-correlated variants are predicted to destabilize the structure and function of the Ca_v1.4α1 channel protein. Overall, the use of structural information to interpret the consequences of in-frame indel variants provides an important adjunct that can improve the diagnosis for individuals with CSNB2.

DNA sequencing is the physical/biochemical process of identifying the location of the four bases (Adenine, Guanine, Cytosine, Thymine) in a DNA strand. As semiconductor technology revolutionized computing, modern DNA sequencing technology (termed Next Generation Sequencing, NGS) revolutionized genomic research. As a result, modern NGS platforms can sequence hundreds of millions of short DNA fragments in parallel. The sequenced DNA fragments, representing the output of NGS platforms, are termed reads. Besides genomic variations, NGS imperfections induce noise in reads. Mapping each read to (the most similar portion of) a reference genome of the same species, i.e., read mapping, is a common critical first step in a diverse set of emerging bioinformatics applications. Mapping represents a search-heavy memory-intensive similarity matching problem, therefore, can greatly benefit from near-memory processing. Intuition suggests using fast associative search enabled by Ternary Content Addressable Memory (TCAM) by construction. However, the excessive energy consumption and lack of support for similarity matching (under NGS and genomic variation induced noise) renders direct application of TCAM infeasible, irrespective of volatility, where only non-volatile TCAM can accommodate the large memory footprint in an area-efficient way. This paper introduces GeNVoM, a scalable, energy-efficient and high-throughput solution. Instead of optimizing an algorithm developed for general-purpose computers or GPUs, GeNVoM rethinks the algorithm and non-volatile TCAM-based accelerator design together from the ground up. Thereby GeNVoM can improve the throughput by up to 3.67×; the energy consumption, by up to 1.36×, when compared to an ASIC baseline, which represents one of the highest-throughput implementations known.

Estimating the effects of variants found in disease driver genes opens the door to personalized therapeutic opportunities. Clinical associations and laboratory experiments can only characterize a tiny fraction of all the available variants, leaving the majority as variants of unknown significance (VUS). In silico methods bridge this gap by providing instant estimates on a large scale, most often based on the numerous genetic differences between species. Despite concerns that these methods may lack reliability in individual subjects, their numerous practical applications over cohorts suggest they are already helpful and have a role to play in genome interpretation when used at the proper scale and context. In this review, we aim to gain insights into the training and validation of these variant effect predicting methods and illustrate representative types of experimental and clinical applications. Objective performance assessments using various datasets that are not yet published indicate the strengths and limitations of each method. These show that cautious use of in silico variant impact predictors is essential for addressing genome interpretation challenges.

Many short-read genome assemblies have been found to be incomplete and contain mis-assemblies. The Vertebrate Genomes Project has been producing new reference genome assemblies with an emphasis on being as complete and error-free as possible, which requires utilizing long reads, long-range scaffolding data, new assembly algorithms, and manual curation. A more thorough evaluation of the recent references relative to prior assemblies can provide a detailed overview of the types and magnitude of improvements.

Here we evaluate new vertebrate genome references relative to the previous assemblies for the same species and, in two cases, the same individuals, including a mammal (platypus), two birds (zebra finch, Anna's hummingbird), and a fish (climbing perch). We find that up to 11% of genomic sequence is entirely missing in the previous assemblies. In the Vertebrate Genomes Project zebra finch assembly, we identify eight new GC- and repeat-rich micro-chromosomes with high gene density. The impact of missing sequences is biased towards GC-rich 5'-proximal promoters and 5' exon regions of protein-coding genes and long non-coding RNAs. Between 26 and 60% of genes include structural or sequence errors that could lead to misunderstanding of their function when using the previous genome assemblies.

Our findings reveal novel regulatory landscapes and protein coding sequences that have been greatly underestimated in previous assemblies and are now present in the Vertebrate Genomes Project reference genomes.

Although the majority of RET alterations are single nucleotide variants (SNV), small deletions and/or insertions have been reported at variable prevalence. No information about the efficacy of RET-specific inhibitors in patients harboring RET indels has been provided.

We present an update on the prevalence of RET indels in medullary thyroid cancer (MTC) and describe the efficacy of selpercatinib in patients with advanced MTC with RET indels.

The MTC tissues of 287 patients were analyzed using an Ion S5 targeted sequencing. The functional role of the reported indels have been evaluated by MutationTaster. Clinical and pathological data of MTC patients harboring a RET indel were collected and analyzed. Two patients with a RET indel were treated with selpercatinib.

Among 178 RET-positive cases, 147 (82.6%) harbored a SNV and 31 (17.4%) a RET in-frame indel. Nine indels were not previously reported and were found to be disease causing by MutationTaster. Patients harboring an indel were found to have an aggressive disease and 2 of them were treated with selpercatinib, experiencing a good response to the treatment.

These data show that RET indels are not infrequent and correlate with an aggressive disease. Two RET indel-positive patients showed a partial response to the treatment with a highly selective RET inhibitor; thus, these RET indels can be considered actionable mutations. In order to not miss these alterations, the analysis of the full gene is recommended.

Human DNA contains several variations, which can affect the structure and normal functioning of a protein. These variations could be single nucleotide polymorphisms (SNPs) or insertion-deletions (InDels). SNPs, as opposed to InDels, are more commonly present in DNA and may cause genetic disorders. In the current study, several bioinformatic tools were used to prioritize the pathogenic variants in the SLITRK1 gene. Out of all of the variants, 16 were commonly predicted to be pathogenic by these tools. All the variants had very low frequency, i.e., <0.0001 in the global population. The secondary structure of all filtered variants was predicted, but no structural change was observed at the site of variation in any variant. Protein stability analysis of these variants was then performed, which determined a decrease in protein stability of 10 of the variants. Amino acid conservation analysis revealed that all the amino acids were highly conserved, indicating their structural and functional importance. Protein 3D structure of wildtype SLITRK1 and all of its variants was predicted using I-TASSER, and the effect of variation on 3D structure of the protein was observed using the Missense3D tool, which presented the probable structural loss in three variants, i.e., Asn529Lys, Leu496Pro and Leu94Phe. The wildtype SLITRK1 protein and these three variants were independently docked with their close interactor protein PTPRD, and remarkable differences were observed in the docking sites of normal and variants, which will ultimately affect the functional activity of the SLITRK1 protein. Previous studies have shown that mutations in SLITRK1 are involved in Tourette syndrome. The present study may assist a molecular geneticist in interpreting the variant pathogenicity in research as well as diagnostic setup.

Discovery of the CRISPR-Cas (clustered regularly interspaced short palindromic repeat, CRISPR-associated) system a decade ago has opened new possibilities in the field of precision medicine. CRISPR-Cas was initially identified in bacteria and archaea to play a protective role against foreign genetic elements during viral infections. The application of this technique for the correction of different mutations found in the Duchenne muscular dystrophy (DMD) gene led to the development of several potential therapeutic approaches for DMD patients. The mutations responsible for Duchenne muscular dystrophy mainly include exon deletions (70% of patients) and point mutations (about 30% of patients). The CRISPR-Cas 9 technology is becoming increasingly precise and is acquiring diverse functions through novel innovations such as base editing and prime editing. However, questions remain about its translation to the clinic. Current research addressing off-target editing, efficient muscle-specific delivery, immune response to nucleases, and vector challenges may eventually lead to the clinical use of the CRISPR-Cas9 technology. In this review, we present recent CRISPR-Cas9 strategies to restore dystrophin expression in vitro and in animal models of DMD.