CircosVCF: circos visualization of whole-genome sequence variations stored in VCF files

Augmenting tomato functional genomics with a genome-wide induced genetic variation resource

Induced mutations accelerate crop improvement by providing novel disease resistance and yield alleles. However, the alleles with no perceptible phenotype but have an altered function remain hidden in mutagenized plants. The whole-genome sequencing (WGS) of mutagenized individuals uncovers the complete spectrum of mutations in the genome. Genome-wide induced mutation resources can improve the targeted breeding of tomatoes and facilitate functional genomics. In this study, we sequenced 132 doubly ethyl methanesulfonate (EMS)-mutagenized lines of tomato and detected approximately 41 million novel mutations and 5.5 million short InDels not present in the parental cultivar. Approximately 97% of the genome had mutations, including the genes, promoters, UTRs, and introns. More than one-third of genes in the mutagenized population had one or more deleterious mutations predicted by Sorting Intolerant From Tolerant (SIFT). Nearly one-fourth of deleterious genes mapped on tomato metabolic pathways modulate multiple pathway steps. In addition to the reported GC>AT transition bias for EMS, our population also had a substantial number of AT>GC transitions. Comparing mutation frequency among synonymous codons revealed that the most preferred codon is the least mutagenic toward EMS. The validation of a potato leaf-like mutation, reduction in carotenoids in ζ-carotene isomerase mutant fruits, and chloroplast relocation loss in phototropin1 mutant validated the mutation discovery pipeline. Our database makes a large repertoire of mutations accessible to functional genomics studies and breeding of tomatoes.

Analysis of phenotypic changes in high temperature and low pH extreme conditions of Alicyclobacillus sendaiensis PA2 related with the cell wall and sporulation genes

The A. sendaiensis PA2 is a polyextremophile bacterium. In this study, we analyze the A. sendaiensis PA2 genome. The genome was assembled and annotated. The A. sendaiensis PA2 genome structure consists of a 2,956,928 bp long chromosome and 62.77% of G + C content. 3056 CDSs were predicted, and 2921 genes were assigned to a putative function. The ANIm and ANIb value resulted in 97.17% and 96.65%, the DDH value was 75.5%, and the value of TETRA (Z-score) was 0.98. Comparative genomic analyses indicated that three systems are enriched in A. sendaiensis PA2. This strain has phenotypic changes in cell wall during batch culture at 65 °C, pH 5.0 and without carbon and nitrogen source. The presence of unique genes of cell wall and sporulation subsystem could be related to the adaptation of A. sendaiensis PA2 to hostile conditions.

A Field Collection of Indigenous Grapevines as a Valuable Repository for Applied Research

The grapevine is an economically important plant, with a historical connection to the development of human culture. Currently, over 6000 accessions are known as individual grapevine varieties, some of which are important to national heritage, valuable for current viticultural practices, and as genetic resources to maintain plasticity under changing climatic conditions, environmental sustainability, and market demands. Recently, the diversity of cultivated grapevines has declined significantly, due to the increased focus of global wine industries on a few major cultivars. Moreover, due to biotic and abiotic stresses, the wild V. vinifera germplasm's genetic diversity has declined, with some varieties on the verge of extinction. Vitis germplasm conservation can be achieved via either in situ (e.g., protected areas) or Ex situ (e.g., field collections, seed banks, and tissue culture collections) methods. This study aims to highlight the importance of Vitis field bank collections. We demonstrate the research done in the Israeli indigenous Vitis vinifera collection. The multi-layer analysis of the varieties enabled the identification of drought stress-resistant varieties, and suggested a mechanism for this resistance through noting the dramatic phenological differences in foliage development between resistant and sensitive varieties. In addition, we show a general characterization of the varieties via major grape characteristics, including bunch and berry shape, as well as their possible utilization based on their aromatic and phenolic profiles.

Genome-wide SNP and InDel analysis of three Philippine mango species inferred from whole-genome sequencing

Background

The Philippines is among the top 10 major exporters of mango worldwide. However, genomic studies of Philippine mangoes remain largely unexplored and lacking. Here, we sequenced the whole genome of the three Philippine mango species, namely, Mangifera odorata (Huani), Mangifera altissima (Paho), and Mangifera indica “Carabao” variety using Illumina HiSeq 2500, to identify and analyze their genome-wide variants (SNPs and InDels).

Results

