Genic Variants Research Articles

Mutation profiling through whole genome sequencing of electron beam-induced black gram (Vigna mungo L. Hepper) mutant

Purpose Black gram (Vigna mungo [L.] Hepper) is an important annual legume with great economic, nutritional and ecological significance. Novel variations through induced mutagenesis can accelerate narrow genetic base-impeded black gram improvement. This is a first study on characterization of genome-wide mutation spectrum induced by electron beam (EB). Materials and methods Black gram genotype ‘Pant U-31’ was irradiated with 400 Gy EB generated in a 10 MeV LINAC. A stable mutant PM-32 (M6) was re-sequenced by combining Illumina (BIOO Scientific, Inc., Austin, TX) and Nanopore Technologies (Oxford, UK). Variants were predicted in reference to the available whole genome scaffold level draft assembly of parent ‘Pant U-31’. Results Genome analysis predicted a total of 76,893 genes of which 58,517 were annotated. The identified variants totaling 728,161, largely comprised (91.56%) of single base substitutions (SBSs) with a transition (Ti) to transversion (Tv) ratio of 1.95. Of the indels constituting 8.44% of total induced variants, insertions accounted for 4.29%, with preponderance of multiple bases (53.63%) and 2–5 bp insertions as the major class (33.71%). Multiple-base deletions (2–5 bases) formed the bulk (31.14%) of the total deletions. The genic variants (2438) with estimated high and moderate effects were located within 1271 predicted genes. A higher number of mutations were observed on chromosomes Vm1 (588) and Vm3 (428) with the highest frequency on chromosome Vm3 (every 0.07 Mb). Conclusions Our study reiterated the mutagenic utility of EB for inducing SBSs and small indels genome-wide. The knowledge gained from SNP-level profiling of EB-induced mutations can expedite comparative mutation breeding studies in legumes.

Skim exome capture genotyping in wheat.

Next-generation sequencing (NGS) technology advancements continue to reduce the cost of high-throughput genome-wide genotyping for breeding and genetics research. Skim sequencing, which surveys the entire genome at low coverage, has become feasible for quantitative trait locus (QTL) mapping and genomic selection in various crops. However, the genome complexity of allopolyploid crops such as wheat (Triticum aestivum L.) still poses a significant challenge for genome-wide genotyping. Targeted sequencing of the protein-coding regions (i.e., exome) reduces sequencing costs compared to whole genome re-sequencing and can be used for marker discovery and genotyping. We developed a method called skim exome capture (SEC) that combines the strengths of these existing technologies and produces targeted genotyping data while decreasing the cost on a per-sample basis compared to traditional exome capture. Specifically, we fragmented genomic DNA using a tagmentation approach, then enriched those fragments for the low-copy genic portion of the genome using commercial wheat exome baits and multiplexed the sequencing at different levels to achieve desired coverage. We demonstrated that for a library of 48 samples, ∼7-8× target coverage was sufficient for high-quality variant detection. For higher multiplexing levels of 528 and 1056 samples per library, we achieved an average coverage of 0.76× and 0.32×, respectively. Combining these lower coverage SEC sequencing data with genotype imputation using a customized wheat practical haplotype graph database that we developed, we identified hundreds of thousands of high-quality genic variants across the genome. The SEC method can be used for high-resolution QTL mapping, genome-wide association studies, genomic selection, and other downstream applications.

Whole-genome sequencing and variant discovery of Citrus reticulata "Kinnow" from Pakistan.

