Pan-genome Studies Research Articles

Bacterial Pangenome: A Review on the Current Strategies, Tools and Applications

The genomic research has grown by the advancement of sequencing technologies, which aim to capture the "totality" of life diversity by progressing from single strain to species or even higher taxa. Out of this scientific progress and inquisitiveness, the term "pangenome" was created. The goal of a pangenome is to represent all of a species genetic variation. It divides the genome into core, accessory and unique genes. Since the introduction of pangenomics, many software programs have been actively used for pangenome analysis. This review focuses on providing a snapshot of useful tools with their use, execution and computational details. In addition to tools for pangenome, the components, approaches, mechanism of variation, applications and challenges associated with pangenome analysis are also discussed in this article. Here, we have also discussed the increasing trend of bacterial pangenome studies with few studies and their significant findings. It has been proved that pangenomics is useful for understanding the clade of bacteria that helps in developing strategies for various application based on biological similarities and differences.

A comprehensive evaluation of the potential of three next-generation short-read-based plant pan-genome construction strategies for the identification of novel non-reference sequence.

Pan-genome studies are important for understanding plant evolution and guiding the breeding of crops by containing all genomic diversity of a certain species. Three short-read-based strategies for plant pan-genome construction include iterative individual, iteration pooling, and map-to-pan. Their performance is very different under various conditions, while comprehensive evaluations have yet to be conducted nowadays. Here, we evaluate the performance of these three pan-genome construction strategies for plants under different sequencing depths and sample sizes. Also, we indicate the influence of length and repeat content percentage of novel sequences on three pan-genome construction strategies. Besides, we compare the computational resource consumption among the three strategies. Our findings indicate that map-to-pan has the greatest recall but the lowest precision. In contrast, both two iterative strategies have superior precision but lower recall. Factors of sample numbers, novel sequence length, and the percentage of novel sequences' repeat content adversely affect the performance of all three strategies. Increased sequencing depth improves map-to-pan's performance, while not affecting the other two iterative strategies. For computational resource consumption, map-to-pan demands considerably more than the other two iterative strategies. Overall, the iterative strategy, especially the iterative pooling strategy, is optimal when the sequencing depth is less than 20X. Map-to-pan is preferable when the sequencing depth exceeds 20X despite its higher computational resource consumption.

Distinct genomic contexts predict gene presence-absence variation in different pathotypes of Magnaporthe oryzae.

Fungi use the accessory gene content of their pangenomes to adapt to their environments. While gene presence-absence variation contributes to shaping accessory gene reservoirs, the genomic contexts that shape these events remain unclear. Since pangenome studies are typically species-wide and do not analyze different populations separately, it is yet to be uncovered whether presence-absence variation patterns and mechanisms are consistent across populations. Fungal plant pathogens are useful models for studying presence-absence variation because they rely on it to adapt to their hosts, and members of a species often infect distinct hosts. We analyzed gene presence-absence variation in the blast fungus, Magnaporthe oryzae (syn. Pyricularia oryzae), and found that presence-absence variation genes involved in host-pathogen and microbe-microbe interactions may drive the adaptation of the fungus to its environment. We then analyzed genomic and epigenomic features of presence-absence variation and observed that proximity to transposable elements, gene GC content, gene length, expression level in the host, and histone H3K27me3 marks were different between presence-absence variation genes and conserved genes. We used these features to construct a model that was able to predict whether a gene is likely to experience presence-absence variation with high precision (86.06%) and recall (92.88%) in M. oryzae. Finally, we found that presence-absence variation genes in the rice and wheat pathotypes of M. oryzae differed in their number and their genomic context. Our results suggest that genomic and epigenomic features of gene presence-absence variation can be used to better understand and predict fungal pangenome evolution. We also show that substantial intra-species variation can exist in these features.

Novel vaccine candidates of Bordetella pertussis identified by reverse vaccinology