The high confidence variants were identified by successfully mapping 93–95% of the quality-filtered reads to the Alphonso and Tommy Atkins mango reference genomes. Using these two currently available mango genomes, most variants were observed in M. odorata (4,353,063 and 4,277,287), followed by M. altissima (3,392,763 and 3,449,917), and lastly, M. indica Carabao (2,755,267 and 2,852,480). Approximately 50, 46, and 38% of the variants were unique in the three Philippine mango genomes. The analysis of variant effects and functional annotation across the three mango species revealed 56,982 variants with high-impact effects mapped onto 37,746 genes, of which 25% were found to be novel. The affected mango genes include those with potential economic importance such as 6945 genes for defense/resistance/immune response, 323 genes for fruit development, and 338 genes for anthocyanin production.

Conclusions

To date, this is the first sequencing effort to comprehensively analyze genome-wide variants essential for the development of genome-wide markers specific to these mango species native to the Philippines. This study provides an important genomic resource that can be used for the genetic improvement of mangoes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-022-00326-3.

Circular RNA hsa_circ_0006091 as a novel biomarker for hepatocellular carcinoma

Circular RNAs (circRNAs) are stable and extensively distributed non-coding RNA molecules that are differentially expressed in liver cancer tissues in the human body. In this study, we aimed to investigate circRNA as a novel candidate biomarker for hepatocellular carcinoma (HCC). For three groups of HCC and neighboring healthy tissues, the differentially expressed circRNAs were identified through high-throughput sequencing analysis. Reverse transcription PCR (RT-PCR) and quantitative polymerase chain reaction (qPCR) were employed for the evaluation of circRNAs that show an elevated expression level in HCC. The obtained results revealed the significantly differential expression of hsa_circ_0006091 in HCC. Then we obtained their target genes through biological analysis, followed by verifying the underlined target genes, and the regulator of G-protein signaling 12 (RGS12) showed an elevated expression level in HCC tissues. Finally, receiver operating characteristic (ROC) curve analysis was conducted on AFP, RGS12, and hsa_circ_0006091, and combined analysis was performed. Furthermore, hsa_circ_0006091 is a novel candidate biomarker for HCC and could improve the diagnostic strategies, prediction, and follow-up of HCC patients. The joint diagnosis of the hsa_circ_0006091&AFP and hsa_circ_0006091&RGS12 has diagnostic significance and can be used as a molecular marker for HCC diagnosis.

Abbreviations: AUC:area under the ROC curve; ROC:Receptor Operating Characteristics; bp:base pair;mRNA:Messenger Ribonucleic acid;ceRNA:Competing endogenous RNA; RT-qPCR: Real time-quantitativen PCR technology; circRNA: circular RNA; HCC:Hepatocellular carcinoma;miRNA:microRNA;KEGG:Kyoto Encyclopedia of Genes and Genomes; RGS12:regulator of G-protein signaling 12; AFP:alpha fetoprotein; ncRNAs:non-coding RNAs; GEO:Gene Expression Omnibus; FDR:false discovery rate

Demographic and ecogeographic factors limit wild grapevine spread at the southern edge of its distribution range

The spatial distribution of plants is constrained by demographic and ecogeographic factors that determine the range and abundance of the species. Wild grapevine (Vitis vinifera ssp. sylvestris) is distributed from Switzerland in the north to Israel in the south. However, little is known about the ecogeographic constraints of this species and its genetic and phenotypic characteristics, especially at the southern edge of its distribution range in the Levant region. In this study, we explore the population structure of southern Levantine wild grapevines and the correlation between demographic and ecogeographic characteristics. Based on our genetic analysis, the wild grapevine populations in this region can be divided into two major subgroups in accordance with a multivariate spatial and ecogeographical clustering model. The identified subpopulations also differ in morphological traits, mainly leaf hairiness which may imply adaptation to environmental stress. The findings suggest that the Upper Jordan River population was spread to the Sea of Galilee area and that a third smaller subpopulation at the south of the Golan Heights may represent a distinguished gene pool or a recent establishment of a new population. A spatial distribution model indicated that distance to water sources, Normalized difference vegetation index, and precipitation are the main environmental factors constraining V. v. sylvestris distribution at its southern distribution range. These factors in addition to limited gene flow between populations prevent further spread of wild grapevines southwards to semi-arid regions.

Muscle allele-specific expression QTLs may affect meat quality traits in Bos indicus

