AFLAP: assembly-free linkage analysis pipeline using k-mers from genome sequencing data

A k-mer-based bulked segregant analysis approach to map seed traits in unphased heterozygous potato genomes

High-throughput sequencing-based methods for bulked segregant analysis (BSA) allow for the rapid identification of genetic markers associated with traits of interest. BSA studies have successfully identified qualitative (binary) and quantitative trait loci (QTLs) using QTL mapping. However, most require population structures that fit the models available and a reference genome. Instead, high-throughput short-read sequencing can be combined with BSA of k-mers (BSA-k-mer) to map traits that appear refractory to standard approaches. This method can be applied to any organism and is particularly useful for species with genomes diverged from the closest sequenced genome. It is also instrumental when dealing with highly heterozygous and potentially polyploid genomes without phased haplotype assemblies and for which a single haplotype can control a trait. Finally, it is flexible in terms of population structure. Here, we apply the BSA-k-mer method for the rapid identification of candidate regions related to seed spot and seed size in diploid potato. Using a mixture of F1 and F2 individuals from a cross between 2 highly heterozygous parents, candidate sequences were identified for each trait using the BSA-k-mer approach. Using parental reads, we were able to determine the parental origin of the loci. Finally, we mapped the identified k-mers to a closely related potato genome to validate the method and determine the genomic loci underlying these sequences. The location identified for the seed spot matches with previously identified loci associated with pigmentation in potato. The loci associated with seed size are novel. Both loci are relevant in future breeding toward true seeds in potato.

Genome-Enabled Insights into Downy Mildew Biology and Evolution

Oomycetes that cause downy mildew diseases are highly specialized, obligately biotrophic phytopathogens that can have major impacts on agriculture and natural ecosystems. Deciphering the genome sequence of these organisms provides foundational tools to study and deploy control strategies against downy mildew pathogens (DMPs). The recent telomere-to-telomere genome assembly of the DMP Peronospora effusa revealed high levels of synteny with distantly related DMPs, higher than expected repeat content, and previously undescribed architectures. This provides a road map for generating similar high-quality genome assemblies for other oomycetes. This review discusses biological insights made using this and other assemblies, including ancestral chromosome architecture, modes of sexual and asexual variation, the occurrence of heterokaryosis, candidate gene identification, functional validation, and population dynamics. We also discuss future avenues of research likely to be fruitful in studies of DMPs and highlight resources necessary for advancing our understanding and ability to forecast and control disease outbreaks.

The genome of the oomycete Peronosclerospora sorghi, a cosmopolitan pathogen of maize and sorghum, is inflated with dispersed pseudogenes

Several species in the oomycete genus Peronosclerospora cause downy mildew on maize and can result in significant yield losses in Asia. Bio-surveillance of these pathogens is a high priority to prevent epidemics on maize in the United States and consequent damage to the US economy. The unresolved taxonomy and dearth of molecular resources for Peronosclerospora spp. hinder these efforts. P. sorghi is a pathogen of sorghum and maize with a global distribution, for which limited diversity has been detected in the southern USA. We characterized the genome, transcriptome, and mitogenome of an isolate, representing the US pathotype 6. The highly homozygous genome was assembled using 10× Genomics linked reads and scaffolded using Hi-C into 13 chromosomes. The total assembled length was 303.2 Mb, larger than any other oomycete previously assembled. The mitogenome was 38 kb, similar in size to other oomycetes, although it had a unique gene order. Nearly 20,000 genes were annotated in the nuclear genome, more than described for other downy mildew causing oomycetes. The 13 chromosomes of P. sorghi were highly syntenic with the 17 chromosomes of Peronospora effusa with conserved centromeric regions and distinct chromosomal fusions. The increased assembly size and gene count of P. sorghi is due to extensive retrotransposition, resulting in putative pseudogenization. Ancestral genes had higher transcript abundance and were enriched for differential expression. This study provides foundational resources for analysis of Peronosclerospora and comparisons to other oomycete genera. Further genomic studies of global Peronosclerospora spp. will determine the suitability of the mitogenome, ancestral genes, and putative pseudogenes for marker development and taxonomic relationships.