Citrus is a source of nutritional and medicinal advantages, cultivated worldwide with major groups of sweet oranges, mandarins, grapefruits, kumquats, lemons and limes. Pakistan produces all major citrus groups with mandarin (Citrus reticulata) being the prominent group that includes local commercial cultivars Feutral's Early, Dancy, Honey, and Kinnow. The present study designed to understand the genetic architecture of this unique variety of Citrus reticulata 'Kinnow.' The whole-genome resequencing and variant calling was performed to map the genomic variability that might be responsible for its particular characteristics like taste, seedlessness, juice content, thickness of peel, and shelf-life. A total of 139,436,350 raw sequence reads were generated with 20.9 Gb data in Fastq format having 98% effectiveness and 0.2% base call error rate. Overall, 3,503,033 SNPs, 176,949 MNPs, 323,287 INS, and 333,083 DEL were identified using the GATK4 variant calling pipeline against Citrus clementina. Furthermore, g:Profiler was applied for annotating the newly found variants, harbor genes/transcripts and their involved pathways. A total of 73,864 transcripts harbors 4,336,352 variants, most of the observed variants were predicted in non-coding regions and 1009 transcripts were found well annotated by different databases. Out of total aforementioned transcripts, 588 involved in biological processes, 234 in molecular functions and 167 transcripts in cellular components. In a nutshell, 18,153 high impact variants and 216 genic variants found in the current study,which may be used after its functional validationfor marker-assisted breeding programs of "Kinnow" topropagate its valued traits for the improvement of contemporary citrus varietiesin the region.

A Particular Focus on the Prevalence of α- and β-Thalassemia in Western Sicilian Population from Trapani Province in the COVID-19 Era

Thalassemia is a Mendelian inherited blood disease caused by α- and β-globin gene mutations, known as one of the major health problems of Mediterranean populations. Here, we examined the distribution of α- and β-globin gene defects in the Trapani province population. A total of 2,401 individuals from Trapani province were enrolled from January 2007 to December 2021, and routine methodologies were used for detecting the α- and β-globin genic variants. Appropriate analysis was also performed. Eight mutations in the α globin gene showed the highest frequency in the sample studied; three of these genetic variants represented the 94% of the total α-thalassemia mutations observed, including the -α3.7 deletion (76%), and the tripling of the α gene (12%) and of the α2 point mutation IVS1-5nt (6%). For the β-globin gene, 12 mutations were detected, six of which constituted 83.4% of the total number of β-thalassemia defects observed, including codon β039 (38%), IVS1.6 T > C (15.6%), IVS1.110 G > A (11.8%), IVS1.1 G > A (11%), IVS2.745 C > G (4%), and IVS2.1 G > A (3%). However, the comparison of these frequencies with those detected in the population of other Sicilian provinces did not demonstrate significant differences, but it contrarily revealed a similitude. The data presented in this retrospective study help provide a picture of the prevalence of defects on the α and β-globin genes in the province of Trapani. The identification of mutations in globin genes in a population is required for carrier screening and for an accurate prenatal diagnosis. It is important and necessary to continue promoting public awareness campaigns and screening programs.

Predicting the insecticide-driven mutations in a crop pest insect: Evidence for multiple polymorphisms of acetylcholinesterase gene with potential relevance for resistance to chemicals.

The silverleaf whitefly Bemisia tabaci (Gennadius, 1889) (Homoptera: Aleyrodidae) is a serious invasive herbivorous insect pest worldwide. The excessive use of pesticides has progressively selected B. tabaci specimens, reducing the effectiveness of the treatments, and ultimately ending in the selection of pesticide-resistant strains. The management of this crop pest has thus become challenging owing to the level of resistance to all major classes of recommended insecticides. Here, we used in silico techniques for detecting sequence polymorphisms in ace1 gene from naturally occurring B. tabaci variants, and monitor the presence and frequency of the detected putative mutations from 30 populations of the silverleaf whitefly from Egypt and Pakistan. We found several point mutations in ace1-type acetylcholinesterase (ace1) in the studied B. tabaci variants naturally occurring in the field. By comparing ace1 sequence data from an organophosphate-susceptible and an organophosphate-resistant strains of B. tabaci to ace1 sequence data retrieved from GenBank for that species and to nucleotide polymorphisms from other arthropods, we identified novel mutations that could potentially influence insecticide resistance. Homology modeling and molecular docking analyses were performed to determine if the mutation-induced changes in form 1 acetylcholinesterase (AChE1) structure could confer resistance to carbamate and organophosphate insecticides. Mutations had small effects on binding energy (ΔGb) interactions between mutant AChE1 and insecticides; they altered the conformation of the peripheral anionic site of AChE1, and modified the enzyme surface, and these changes have potential effects on the target-site sensitivity. Altogether, the results from this study provide information on genic variants of B. tabaci ace1 for future monitoring insecticide resistance development and report a potential case of environmentally driven gene variations.

