TOGGLE: toolbox for generic NGS analyses

Genetic control of the leaf ionome in pearl millet and correlation with root and agromorphological traits

Pearl millet (Pennisetum glaucum) thrives in arid and nutrient-poor environments, establishing its role as a crucial cereal crop for food security in sub-Saharan Africa. Despite its remarkable adaptability, its yields remain below genetic potential, primarily due to limited water and nutrient availability. In this study, we conducted ionomic profiling and genome-wide association studies (GWAS) in field conditions across two growing seasons to unravel the genetic basis of nutrient acquisition in pearl millet. Soil ion content analyses revealed significant differences in nutrient distribution between field sites, while certain ions, such as phosphorus (P) and zinc (Zn), consistently displayed stratified accumulation patterns across years, suggesting stable depth-dependent trends. Evaluation of a genetically diverse panel of inbred lines revealed substantial variation in leaf ion concentrations, with high heritability estimates. Correlations between leaf ion content and root anatomical or agromorphological traits highlighted the intricate interplay between genetic and environmental factors shaping leaf ion accumulation. These analyses also uncovered potential trade-offs in nutrient acquisition strategies. GWAS identified genomic regions associated with leaf ion concentrations, and the integration of genetic and gene expression data facilitated the identification of candidate genes implicated in ion transport and homeostasis. Our findings provide valuable insights into the genetic regulation of nutrient acquisition in pearl millet, offering potential targets for breeding nutrient-efficient and climate-resilient varieties. This study underscores the importance of integrating genetic, physiological, and root architectural traits to enhance agricultural productivity and sustainability in resource-constrained environments.

Early transcriptomic responses of rice leaves to herbivory by Spodoptera frugiperda

During herbivory, chewing insects deposit complex oral secretions (OS) onto the plant wound. Understanding how plants respond to the different cues of herbivory remains an active area of research. In this study, we used an herbivory-mimick experiment to investigate the early transcriptional response of rice plants leaves to wounding, OS, and OS microbiota from Spodoptera frugiperda larvae. Wounding induced a massive early response associated to hormones such as jasmonates. This response switched drastically upon OS treatment indicating the activation of OS specific pathways. When comparing native and dysbiotic OS treatments, we observed few gene regulation. This suggests that in addition to wounding the early response in rice is mainly driven by the insect compounds of the OS rather than microbial. However, microbiota affected genes encoding key phytohormone synthesis enzymes, suggesting an additional modulation of plant response by OS microbiota.

Rapid evolution of an RNA virus to escape recognition by a rice nucleotide-binding and leucine-rich repeat domain immune receptor

Viral diseases are a major limitation for crop production, and their control is crucial for sustainable food supply. We investigated by a combination of functional genetics and experimental evolution the resistance of rice to the rice yellow mottle virus (RYMV), which is among the most devastating rice pathogens in Africa, and the mechanisms underlying the extremely fast adaptation of the virus to its host. We found that the RYMV3 gene that protects rice against the virus codes for a nucleotide-binding and leucine-rich repeat domain immune receptor (NLRs) from the Mla-like clade of NLRs. RYMV3 detects the virus by forming a recognition complex with the viral coat protein (CP). The virus escapes efficiently from detection by mutations in its CP, some of which interfere with the formation of the recognition complex. This study establishes that NLRs also confer in monocotyledonous plants immunity to viruses, and reveals an unexpected functional diversity for NLRs of the Mla clade that were previously only known as fungal disease resistance proteins. In addition, it provides precise insight into the mechanisms by which viruses adapt to plant immunity and gives important knowledge for the development of sustainable resistance against viral diseases of cereals.

Genomic and Bioinformatic Resources for Perennial Fruit Species

In the post-genomic era, data management and development of bioinformatic tools are critical for the adequate exploitation of genomics data. In this review, we address the actual situation for the subset of crops represented by the perennial fruit species. The agronomical singularity of these species compared to plant and crop model species provides significant challenges on the implementation of good practices generally not addressed in other species. Studies are usually performed over several years in non-controlled environments, usage of rootstock is common, and breeders heavily rely on vegetative propagation. A reference genome is now available for all the major species as well as many members of the economically important genera for breeding purposes. Development of pangenome for these species is beginning to gain momentum which will require a substantial effort in term of bioinformatic tool development. The available tools for genome annotation and functional analysis will also be presented.