Ancestral Chromosomes for Family Peronosporaceae Inferred from a Telomere-to-Telomere Genome Assembly of Peronospora effusa

Downy mildew disease of spinach, caused by the oomycete Peronospora effusa, causes major losses to spinach production. In this study, the 17 chromosomes of P. effusa were assembled telomere-to-telomere, using Pacific Biosciences high-fidelity reads. Of these, 16 chromosomes are complete and gapless; chromosome 15 contains one gap bridging the nucleolus organizer region. This is the first telomere-to-telomere genome assembly for an oomycete. Putative centromeric regions were identified on all chromosomes. This new assembly enables a reevaluation of the genomic composition of Peronospora spp.; the assembly was almost double the size and contained more repeat sequences than previously reported for any Peronospora species. Genome fragments consistently underrepresented in six previously reported assemblies of P. effusa typically encoded repeats. Some genes annotated as encoding effectors were organized into multigene clusters on several chromosomes. Putative effectors were annotated on 16 of the 17 chromosomes. The intergenic distances between annotated genes were consistent with compartmentalization of the genome into gene-dense and gene-sparse regions. Genes encoding putative effectors were enriched in gene-sparse regions. The near-gapless assembly revealed apparent horizontal gene transfer from Ascomycete fungi. Gene order was highly conserved between P. effusa and the genetically oriented assembly of the oomycete Bremia lactucae; high levels of synteny were also detected with Phytophthora sojae. Extensive synteny between phylogenetically distant species suggests that many other oomycete species may have similar chromosome organization. Therefore, this assembly provides the foundation for genomic analyses of diverse oomycetes.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

HGGA: hierarchical guided genome assembler

BACKGROUND

De novo genome assembly typically produces a set of contigs instead of the complete genome. Thus additional data such as genetic linkage maps, optical maps, or Hi-C data is needed to resolve the complete structure of the genome. Most of the previous work uses the additional data to order and orient contigs.

RESULTS

Here we introduce a framework to guide genome assembly with additional data. Our approach is based on clustering the reads, such that each read in each cluster originates from nearby positions in the genome according to the additional data. These sets are then assembled independently and the resulting contigs are further assembled in a hierarchical manner. We implemented our approach for genetic linkage maps in a tool called HGGA.

CONCLUSIONS

Our experiments on simulated and real Pacific Biosciences long reads and genetic linkage maps show that HGGA produces a more contiguous assembly with less contigs and from 1.2 to 9.8 times higher NGA50 or N50 than a plain assembly of the reads and 1.03 to 6.5 times higher NGA50 or N50 than a previous approach integrating genetic linkage maps with contig assembly. Furthermore, also the correctness of the assembly remains similar or improves as compared to an assembly using only the read data.

High-resolution genetic dissection of the major QTL for tipburn resistance in lettuce, Lactuca sativa

Tipburn is an important physiological disorder of lettuce, Lactuca sativa L., related to calcium deficiency that can result in leaf necrosis and unmarketable crops. The major quantitative trait locus (QTL), qTPB5.2, can account for up to 70% of the phenotypic variance for tipburn incidence in the field. This QTL was genetically dissected to identify candidate genes for tipburn by creating lines with recombination events within the QTL and assessing their resistance to tipburn. By comparing lines with contrasting haplotypes, the genetic region was narrowed down to ∼877 Kb that was associated with a reduction of tipburn by ∼60%. Analysis of the lettuce reference genome sequence revealed 12 genes in this region, one of which is a calcium transporter with a single nucleotide polymorphism in an exon between haplotypes with contrasting phenotypes. RNA-seq analysis of recombinants revealed two genes that were differentially expressed between contrasting haplotypes consistent with the tipburn phenotype. One encodes a Teosinte branched1/Cycloidea/Proliferating Cell factor transcription factor; however, differential expression of the calcium transporter was detected. The phenotypic data indicated that there is a second region outside of the ∼877 Kb region but within the QTL, at which a haplotype from the susceptible parent decreased tipburn by 10-20%. A recombinant line was identified with beneficial haplotypes in each region from both parents that showed greater tipburn resistance than the resistant parent; this line could be used as the foundation for breeding cultivars with more resistance than is currently available.