Single nucleotide polymorphisms (SNPs) located in transcript sequences showing allele-specific expression (ASE SNPs) were previously identified in the Longissimus thoracis muscle of a Nelore (Bos indicus) population consisting of 190 steers. Given that the allele-specific expression pattern may result from cis-regulatory SNPs, called allele-specific expression quantitative trait loci (aseQTLs), in this study, we searched for aseQTLs in a window of 1 Mb upstream and downstream from each ASE SNP. After this initial analysis, aiming to investigate variants with a potential regulatory role, we further screened our aseQTL data for sequence similarity with transcription factor binding sites and microRNA (miRNA) binding sites. These aseQTLs were overlapped with methylation data from reduced representation bisulfite sequencing (RRBS) obtained from 12 animals of the same population. We identified 1134 aseQTLs associated with 126 different ASE SNPs. For 215 aseQTLs, one allele potentially affected the affinity of a muscle-expressed transcription factor to its binding site. 162 aseQTLs were predicted to affect 149 miRNA binding sites, from which 114 miRNAs were expressed in muscle. Also, 16 aseQTLs were methylated in our population. Integration of aseQTL with GWAS data revealed enrichment for traits such as meat tenderness, ribeye area, and intramuscular fat . To our knowledge, this is the first report of aseQTLs identification in bovine muscle. Our findings indicate that various cis-regulatory and epigenetic mechanisms can affect multiple variants to modulate the allelic expression. Some of the potential regulatory variants described here were associated with the expression pattern of genes related to interesting phenotypes for livestock. Thus, these variants might be useful for the comprehension of the genetic control of these phenotypes.

A study of normal copy number variations in Israeli population

The population of Israel is ethnically diverse, and individuals from different ethnic groups share specific genetic variations. These variations, which have been passed on from common ancestors, are usually reported in public databases as rare variants. Here, we aimed to identify ethnicity-based benign copy number variants (CNVs) and generate the first Israeli CNV database. We applied a data-mining approach to the results of 10,193 chromosomal microarray tests, of which 2150 tests were from individuals of 13 common ethnic backgrounds (n ≥ 10). We found 165 CNV regions (> 50 kbp) that are unique to specific ethnic groups (uCNVRs). The frequency of more than 19% of these uCNVRs is between 1 and 20% of the common ethnic origin, while their frequency in the overall cohort is between 0.5 and 1.6%. Of these 165 uCNVRs, 98 are reported as variants of unknown significance or as not available in dbVar; we postulate that these uCNVRs should be annotated as either "likely benign" or "benign". The ethnic-specific CNVs extracted in this study will allow geneticists to distinguish between relevant pathogenic genomic aberrations and benign ethnicity-related variations, thus preventing variant misinterpretation that may lead to unnecessary pregnancy terminations.

Review on the Computational Genome Annotation of Sequences Obtained by Next-Generation Sequencing

Next-Generation Sequencing (NGS) has made it easier to obtain genome-wide sequence data and it has shifted the research focus into genome annotation. The challenging tasks involved in annotation rely on the currently available tools and techniques to decode the information contained in nucleotide sequences. This information will improve our understanding of general aspects of life and evolution and improve our ability to diagnose genetic disorders. Here, we present a summary of both structural and functional annotations, as well as the associated comparative annotation tools and pipelines. We highlight visualization tools that immensely aid the annotation process and the contributions of the scientific community to the annotation. Further, we discuss quality-control practices and the need for re-annotation, and highlight the future of annotation.

High-Throughput Genotyping of Resilient Tomato Landraces to Detect Candidate Genes Involved in the Response to High Temperatures

The selection of tolerant varieties is a powerful strategy to ensure highly stable yield under elevated temperatures. In this paper, we report the phenotypic and genotypic characterization of 10 tomato landraces to identify the best performing under high temperatures. The phenotyping of five yield-related traits allowed us to select one genotype that exhibits highly stable yield performances in different environmental conditions. Moreover, a Genotyping-by-Sequencing approach allowed us to explore the genetic variability of the tested genotypes. The high and stable yielding landrace E42 was the most polymorphic one, with ~49% and ~47% private SNPs and InDels, respectively. The effect of 26,113 mutations on proteins’ structure was investigated and it was discovered that 37 had a high impact on the structure of 34 proteins of which some are putatively involved in responses to high temperatures. Additionally, 129 polymorphic sequences aligned against tomato wild species genomes revealed the presence in the genotype E42 of several introgressed regions deriving from S. pimpinellifolium. The position on the tomato map of genes affected by moderate and high impact mutations was also compared with that of known markers/QTLs (Quantitative Trait Loci) associated with reproductive and yield-related traits. The candidate genes/QTLs regulating heat tolerance in the selected landrace E42 could be further investigated to better understand the genetic mechanisms controlling traits for high and stable yield trait under high temperatures.

Different Phenotypes, Similar Genomes: Three Newly Sequenced Fusarium fujikuroi Strains Induce Different Symptoms in Rice Depending on Temperature