Population Genome Sequencing of the Scab Fungal Species Venturia inaequalis, Venturia pirina, Venturia aucupariae and Venturia asperata

The Venturia genus comprises fungal species that are pathogens on Rosaceae host plants, including V. inaequalis and V. asperata on apple, V. aucupariae on sorbus and V. pirina on pear. Although the genetic structure of V. inaequalis populations has been investigated in detail, genomic features underlying these subdivisions remain poorly understood. Here, we report whole genome sequencing of 87 Venturia strains that represent each species and each population within V. inaequalis. We present a PacBio genome assembly for the V. inaequalis EU-B04 reference isolate. The size of selected genomes was determined by flow cytometry, and varied from 45 to 93 Mb. Genome assemblies of V. inaequalis and V. aucupariae contain a high content of transposable elements (TEs), most of which belong to the Gypsy or Copia LTR superfamilies and have been inactivated by Repeat-Induced Point mutations. The reference assembly of V. inaequalis presents a mosaic structure of GC-equilibrated regions that mainly contain predicted genes and AT-rich regions, mainly composed of TEs. Six pairs of strains were identified as clones. Single-Nucleotide Polymorphism (SNP) analysis between these clones revealed a high number of SNPs that are mostly located in AT-rich regions due to misalignments and allowed determining a false discovery rate. The availability of these genome sequences is expected to stimulate genetics and population genomics research of Venturia pathogens. Especially, it will help understanding the evolutionary history of Venturia species that are pathogenic on different hosts, a history that has probably been substantially influenced by TEs.

Managing genomic variant calling workflows with Swift/T

Bioinformatics research is frequently performed using complex workflows with multiple steps, fans, merges, and conditionals. This complexity makes management of the workflow difficult on a computer cluster, especially when running in parallel on large batches of data: hundreds or thousands of samples at a time. Scientific workflow management systems could help with that. Many are now being proposed, but is there yet the “best” workflow management system for bioinformatics? Such a system would need to satisfy numerous, sometimes conflicting requirements: from ease of use, to seamless deployment at peta- and exa-scale, and portability to the cloud. We evaluated Swift/T as a candidate for such role by implementing a primary genomic variant calling workflow in the Swift/T language, focusing on workflow management, performance and scalability issues that arise from production-grade big data genomic analyses. In the process we introduced novel features into the language, which are now part of its open repository. Additionally, we formalized a set of design criteria for quality, robust, maintainable workflows that must function at-scale in a production setting, such as a large genomic sequencing facility or a major hospital system. The use of Swift/T conveys two key advantages. (1) It operates transparently in multiple cluster scheduling environments (PBS Torque, SLURM, Cray aprun environment, etc.), thus a single workflow is trivially portable across numerous clusters. (2) The leaf functions of Swift/T permit developers to easily swap executables in and out of the workflow, which makes it easy to maintain and to request resources optimal for each stage of the pipeline. While Swift/T’s data-level parallelism eliminates the need to code parallel analysis of multiple samples, it does make debugging more difficult, as is common for implicitly parallel code. Nonetheless, the language gives users a powerful and portable way to scale up analyses in many computing architectures. The code for our implementation of a variant calling workflow using Swift/T can be found on GitHub at https://github.com/ncsa/Swift-T-Variant-Calling, with full documentation provided at http://swift-t-variant-calling.readthedocs.io/en/latest/.

Development and evaluation of a genome-wide Coffee 8.5K SNP array and its application for high-density genetic mapping and for investigating the origin of Coffea arabica L