Whooping cough is a disease caused by Bordetella pertussis, whose morbidity has increased, motivating the improvement of current vaccines. Reverse vaccinology is a strategy that helps identify proteins with good characteristics fast and with fewer resources. In this work, we applied reverse vaccinology to study the B. pertussis proteome and pangenome with several in-silico tools. We analyzed the B. pertussis Tohama I proteome with NERVE software and compared 234 proteins with B. parapertussis, B. bronchiseptica, and B. holmessi. VaxiJen was used to calculate an antigenicity value; our threshold was 0.6, selecting 84 proteins. The candidates were depurated and grouped in eight family proteins to select representative candidates, according to bibliographic information and their immunological response predicted with ABCpred, Bcepred, IgPred, and C-ImmSim. Additionally, a pangenome study was conducted with 603 B. pertussis strains and PanRV software, identifying 3421 core proteins that were analyzed to select the best candidates. Finally, we selected 15 proteins from the proteome study and seven proteins from the pangenome analysis as good vaccine candidates.

Pangenome and genomic signatures linked to the dominance of the lineage-4 of Mycobacterium tuberculosis isolated from extrapulmonary tuberculosis patients in western Ethiopia.

The lineage 4 (L4) of Mycobacterium tuberculosis (MTB) is not only globally prevalent but also locally dominant, surpassing other lineages, with lineage 2 (L2) following in prevalence. Despite its widespread occurrence, factors influencing the expansion of L4 and its sub-lineages remain poorly understood both at local and global levels. Therefore, this study aimed to conduct a pan-genome and identify genomic signatures linked to the elevated prevalence of L4 sublineages among extrapulmonary TB (EPTB) patients in western Ethiopia. A cross-sectional study was conducted at an institutional level involving confirmed cases of extrapulmonary tuberculosis (EPTB) patients from August 5, 2018, to December 30, 2019. A total of 75 MTB genomes, classified under lineage 4 (L4), were used for conducting pan-genome and genome-wide association study (GWAS) analyses. After a quality check, variants were identified using MTBseq, and genomes were de novo assembled using SPAdes. Gene prediction and annotation were performed using Prokka. The pan-genome was constructed using GET_HOMOLOGUES, and its functional analysis was carried out with the Bacterial Pan-Genome Analysis tool (BPGA). For GWAS analysis, Scoary was employed with Benjamini-Hochberg correction, with a significance threshold set at p-value ≤ 0.05. The analysis revealed a total of 3,270 core genes, predominantly associated with orthologous groups (COG) functions, notably in the categories of '[R] General function prediction only' and '[I] Lipid transport and metabolism'. Conversely, functions related to '[N] Cell motility' and '[Q] Secondary metabolites biosynthesis, transport, and catabolism' were primarily linked to unique and accessory genes. The pan-genome of MTB L4 was found to be open. Furthermore, the GWAS study identified genomic signatures linked to the prevalence of sublineages L4.6.3 and L4.2.2.2. Apart from host and environmental factors, the sublineage of L4 employs distinct virulence factors for successful dissemination in western Ethiopia. Given that the functions of these newly identified genes are not well understood, it is advisable to experimentally validate their roles, particularly in the successful transmission of specific L4 sublineages over others.

Comparative genomics profiling revealed multi-stress responsive roles of the CC-NBS-LRR genes in three mango cultivars.

The nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family is the largest group of disease resistance (R) genes in plants and is active in response to viruses, bacteria, and fungi usually involved in effector-triggered immunity (ETI). Pangenome-wide studies allow researchers to analyze the genetic diversity of multiple species or their members simultaneously, providing a comprehensive understanding of the evolutionary relationships and diversity present among them. The draft pan-genome of three Mangifera indica cultivars (Alphonso, Hong Xiang Ya, and Tommy atkins) was constructed and Presence/absence variants (PAVs) were filtered through the ppsPCP pipeline. As a result, 2823 genes and 5907 PAVs from H. Xiang Ya, and 1266 genes and 2098 PAVs from T. atkins were added to the reference genome. For the identification of CC-NBS-LRR (CNL) genes in these mango cultivars, this draft pan-genome study has successfully identified 47, 27, and 36 members in Alphonso, H. Xiang Ya, and T. atkins respectively. The phylogenetic analysis divided MiCNL proteins into four distinct subgroups. All MiCNL genes are unevenly distributed on chromosomes. Both tandem and segmental duplication events played a significant role in the expansion of the CNL gene family. These genes contain cis-elements related to light, stress, hormone, and development. The analysis of protein-protein interactions (PPI) revealed that MiCNL proteins interacted with other defense-responsive proteins. Gene Ontology (GO) analysis indicated that MiCNL genes play a role in defense mechanisms within the organism. The expression level of the identified genes in fruit peel was observed under disease and cold stress which showed that Mi_A_CNL13 and 14 were up-regulated while Mi_A_CNL15, 25, 30, 31, and 40 were down-regulated in disease stress. On the other hand, Mi_A_CNL2, 14, 41, and 45 were up-regulated and Mi_A_CNL47 is down-regulated in cold stress. Subsequently, the Random Forest (RF) classifier was used to assess the multi-stress response of MiCNLs. It was found that Mi_A_CNL14 is a gene that responds to multiple stress conditions. The CNLs have similar protein structures which show that they are involved in the same function. The above findings provide a foundation for a deeper understanding of the functional characteristics of the mango CNL gene family.

Triticeae crop genome biology: an endless frontier.

Triticeae, the wheatgrass tribe, includes several major cereal crops and their wild relatives. Major crops within the Triticeae are wheat, barley, rye, and oat, which are important for human consumption, animal feed, and rangeland protection. Species within this tribe are known for their large genomes and complex genetic histories. Powered by recent advances in sequencing technology, researchers worldwide have made progress in elucidating the genomes of Triticeae crops. In addition to assemblies of high-quality reference genomes, pan-genome studies have just started to capture the genomic diversities of these species, shedding light on our understanding of the genetic basis of domestication and environmental adaptation of Triticeae crops. In this review, we focus on recent signs of progress in genome sequencing, pan-genome analyses, and resequencing analysis of Triticeae crops. We also propose future research avenues in Triticeae crop genomes, including identifying genome structure variations, the association of genomic regions with desired traits, mining functions of the non-coding area, introgression of high-quality genes from wild Triticeae resources, genome editing, and integration of genomic resources.

Deciphering the complex cotton genome for improving fiber traits and abiotic stress resilience in sustainable agriculture.

Understanding the complex cotton genome is of paramount importance in devising a strategy for sustainable agriculture. Cotton is probably the most economically important cash crop known for its cellulose-rich fiber content. The cotton genome has become an ideal model for deciphering polyploidization due to its polyploidy, setting it apart from other major crops. However, the main challenge in understanding the functional and regulatory functions of many genes in cotton is still the complex cotton polyploidy genome, which is not limited to a single role. Cotton production is vulnerable to the sensitive effects of climate change, which can alter or aggravate soil, pests, and diseases. Thus, conventional plant breeding coupled with advanced technologies has led to substantial progress being made in cotton production. In the frontier areas of genomics research, cotton genomics has gained momentum accomplished by robust high-throughput sequencing platforms combined with novel computational tools to make the cotton genome more tractable. Advances in long-read sequencing have allowed for the generation of the complete set of cotton gene transcripts giving incisive scientific knowledge in cotton improvement. In contrast, the integration of the latest sequencing platforms has been used to generate multiple high-quality reference genomes in diploid and tetraploid cotton. While pan-genome and 3D genomic studies are still in the early stages in cotton, it is anticipated that rapid advances in sequencing, assembly algorithms, and analysis pipelines will have a greater impact on advanced cotton research. This review article briefly compiles substantial contributions in different areas of the cotton genome, which include genome sequencing, genes, and their molecular regulatory networks in fiber development and stress tolerance mechanism. This will greatly help us in understanding the robust genomic organization which in turn will help unearth candidate genes for functionally important agronomic traits.