Bakanae, caused by the hemibiotrophic fungus Fusarium fujikuroi, is one of the most important diseases of rice and is attributed to up to 75% of losses, depending on the strain and environmental conditions. Some strains cause elongation and thin leaves, whereas others induce stunting and chlorotic seedlings. Differences in symptoms are attributed to genetic differences in the strains. F. fujikuroi strains Augusto2, CSV1, and I1.3 were sequenced with Illumina MiSeq, and pathogenicity trials were conducted on rice cultivar Galileo, which is susceptible to bakanae. By performing gene prediction, single nucleotide polymorphism (SNP) calling, and structural variant analysis with a reference genome, we show how an extremely limited number of polymorphisms in genes not commonly associated with bakanae disease can cause strong differences in phenotype. CSV1 and Augusto2 were particularly close, with only 21,887 SNPs between them, but they differed in virulence, reaction to temperature, induced symptoms, colony morphology and color, growth speed, fumonisin, and gibberellin production. Genes potentially involved in the shift in phenotype were identified. Furthermore, we show how temperature variation may result in different symptoms even in rice plants inoculated with the same F. fujikuroi strain. Moreover, all of the F. fujikuroi strains became more virulent at higher temperatures. Significant differences were likewise observed in gibberellic acid production and in the expression of both fungal and plant gibberellin biosynthetic genes.

Genome Editing in Patient iPSCs Corrects the Most Prevalent USH2A Mutations and Reveals Intriguing Mutant mRNA Expression Profiles

Inherited retinal dystrophies (IRDs) are characterized by progressive photoreceptor degeneration and vision loss. Usher syndrome (USH) is a syndromic IRD characterized by retinitis pigmentosa (RP) and hearing loss. USH is clinically and genetically heterogeneous, and the most prevalent causative gene is USH2A. USH2A mutations also account for a large number of isolated autosomal recessive RP (arRP) cases. This high prevalence is due to two recurrent USH2A mutations, c.2276G>T and c.2299delG. Due to the large size of the USH2A cDNA, gene augmentation therapy is inaccessible. However, CRISPR/Cas9-mediated genome editing is a viable alternative. We used enhanced specificity Cas9 of Streptococcus pyogenes (eSpCas9) to successfully achieve seamless correction of the two most prevalent USH2A mutations in induced pluripotent stem cells (iPSCs) of patients with USH or arRP. Our results highlight features that promote high target efficacy and specificity of eSpCas9. Consistently, we did not identify any off-target mutagenesis in the corrected iPSCs, which also retained pluripotency and genetic stability. Furthermore, analysis of USH2A expression unexpectedly identified aberrant mRNA levels associated with the c.2276G>T and c.2299delG mutations that were reverted following correction. Taken together, our efficient CRISPR/Cas9-mediated strategy for USH2A mutation correction brings hope for a potential treatment for USH and arRP patients.

Overview of Genomic Tools for Circular Visualization in the Next-generation Genomic Sequencing Era

After human genome sequencing and rapid changes in genome sequencing methods, we have entered into the era of rapidly accumulating genome-sequencing data. This has derived the development of several types of methods for representing results of genome sequencing data. Circular genome visual-ization tools are also critical in this area as they provide rapid interpretation and simple visualization of overall data. In the last 15 years, we have seen rapid changes in circular visualization tools after the de-velopment of the circos tool with 1-2 tools published per year. Herein we have summarized and revisited all these tools until the third quarter of 2018.

Genome Sequence, Assembly and Characterization of Two Metschnikowia fructicola Strains Used as Biocontrol Agents of Postharvest Diseases

The yeast Metschnikowia fructicola was reported as an efficient biological control agent of postharvest diseases of fruits and vegetables, and it is the bases of the commercial formulated product "Shemer." Several mechanisms of action by which M. fructicola inhibits postharvest pathogens were suggested including iron-binding compounds, induction of defense signaling genes, production of fungal cell wall degrading enzymes and relatively high amounts of superoxide anions. We assembled the whole genome sequence of two strains of M. fructicola using PacBio and Illumina shotgun sequencing technologies. Using the PacBio, a high-quality draft genome consisting of 93 contigs, with an estimated genome size of approximately 26 Mb, was obtained. Comparative analysis of M. fructicola proteins with the other three available closely related genomes revealed a shared core of homologous proteins coded by 5,776 genes. Comparing the genomes of the two M. fructicola strains using a SNP calling approach resulted in the identification of 564,302 homologous SNPs with 2,004 predicted high impact mutations. The size of the genome is exceptionally high when compared with those of available closely related organisms, and the high rate of homology among M. fructicola genes points toward a recent whole-genome duplication event as the cause of this large genome. Based on the assembled genome, sequences were annotated with a gene description and gene ontology (GO term) and clustered in functional groups. Analysis of CAZymes family genes revealed 1,145 putative genes, and transcriptomic analysis of CAZyme expression levels in M. fructicola during its interaction with either grapefruit peel tissue or Penicillium digitatum revealed a high level of CAZyme gene expression when the yeast was placed in wounded fruit tissue.