Coffee species such as Coffea canephora P. (Robusta) and C. arabica L. (Arabica) are important cash crops in tropical regions around the world. C. arabica is an allotetraploid (2n = 4x = 44) originating from a hybridization event of the two diploid species C. canephora and C. eugenioides (2n = 2x = 22). Interestingly, these progenitor species harbour a greater level of genetic variability and are an important source of genes to broaden the narrow Arabica genetic base. Here, we describe the development, evaluation and use of a single-nucleotide polymorphism (SNP) array for coffee trees. A total of 8580 unique and informative SNPs were selected from C. canephora and C. arabica sequencing data, with 40% of the SNP located in annotated genes. In particular, this array contains 227 markers associated to 149 genes and traits of agronomic importance. Among these, 7065 SNPs (~82.3%) were scorable and evenly distributed over the genome with a mean distance of 54.4 Kb between markers. With this array, we improved the Robusta high-density genetic map by adding 1307 SNP markers, whereas 945 SNPs were found segregating in the Arabica mapping progeny. A panel of C. canephora accessions was successfully discriminated and over 70% of the SNP markers were transferable across the three species. Furthermore, the canephora-derived subgenome of C. arabica was shown to be more closely related to C. canephora accessions from northern Uganda than to other current populations. These validated SNP markers and high-density genetic maps will be useful to molecular genetics and for innovative approaches in coffee breeding.

Origin of the Aromatic Group of Cultivated Rice (Oryza sativa L.) Traced to the Indian Subcontinent

The aromatic group of Asian cultivated rice is a distinct population with considerable genetic diversity on the Indian subcontinent and includes the popular Basmati types characterized by pleasant fragrance. Genetic and phenotypic associations with other cultivated groups are ambiguous, obscuring the origin of the aromatic population. From analysis of genome-wide diversity among over 1,000 wild and cultivated rice accessions, we show that aromatic rice originated in the Indian subcontinent from hybridization between a local wild population and examples of domesticated japonica that had spread to the region from their own center of origin in East Asia. Most present-day aromatic accessions have inherited their cytoplasm along with 29-47% of their nuclear genome from the local Indian rice. We infer that the admixture occurred 4,000-2,400 years ago, soon after japonica rice reached the region. We identify aus as the original crop of the Indian subcontinent, indica and japonica as later arrivals, and aromatic a specific product of local agriculture. These results prompt a reappraisal of our understanding of the emergence and development of rice agriculture in the Indian subcontinent.

The Rise and Fall of African Rice Cultivation Revealed by Analysis of 246 New Genomes

African rice (Oryza glaberrima) was domesticated independently from Asian rice. The geographical origin of its domestication remains elusive. Using 246 new whole-genome sequences, we inferred the cradle of its domestication to be in the Inner Niger Delta. Domestication was preceded by a sharp decline of most wild populations that started more than 10,000 years ago. The wild population collapse occurred during the drying of the Sahara. This finding supports the hypothesis that depletion of wild resources in the Sahara triggered African rice domestication. African rice cultivation strongly expanded 2,000 years ago. During the last 5 centuries, a sharp decline of its cultivation coincided with the introduction of Asian rice in Africa. A gene, PROG1, associated with an erect plant architecture phenotype, showed convergent selection in two rice cultivated species, Oryza glaberrima from Africa and Oryza sativa from Asia. In contrast, a shattering gene, SH5, showed selection signature during African rice domestication, but not during Asian rice domestication. Overall, our genomic data revealed a complex history of African rice domestication influenced by important climatic changes in the Saharan area, by the expansion of African agricultural society, and by recent replacement by another domesticated species.

The Number of Target Molecules of the Amplification Step Limits Accuracy and Sensitivity in Ultradeep-Sequencing Viral Population Studies

The invention of next-generation sequencing (NGS) techniques marked the coming of a new era in the detection of the genetic diversity of intrahost viral populations. A good understanding of the genetic structure of these populations requires, first, the ability to identify the different isolates or variants and, second, the ability to accurately quantify them. However, the initial amplification step of NGS studies can impose potential quantitative biases, modifying the variant relative frequencies. In particular, the number of target molecules (NTM) used during the amplification step is vastly overlooked although of primary importance, as it sets the limit of the accuracy and sensitivity of the sequencing procedure. In the present article, we investigated quantitative biases in an NGS study of populations of a multipartite single-stranded DNA (ssDNA) virus at different steps of the procedure. We studied 20 independent populations of the ssDNA virus faba bean necrotic stunt virus (FBNSV) in two host plants, Vicia faba and Medicago truncatula FBNSV is a multipartite virus composed of eight genomic segments, whose specific and host-dependent relative frequencies are defined as the "genome formula." Our results show a significant distortion of the FBNSV genome formula after the amplification and sequencing steps. We also quantified the genetic bottleneck occurring at the amplification step by documenting the NTM of two genomic segments of FBNSV. We argue that the NTM must be documented and carefully considered when determining the sensitivity and accuracy of data from NGS studies.IMPORTANCE The advent of next-generation sequencing (NGS) techniques now enables study of the genetic diversity of viral populations. A good understanding of the genetic structure of these populations first requires the ability to identify the different isolates or variants and second requires the ability to accurately quantify them. Prior to sequencing, viral genomes need to be amplified, a step that potentially imposes quantitative biases and modifies the viral population structure. In particular, the number of target molecules (NTM) used during the amplification step is of primary importance, as it sets the limit of the accuracy and sensitivity of the sequencing procedure. In this work, we used 20 replicated populations of the multipartite faba bean necrotic stunt virus (FBNSV) to estimate the various limitations of ultradeep-sequencing studies performed on intrahost viral populations. We report quantitative biases during rolling-circle amplification and the NTM of two genomic segments of FBNSV.