Chloroplast Pan-Genomes and Comparative Transcriptomics Reveal Genetic Variation and Temperature Adaptation in the Cucumber.

Although whole genome sequencing, genetic variation mapping, and pan-genome studies have been done on a large group of cucumber nuclear genomes, organelle genome information is largely unclear. As an important component of the organelle genome, the chloroplast genome is highly conserved, which makes it a useful tool for studying plant phylogeny, crop domestication, and species adaptation. Here, we have constructed the first cucumber chloroplast pan-genome based on 121 cucumber germplasms, and investigated the genetic variations of the cucumber chloroplast genome through comparative genomic, phylogenetic, haplotype, and population genetic structure analysis. Meanwhile, we explored the changes in expression of cucumber chloroplast genes under high- and low-temperature stimulation via transcriptome analysis. As a result, a total of 50 complete chloroplast genomes were successfully assembled from 121 cucumber resequencing data, ranging in size from 156,616-157,641 bp. The 50 cucumber chloroplast genomes have typical quadripartite structures, consisting of a large single copy (LSC, 86,339-86,883 bp), a small single copy (SSC, 18,069-18,363 bp), and two inverted repeats (IRs, 25,166-25,797 bp). Comparative genomic, haplotype, and population genetic structure results showed that there is more genetic variation in Indian ecotype cucumbers compared to other cucumber cultivars, which means that many genetic resources remain to be explored in Indian ecotype cucumbers. Phylogenetic analysis showed that the 50 cucumber germplasms could be classified into 3 types: East Asian, Eurasian + Indian, and Xishuangbanna + Indian. The transcriptomic analysis showed that matK were significantly up-regulated under high- and low-temperature stresses, further demonstrating that cucumber chloroplasts respond to temperature adversity by regulating lipid metabolism and ribosome metabolism. Further, accD has higher editing efficiency under high-temperature stress, which may contribute to the heat tolerance. These studies provide useful insight into genetic variation in the chloroplast genome, and established the foundation for exploring the mechanisms of temperature-stimulated chloroplast adaptation.

Wide-genome selection of lactic acid bacteria harboring genes that promote the elimination of antinutritional factors.

Anti-nutritional factors (ANFs) substances in plant products, such as indigestible non-starchy polysaccharides (α-galactooligosaccharides, α-GOS), phytate, tannins, and alkaloids can impede the absorption of many critical nutrients and cause major physiological disorders. To enhance silage quality and its tolerance threshold for humans as well as other animals, ANFs must be reduced. This study aims to identify and compare the bacterial species/strains that are potential use for industrial fermentation and ANFs reduction. A pan-genome study of 351 bacterial genomes was performed, and binary data was processed to quantify the number of genes involved in the removal of ANFs. Among four pan-genomes analysis, all 37 tested Bacillus subtilis genomes had one phytate degradation gene, while 91 out of 150 Enterobacteriacae genomes harbor at least one genes (maximum three). Although, no gene encoding phytase detected in genomes of Lactobacillus and Pediococcus species, they have genes involving indirectly in metabolism of phytate-derivatives to produce Myo-inositol, an important compound in animal cells physiology. In contrast, genes related to production of lectin, tannase and saponin degrading enzyme did not include in genomes of B. subtilis and Pediococcus species. Our findings suggest a combination of bacterial species and/or unique strains in fermentation, for examples, two Lactobacillus strains (DSM 21115 and ATCC 14869) with B. subtilis SRCM103689, would maximize the efficiency in reducing the ANFs concentration. In conclusion, this study provides insights into bacterial genomes analysis for maximizing nutritional value in plant-based food. Further investigations of gene numbers and repertories correlated to metabolism of different ANFs will help clarifying the efficiency of time consuming and food qualities.