Deleterious mutations are characterized by higher genomic heterozygosity than other genic variants in plant genomes.

Deleterious mutations can reduce the fitness of crop varieties, which limits the plant breeding efficacy. While crop deleterious mutations have been extensively examined, most studies focused on one specific crop with different analyzing methods, which hinders unveiling shared genomic characteristics of deleterious mutations across diverse crops. Here we used standardized approaches to characterize the deleterious mutations in genomes of domesticated inbreeding (i.e., rice, soybean, and tomato) and clonally propagated crops (i.e., grape and pineapple). We found that deleterious mutations are commonly targeted by purifying selection, and are over-presented in a nearly fixed derived allele frequency in the course of plant domestication. Further, a generally negative correlation between genetic load and the artificial selection strength is observed. Importantly, we consistently uncovered the higher derived genomic heterozygosity for deleterious mutations compared to other genic variants. This study broadens our understanding of the evolution of deleterious mutations in plant genomes.

The genetic heterogeneity of colorectal cancer predisposition - guidelines for gene discovery.

BackgroundColorectal cancer (CRC) is a cumulative term applied to a clinically and genetically heterogeneous group of neoplasms that occur in the bowel. Based on twin studies, up to 45 % of the CRC cases may involve a heritable component. Yet, only in 5–10 % of these cases high-penetrant germline mutations are found (e.g. mutations in APC and DNA mismatch repair genes) that result in a familial aggregation and/or an early onset of the disease. Genome-wide association studies have revealed that another ~5 % of the CRC cases may be explained by a cumulative effect of low-penetrant risk factors. Recent attempts to identify novel genetic factors using whole exome and whole genome sequencing has proven to be difficult since the remaining, yet to be discovered, high penetrant CRC predisposing genes appear to be rare. In addition, most of the moderately penetrant candidate genes identified so far have not been confirmed in independent cohorts. Based on literature examples, we here discuss how careful patient and cohort selection, candidate gene and variant selection, and corroborative evidence may be employed to facilitate the discovery of novel CRC predisposing genes.ConclusionsThe picture emerges that the genetic predisposition to CRC is heterogeneous, involving complex interplays between common and rare (inter)genic variants with different penetrances. It is anticipated, however, that the use of large clinically well-defined patient and control datasets, together with improved functional and technical possibilities, will yield enough power to unravel this complex interplay and to generate accurate individualized estimates for the risk to develop CRC.

Targeted exon sequencing fails to identify rare coding variants with large effect in rheumatoid arthritis.

IntroductionAlthough it has been suggested that rare coding variants could explain the substantial missing heritability, very few sequencing studies have been performed in rheumatoid arthritis (RA). We aimed to identify novel functional variants with rare to low frequency using targeted exon sequencing of RA in Korea.MethodsWe analyzed targeted exon sequencing data of 398 genes selected from a multifaceted approach in Korean RA patients (n = 1,217) and controls (n = 717). We conducted a single-marker association test and a gene-based analysis of rare variants. For meta-analysis or enrichment tests, we also used ethnically matched independent samples of Korean genome-wide association studies (GWAS) (n = 4,799) or immunochip data (n = 4,722).ResultsAfter stringent quality control, we analyzed 10,588 variants of 398 genes from 1,934 Korean RA case controls. We identified 13 nonsynonymous variants with nominal association in single-variant association tests. In a meta-analysis, we did not find any novel variant with genome-wide significance for RA risk. Using a gene-based approach, we identified 17 genes with nominal burden signals. Among them, VSTM1 showed the greatest association with RA (P = 7.80 × 10−4). In the enrichment test using Korean GWAS, although the significant signal appeared to be driven by total genic variants, we found no evidence for enriched association of coding variants only with RA.ConclusionsWe were unable to identify rare coding variants with large effect to explain the missing heritability for RA in the current targeted resequencing study. Our study raises skepticism about exon sequencing of targeted genes for complex diseases like RA.Electronic supplementary materialThe online version of this article (doi:10.1186/s13075-014-0447-7) contains supplementary material, which is available to authorized users.