GenomeVIP: a cloud platform for genomic variant discovery and interpretation

Identifying genomic variants is a fundamental first step toward the understanding of the role of inherited and acquired variation in disease. The accelerating growth in the corpus of sequencing data that underpins such analysis is making the data-download bottleneck more evident, placing substantial burdens on the research community to keep pace. As a result, the search for alternative approaches to the traditional “download and analyze” paradigm on local computing resources has led to a rapidly growing demand for cloud-computing solutions for genomics analysis. Here, we introduce the Genome Variant Investigation Platform (GenomeVIP), an open-source framework for performing genomics variant discovery and annotation using cloud- or local high-performance computing infrastructure. GenomeVIP orchestrates the analysis of whole-genome and exome sequence data using a set of robust and popular task-specific tools, including VarScan, GATK, Pindel, BreakDancer, Strelka, and Genome STRiP, through a web interface. GenomeVIP has been used for genomic analysis in large-data projects such as the TCGA PanCanAtlas and in other projects, such as the ICGC Pilots, CPTAC, ICGC-TCGA DREAM Challenges, and the 1000 Genomes SV Project. Here, we demonstrate GenomeVIP's ability to provide high-confidence annotated somatic, germline, and de novo variants of potential biological significance using publicly available data sets.

Mutations in Rice yellow mottle virus Polyprotein P2a Involved in RYMV2 Gene Resistance Breakdown

Rice yellow mottle virus (RYMV) is one of the major diseases of rice in Africa. The high resistance of the Oryza glaberrima Tog7291 accession involves a null allele of the RYMV2 gene, whose ortholog in Arabidopsis, CPR5, is a transmembrane nucleoporin involved in effector-triggered immunity. To optimize field deployment of the RYMV2 gene and improve its durability, which is often a weak point in varietal resistance, we analyzed its efficiency toward RYMV isolates representing the genetic diversity of the virus and the molecular basis of resistance breakdown. Tog7291 resistance efficiency was highly variable depending on the isolate used, with infection rates ranging from 0 to 98% of plants. Back-inoculation experiments indicated that infection cases were not due to an incomplete resistance phenotype but to the emergence of resistance-breaking (RB) variants. Interestingly, the capacity of the virus to overcome Tog7291 resistance is associated with a polymorphism at amino-acid 49 of the VPg protein which also affects capacity to overcome the previously studied RYMV1 resistance gene. This polymorphism appeared to be a main determinant of the emergence of RB variants. It acts independently of the resistance gene and rather reflects inter-species adaptation with potential consequences for the durability of resistance. RB mutations were identified by full-length or partial sequencing of the RYMV genome in infected Tog7291 plants and were validated by directed mutagenesis of an infectious viral clone. We found that Tog7291 resistance breakdown involved mutations in the putative membrane anchor domain of the polyprotein P2a. Although the precise effect of these mutations on rice/RYMV interaction is still unknown, our results offer a new perspective for the understanding of RYMV2 mediated resistance mechanisms. Interestingly, in the susceptible IR64 variety, RB variants showed low infectivity and frequent reversion to the wild-type genotype, suggesting that Tog7291 resistance breakdown is associated with a major loss of viral fitness in normally susceptible O. sativa varieties. Despite the high frequency of resistance breakdown in controlled conditions, this loss of fitness is an encouraging element with regards to RYMV2 resistance durability.