FrangiPANe, a tool for creating a panreference using left behind reads.

We present here FrangiPANe, a pipeline developed to build panreference using short reads through a map-then-assemble strategy. Applying it to 248 African rice genomes using an improved CG14 reference genome, we identified an average of 8 Mb of new sequences and 5290 new contigs per individual. In total, 1.4 G of new sequences, consisting of 1 306 676 contigs, were assembled. We validated 97.7% of the contigs of the TOG5681 cultivar individual assembly from short reads on a newly long reads genome assembly of the same TOG5681 cultivar. FrangiPANe also allowed the anchoring of 31.5% of the new contigs within the CG14 reference genome, with a 92.5% accuracy at 2kb span. We annotated in addition 3252 new genes absent from the reference. FrangiPANe was developed as a modular and interactive application to simplify the construction of a panreference using the map-then-assemble approach. It is available as a Docker image containing (i) a Jupyter notebook centralizing codes, documentation and interactive visualization of results, (ii) python scripts and (iii) all the software and libraries requested for each step of the analysis. We foreseen our approach will help leverage large-scale illumina dataset for pangenome studies in GWAS or detection of selection.

Efficient detection and assembly of non-reference DNA sequences with synthetic long reads.

Recent pan-genome studies have revealed an abundance of DNA sequences in human genomes that are not present in the reference genome. A lion’s share of these non-reference sequences (NRSs) cannot be reliably assembled or placed on the reference genome. Improvements in long-read and synthetic long-read (aka linked-read) technologies have great potential for the characterization of NRSs. While synthetic long reads require less input DNA than long-read datasets, they are algorithmically more challenging to use. Except for computationally expensive whole-genome assembly methods, there is no synthetic long-read method for NRS detection. We propose a novel integrated alignment-based and local assembly-based algorithm, Novel-X, that uses the barcode information encoded in synthetic long reads to improve the detection of such events without a whole-genome de novo assembly. Our evaluations demonstrate that Novel-X finds many non-reference sequences that cannot be found by state-of-the-art short-read methods. We applied Novel-X to a diverse set of 68 samples from the Polaris HiSeq 4000 PGx cohort. Novel-X discovered 16 691 NRS insertions of size > 300 bp (total length 18.2 Mb). Many of them are population specific or may have a functional impact.

High-quality genome assembly and pan-genome studies facilitate genetic discovery in mung bean and its improvement

Mung bean is an economically important legume crop species that is used as a food, consumed as a vegetable, and used as an ingredient and even as a medicine. To explore the genomic diversity of mung bean, we assembled a high-quality reference genome (Vrad_JL7) that was ∼479.35 Mb in size, with a contig N50 length of 10.34 Mb. A total of 40,125 protein-coding genes were annotated, representing ∼96.9% of the genetic region. We also sequenced 217 accessions, mainly landraces and cultivars from China, and identified 2,229,343 high-quality single-nucleotide polymorphisms (SNPs). Population structure revealed that the Chinese accessions diverged into two groups and were distinct from non-Chinese lines. Genetic diversity analysis based on genomic data from 750 accessions in 23 countries supported the hypothesis that mung bean was first domesticated in south Asia and introduced to east Asia probably through the Silk Road. We constructed the first pan-genome of mung bean germplasm and assembled 287.73 Mb of non-reference sequences. Among the genes, 83.1% were core genes and 16.9% were variable. Presence/absence variation (PAV) events of nine genes involved in the regulation of the photoperiodic flowering pathway were identified as being under selection during the adaptation process to promote early flowering in the spring. Genome-wide association studies (GWASs) revealed 2,912 SNPs and 259 gene PAV events associated with 33 agronomic traits, including a SNP in the coding region of the SWEET10 homolog (jg24043) involved in crude starch content and a PAV event in a large fragment containing 11 genes for color-related traits. This high-quality reference genome and pan-genome will provide insights into mung bean breeding.