Annotation of functional variation within non-MHC MS susceptibility loci through bioinformatics analysis.

There is a strong and complex genetic component to multiple sclerosis (MS). In addition to variation in the major histocompatibility complex (MHC) region on chromosome 6p21.3, 110 non-MHC susceptibility variants have been identified in Northern Europeans, thus far. The majority of the MS-associated genes are immune related; however, similar to most other complex genetic diseases, the causal variants and biological processes underlying pathogenesis remain largely unknown. We created a comprehensive catalog of putative functional variants that reside within linkage disequilibrium regions of the MS-associated genic variants to guide future studies. Bioinformatics analyses were also conducted using publicly available resources to identify plausible pathological processes relevant to MS and functional hypotheses for established MS-associated variants.

Adaptive phylogeography: functional divergence between haemoglobins derived from different glacial refugia in the bank vole.

Over the years, researchers have used presumptively neutral molecular variation to infer the origins of current species' distributions in northern latitudes (especially Europe). However, several reported examples of genic and chromosomal replacements suggest that end-glacial colonizations of particular northern areas may have involved genetic input from different source populations at different times, coupled with competition and selection. We investigate the functional consequences of differences between two bank vole (Clethrionomys glareolus) haemoglobins deriving from different glacial refugia, one of which partially replaced the other in Britain during end-glacial climate warming. This allows us to examine their adaptive divergence and hence a possible role of selection in the replacement. We determine the amino acid substitution Ser52Cys in the major expressed β-globin gene as the allelic difference. We use structural modelling to reveal that the protein environment renders the 52Cys thiol a highly reactive functional group and we show its reactivity in vitro. We demonstrate that possessing the reactive thiol in haemoglobin increases the resistance of bank vole erythrocytes to oxidative stress. Our study thus provides striking evidence for physiological differences between products of genic variants that spread at the expense of one another during colonization of an area from different glacial refugia.

Genetics of Mitochondrial Myopathies

Pearson syndrome. Clinical spectrum of mitochondrial disorders is broad, so consensus diagnostic criteria for mitochondrial disorders have been proposed for children and adults. Sequence map of human mitochondrial genome with its normal and patho­ genic variants is publicly available. In this review, we will discuss genetic features of mitochondrial myopathies as well as their key findings. Representative syndromes will be described in detail and the less common entities will be summarized in Table 1.

Polimorfismos genéticos da kappa-caseína e da beta-lactoglobulina e produção de leite em bovinos

As variantes gênicas da beta-lactoglobulina (β-LG) e da kappa-caseína (κ-CN) bovinas são associadas à produção, qualidade e características de processamento do leite. O objetivo deste trabalho foi analisar as frequências dos genótipos AA, AB e BB, por meio da técnica de PCR-RFLP, da β-LG e da κ-CN bovinas, e suas associações com a produção de leite (kg leite/dia) em bovinos das raças Girolanda, Holandesa e Jersey. Para a κ-CN, a frequência do genótipo AA foi maior nos animais das raças Holandesa (37%) e Girolanda (63%). Na raça Jersey, houve predomínio do genótipo BB (60%). Para a β-LG, o genótipo AB foi o mais encontrado nas raças Girolanda (54%) e Holandesa (58%), enquanto nos animais da raça Jersey houve predomínio do genótipo BB (45%). Houve associação do alelo B da κ-CN com maior produtividade leiteira nas raças Girolanda e Holandesa, e do alelo A da β-LG com maior produtividade de leite na raça Jersey. As variantes genéticas da κ-CN podem ser usadas como marcadores na seleção para a produtividade leiteira nas raças Girolanda e Holandesa. Para a raça Jersey, as variantes da β-LG seriam mais adequadas para essa seleção.