Comparative Genomic Analysis of Vibrio cincinnatiensis Provides Insights into Genetic Diversity, Evolutionary Dynamics, and Pathogenic Traits of the Species.

Vibrio cincinnatiensis is a poorly understood pathogenic Vibrio species, and the underlying mechanisms of its genetic diversity, genomic plasticity, evolutionary dynamics, and pathogenicity have not yet been comprehensively investigated. Here, a comparative genomic analysis of V. cincinnatiensis was constructed. The open pan-genome with a flexible gene repertoire exhibited genetic diversity. The genomic plasticity and stability were characterized by the determinations of diverse mobile genetic elements (MGEs) and barriers to horizontal gene transfer (HGT), respectively. Evolutionary divergences were exhibited by the difference in functional enrichment and selective pressure between the different components of the pan-genome. The evolution on the Chr I and Chr II core genomes was mainly driven by purifying selection. Predicted essential genes in V. cincinnatiensis were mainly found in the core gene families on Chr I and were subject to stronger evolutionary constraints. We identified diverse virulence-related elements, including the gene clusters involved in encoding flagella, secretion systems, several pili, and scattered virulence genes. Our results indicated the pathogenic potential of V. cincinnatiensis and highlighted that HGT events from other Vibrio species promoted pathogenicity. This pan-genome study provides comprehensive insights into this poorly understood species from the genomic perspective.

Genotype to ecotype in niche environments: adaptation of Arthrobacter to carbon availability and environmental conditions

Niche environmental conditions influence both the structure and function of microbial communities and the cellular function of individual strains. The terrestrial subsurface is a dynamic and diverse environment that exhibits specific biogeochemical conditions associated with depth, resulting in distinct environmental niches. Here, we present the characterization of seven distinct strains belonging to the genus Arthrobacter isolated from varying depths of a single sediment core and associated groundwater from an adjacent well. We characterized genotype and phenotype of each isolate to connect specific cellular functions and metabolisms to ecotype. Arthrobacter isolates from each ecotype demonstrated functional and genomic capacities specific to their biogeochemical conditions of origin, including laboratory-demonstrated characterization of salinity tolerance and optimal pH, and genes for utilization of carbohydrates and other carbon substrates. Analysis of the Arthrobacter pangenome revealed that it is notably open with a volatile accessory genome compared to previous pangenome studies on other genera, suggesting a high potential for adaptability to environmental niches.

A donor-DNA-free CRISPR/Cas-based approach to gene knock-up in rice.

Structural variations (SVs), such as inversion and duplication, contribute to important agronomic traits in crops1. Pan-genome studies revealed that SVs were a crucial and ubiquitous force driving genetic diversification2-4. Although genome editing can effectively create SVs in plants and animals5-8, the potential of designed SVs in breeding has been overlooked. Here, we show that new genes and traits can be created in rice by designed large-scale genomic inversion or duplication using CRISPR/Cas9. A 911 kb inversion on chromosome 1 resulted in a designed promoter swap between CP12 and PPO1, and a 338 kb duplication between HPPD and Ubiquitin2 on chromosome 2 created a novel gene cassette at the joint, promoterUbiquitin2::HPPD. Since the original CP12 and Ubiquitin2 genes were highly expressed in leaves, the expression of PPO1 and HPPD in edited plants with homozygous SV alleles was increased by tens of folds and conferred sufficient herbicide resistance in field trials without adverse effects on other important agronomic traits. CRISPR/Cas-based genome editing for gene knock-ups has been generally considered very difficult without inserting donor DNA as regulatory elements. Our study challenges this notion by providing a donor-DNA-free strategy, thus greatly expanding the utility of CRISPR/Cas in plant and animal improvements.

Towards a better understanding of the bacterial pan-genome

The bacterial pan-genome is a relatively new concept that refers to the number of genes observed in a given set of bacterial genome sequences, either at the intra- or inter-species level. Determining the pan-genome of a given species of bacteria using a large number of strains allows one to compare multiple genes and to determine evolutionary links between isolates. This information can help to determine population structure, diversity in terms of prevalence in a given environment and pathogenicity of microorganisms. Within this review, we explain the most important issues related to pan-genome studies. We also include a brief description of some selected bacterial pan-genomes. Finally, we propose an easy-toperform workflow to study bacterial pan-genomes that will facilitate nonexperts in a pan-genome-based investigation.

The pangenome of banana highlights differences between genera and genomes.

Banana (Musaceae family) has a complex genetic history and includes a genus Musa with a variety of cultivated clones with edible fruits, Ensete species that are grown for their edible corm, and monospecific Musella whose generic status has been questioned. The most commonly exported banana cultivars belong to Cavendish, a subgroup of Musa triploid cultivars, which is under threat by fungal pathogens, though there are also related species M. balbisiana Colla (B genome), M. textilis Née (T genome), and M. schizocarpa N. W. Simmonds (S genome), along with hybrids of these genomes, which potentially host genes of agronomic interest. Here we present the first cross-genus pangenome of banana, which contains representatives of the Musa and Ensete genera. Clusters based on gene presence-absence variation (PAV) clearly separate Musa and Ensete, while Musa is split further based on species. These results present the first pangenome study across genus boundaries and identifies genes that differentiate between Musaceae species, information that may support breeding programs in these crops.

Genomic region associated with pod color variation in pea (Pisum sativum).

Pea (Pisum sativum) was chosen as the research material by Gregor Mendel to discover the laws of inheritance. Out of seven traits studied by Mendel, genes controlling three traits including pod shape, pod color, and flower position have not been identified to date. With the aim of identifying the genomic region controlling pod color, we determined the genome sequence of a pea line with yellow pods. Genome sequence reads obtained using a Nanopore sequencing technology were assembled into 117,981 contigs (3.3 Gb), with an N50 value of 51.2 kb. A total of 531,242 potential protein-coding genes were predicted, of which 519,349 (2.8 Gb) were located within repetitive sequences (2.8 Gb). The assembled sequences were ordered using a reference as a guide to build pseudomolecules. Subsequent genetic and association analyses led to the identification of a genomic region that controls pea pod color. DNA sequences at this genomic location and transcriptome profiles of green and yellow pod lines were analyzed, and genes encoding 3' exoribonucleases were selected as potential candidates controlling pod color. The results presented in this study are expected to accelerate pan-genome studies in pea and facilitate the identification of the gene controlling one of the traits studied by Mendel.

Long-read sequence assembly: a technical evaluation in barley.

Sequence assembly of large and repeat-rich plant genomes has been challenging, requiring substantial computational resources and often several complementary sequence assembly and genome mapping approaches. The recent development of fast and accurate long-read sequencing by circular consensus sequencing (CCS) on the PacBio platform may greatly increase the scope of plant pan-genome projects. Here, we compare current long-read sequencing platforms regarding their ability to rapidly generate contiguous sequence assemblies in pan-genome studies of barley (Hordeum vulgare). Most long-read assemblies are clearly superior to the current barley reference sequence based on short-reads. Assemblies derived from accurate long reads excel in most metrics, but the CCS approach was the most cost-effective strategy for assembling tens of barley genomes. A downsampling analysis indicated that 20-fold CCS coverage can yield very good sequence assemblies, while even five-fold CCS data may capture the complete sequence of most genes. We present an updated reference genome assembly for barley with near-complete representation of the repeat-rich intergenic space. Long-read assembly can underpin the construction of accurate and complete sequences of multiple genomes of a species to build pan-genome infrastructures in Triticeae crops and their wild relatives.