Gene-Centric Genomewide Association Study via Entropy

Genes are the functional units in most organisms. Compared to genetic variants located outside genes, genic variants are more likely to affect disease risk. The development of the human HapMap project provides an unprecedented opportunity for genetic association studies at the genomewide level for elucidating disease etiology. Currently, most association studies at the single-nucleotide polymorphism (SNP) or the haplotype level rely on the linkage information between SNP markers and disease variants, with which association findings are difficult to replicate. Moreover, variants in genes might not be sufficiently covered by currently available methods. In this article, we present a gene-centric approach via entropy statistics for a genomewide association study to identify disease genes. The new entropy-based approach considers genic variants within one gene simultaneously and is developed on the basis of a joint genotype distribution among genetic variants for an association test. A grouping algorithm based on a penalized entropy measure is proposed to reduce the dimension of the test statistic. Type I error rates and power of the entropy test are evaluated through extensive simulation studies. The results indicate that the entropy test has stable power under different disease models with a reasonable sample size. Compared to single SNP-based analysis, the gene-centric approach has greater power, especially when there is more than one disease variant in a gene. As the genomewide genic SNPs become available, our entropy-based gene-centric approach would provide a robust and computationally efficient way for gene-based genomewide association study.

Gene-centric characteristics of genome-wide association studies.

BackgroundThe high-throughput genotyping chips have contributed greatly to genome-wide association (GWA) studies to identify novel disease susceptibility single nucleotide polymorphisms (SNPs). The high-density chips are designed using two different SNP selection approaches, the direct gene-centric approach, and the indirect quasi-random SNPs or linkage disequilibrium (LD)-based tagSNPs approaches. Although all these approaches can provide high genome coverage and ascertain variants in genes, it is not clear to which extent these approaches could capture the common genic variants. It is also important to characterize and compare the differences between these approaches.Methodology/Principal FindingsIn our study, by using both the Phase II HapMap data and the disease variants extracted from OMIM, a gene-centric evaluation was first performed to evaluate the ability of the approaches in capturing the disease variants in Caucasian population. Then the distribution patterns of SNPs were also characterized in genic regions, evolutionarily conserved introns and nongenic regions, ontologies and pathways. The results show that, no mater which SNP selection approach is used, the current high-density SNP chips provide very high coverage in genic regions and can capture most of known common disease variants under HapMap frame. The results also show that the differences between the direct and the indirect approaches are relatively small. Both have similar SNP distribution patterns in these gene-centric characteristics.Conclusions/SignificanceThis study suggests that the indirect approaches not only have the advantage of high coverage but also are useful for studies focusing on various functional SNPs either in genes or in the conserved regions that the direct approach supports. The study and the annotation of characteristics will be helpful for designing and analyzing GWA studies that aim to identify genetic risk factors involved in common diseases, especially variants in genes and conserved regions.

A gene-centric approach to genome-wide association studies

Genic variants are more likely to alter gene function and affect disease risk than those that occur outside genes. Variants in genes, however, might not be sufficiently covered by the existing approaches to genome-wide association studies. Our analysis of the HapMap ENCODE data indicates that this concern is valid, and that an alternative approach that focuses on genic variants provides a more complete coverage of functionally important regions and a greater genotyping efficiency. We therefore argue that resources should be developed to make gene-centric genome-wide association studies feasible.

Molecular heterosis for heat-sensitive enzyme alleles.

Heat denaturation studies were carried out and revealed hidden genic variants with the same net charge at the Octanol dehydrogenase-1 locus in Drosophila pseudoobscura. Studies of several genetic crosses between strains with different heat-sensitivity alleles showed that the F(1) retained more in vitro enzyme activity after being heat-treated for a specified amount of time at a given temperature than the heat-resistant parent. We call this phenomenon "heterosis for heat-stability of enzyme activity" and discuss its possible molecular mechanism, its relation to maintenance of genetic variation in natural populations, and its bearing on the "classical" and "balance" hypotheses.