Complete Chromosome Research Articles

De novo Phased Genome Assembly, Annotation and Population Genotyping of Alectoris Chukar

The Alectoris Chukar (chukar) is the most geographically widespread partridge species in the world, demonstrating exceptional adaptability to diverse ecological environments. However, the scarcity of genetic resources for chukar has hindered research into its adaptive evolution and molecular breeding. In this study, we have sequenced and assembled a high-quality, phased chukar genome that consists of 31 pairs of relatively complete diploid chromosomes. Our BUSCO analysis reported a high completeness score of 96.8% and 96.5%, with respect to universal single-copy orthologs and a low duplication rate (0.3% and 0.5%) for two assemblies. Through resequencing and population genomic analyses of six subspecies, we have curated invaluable genotype data that underscores the adaptive evolution of chukar in response to both arid and high-altitude environments. These data will significantly contribute to research on how chukars adaptively evolve to cope with desertification and alpine climates.

The Genome of Plasmodium gonderi: Insights into the Evolution of Human Malaria Parasites.

Plasmodium species causing malaria in humans are not monophyletic, sharing common ancestors with nonhuman primate parasites. Plasmodium gonderi is one of the few known Plasmodium species infecting African old-world monkeys that are not found in apes. This study reports a de novo assembled P. gonderi genome with complete chromosomes. The P. gonderi genome shares codon usage, syntenic blocks, and other characteristics with the human parasites Plasmodium ovale s.l. and Plasmodium malariae, also of African origin, and the human parasite Plasmodium vivax and species found in nonhuman primates from Southeast Asia. Using phylogenetically aware methods, newly identified syntenic blocks were found enriched with conserved metabolic genes. Regions outside those blocks harbored genes encoding proteins involved in the vertebrate host-Plasmodium relationship undergoing faster evolution. Such genome architecture may have facilitated colonizing vertebrate hosts. Phylogenomic analyses estimated the common ancestor between P. vivax and an African ape parasite P. vivax-like, within the Asian nonhuman primates parasites clade. Time estimates incorporating P. gonderi placed the P. vivax and P. vivax-like common ancestor in the late Pleistocene, a time of active migration of hominids between Africa and Asia. Thus, phylogenomic and time-tree analyses are consistent with an Asian origin for P. vivax and an introduction of P. vivax-like into Africa. Unlike other studies, time estimates for the clade with Plasmodium falciparum, the most lethal human malaria parasite, coincide with their host species radiation, African hominids. Overall, the newly assembled genome presented here has the quality to support comparative genomic investigations in Plasmodium.

Natural variation in the zinc-finger-encoding exon of Prdm9 affects hybrid sterility phenotypes in mice.

PRDM9-mediated reproductive isolation was first described in the progeny of Mus musculus musculus (MUS) PWD/Ph and Mus musculus domesticus (DOM) C57BL/6J inbred strains. These male F1 hybrids fail to complete chromosome synapsis and arrest meiosis at prophase I, due to incompatibilities between the Prdm9 gene and hybrid sterility locus Hstx2. We identified 14 alleles of Prdm9 in exon 12, encoding the DNA-binding domain of the PRDM9 protein in outcrossed wild mouse populations from Europe, Asia, and the Middle East, 8 of which are novel. The same allele was found in all mice bearing introgressed t-haplotypes encompassing Prdm9. We asked whether 7 novel Prdm9 alleles in MUS populations and the t-haplotype allele in 1 MUS and 3 DOM populations induce Prdm9-mediated reproductive isolation. The results show that only combinations of the dom2 allele of DOM origin and the MUS msc1 allele ensure complete infertility of intersubspecific hybrids in outcrossed wild populations and inbred mouse strains examined so far. The results further indicate that MUS mice may share the erasure of PRDM9msc1 binding motifs in populations with different Prdm9 alleles, which implies that erased PRDM9 binding motifs may be uncoupled from their corresponding Prdm9 alleles at the population level. Our data corroborate the model of Prdm9-mediated hybrid sterility beyond inbred strains of mice and suggest that sterility alleles of Prdm9 may be rare.

Tracking Karyotype Changes in Treatment-Induced Drug-Resistant Evolution.

Karyotype coding, which encompasses the complete chromosome sets and their topological genomic relationships within a given species, encodes system-level information that organizes and preserves genes' function, and determines the macroevolution of cancer. This new recognition emphasizes the crucial role of karyotype characterization in cancer research. To advance this cancer cytogenetic/cytogenomic concept and its platforms, this study outlines protocols for monitoring the karyotype landscape during treatment-induced rapid drug resistance in cancer. It emphasizes four key perspectives: combinational analyses of phenotype and karyotype, a focus on the entire evolutionary process through longitudinal analysis, a comparison of whole landscape dynamics by including various types of NCCAs (including genome chaos), and the use of the same process to prioritize different genomic scales. This protocol holds promise for studying numerous evolutionary aspects of cancers, and it further enhances the power of karyotype analysis in cancer research.

Telomere-to-Telomere Genome Assembly of Tibetan Medicinal Mushroom Ganoderma leucocontextum and the First Copia Centromeric Retrotransposon in Macro-Fungi Genome.

A complete telomere-to-telomere (T2T) genome has been a longstanding goal in the field of genomic research. By integrating high-coverage and precise long-read sequencing data using multiple assembly strategies, we present here the first T2T gap-free genome assembly of Ganoderma leucocontextum strain GL72, a Tibetan medicinal mushroom. The T2T genome, with a size of 46.69 Mb, consists 13 complete nuclear chromosomes and typical telomeric repeats (CCCTAA)n were detected at both ends of 13 chromosomes. The high mapping rate, uniform genome coverage, a complete BUSCOs of 99.7%, and base accuracy exceeding 99.999% indicate that this assembly represents the highest level of completeness and quality. Regions characterized by distinct structural attributes, including highest Hi-C interaction intensity, high repeat content, decreased gene density, low GC content, and minimal or no transcription levels across all chromosomes may represent potential centromeres. Sequence analysis revealed the first Copia centromeric retrotransposon in macro-fungi genome. Phylogenomic analysis identified that G. leucocontextum and G. tsugae diverged from the other Ganoderma species approximately 9.8-17.9 MYA. The prediction of secondary metabolic clusters confirmed the capability of this fungus to produce a substantial quantity of metabolites. This T2T gap-free genome will contribute to the genomic 'dark matter' elucidation and server as a great reference for genetics, genomics, and evolutionary studies of G. leucocontextum.

Development of a novel streamlined workflow (AACRE) and database (inCREDBle) for genomic analysis of carbapenem-resistant Enterobacterales.

In response to the threat of increasing antimicrobial resistance, we must increase the amount of available high-quality genomic data gathered on antibiotic-resistant bacteria. To this end, we developed an integrated pipeline for high-throughput long-read sequencing, assembly, annotation and analysis of bacterial isolates and used it to generate a large genomic data set of carbapenemase-producing Enterobacterales (CPE) isolates collected in Spain. The set of 461 isolates were sequenced with a combination of both Illumina and Oxford Nanopore Technologies (ONT) DNA sequencing technologies in order to provide genomic context for chromosomal loci and, most importantly, structural resolution of plasmids, important determinants for transmission of antimicrobial resistance. We developed an informatics pipeline called Assembly and Annotation of Carbapenem-Resistant Enterobacteriaceae (AACRE) for the full assembly and annotation of the bacterial genomes and their complement of plasmids. To explore the resulting genomic data set, we developed a new database called inCREDBle that not only stores the genomic data, but provides unique ways to filter and compare data, enabling comparative genomic analyses at the level of chromosomes, plasmids and individual genes. We identified a new sequence type, ST5000, and discovered a genomic locus unique to ST15 that may be linked to its increased spread in the population. In addition to our major objective of generating a large regional data set, we took the opportunity to compare the effects of sample quality and sequencing methods, including R9 versus R10 nanopore chemistry, on genome assembly and annotation quality. We conclude that converting short-read and hybrid microbial sequencing and assembly workflows to the latest nanopore chemistry will further reduce processing time and cost, truly enabling the routine monitoring of resistance transmission patterns at the resolution of complete chromosomes and plasmids.

A chromosome-level genome assembly of Zasmidium syzygii isolated from banana leaves.

Accurate taxonomic classification of samples from infected host material is essential for disease diagnostics and genome analyses. Despite the importance, diagnosis of fungal pathogens causing banana leaf diseases remains challenging. Foliar diseases of bananas are mainly caused by 3 Pseudocercospora species, of which the most predominant causal agent is Pseudocercospora fijiensis. Here, we sequenced and assembled four fungal isolates obtained from necrotic banana leaves in Bohol (Philippines) and obtained a high-quality genome assembly for one of these isolates. The samples were initially identified as P. fijiensis using PCR diagnostics; however, the assembly size was consistently 30 Mb smaller than expected. Based on the internal transcribed spacer (ITS) sequences, we identified the samples as Zasmidium syzygii (98.7% identity). The high-quality Zasmidium syzygii assembly is 42.5 Mb in size, comprising 16 contigs, of which 11 are most likely complete chromosomes. The genome contains 98.6% of the expected single-copy BUSCO genes and contains 14,789 genes and 10.3% repeats. The 3 short-read assemblies are less continuous but have similar genome sizes (40.4-42.4 Mb) and contain between 96.5 and 98.4% BUSCO genes. All 4 isolates have identical ITS sequences and are distinct from Zasmidium isolates that were previously sampled from banana leaves. We thus report the first continuous genome assembly of a member of the Zasmidium genus, forming an essential resource for further analysis to enhance our understanding of the diversity of pathogenic fungal isolates as well as fungal diversity.

Tigecycline-resistant Klebsiella pneumoniae strains from sewage in Norway carry heavy-metal resistance genes encoding conjugative plasmids

Tigecycline is a last-resort antibiotic used for treatment of infections with carbapenem-resistant Klebsiella pneumoniae. The aim of the study was to understand the genetic mechanism of resistance and the genetic context of resistance genes in two tigecycline-resistant K. pneumoniae strains isolated from sewage in Bergen, Norway. Complete genome sequencing of the two strains was accomplished using a combination of short-read Illumina MiSeq-based and long-read Oxford Nanopore MinION-based sequencing. Conjugation experiments were performed using filter mating and a green fluorescent protein (GFP)-tagged Escherichia coli strain. The complete genome sequences of strain K6-320.1 and strain K7-325 were assembled into two contigs for each strain, one contig representing the complete circular chromosomes of 5 223 440 bp (K6-320.1) and 5 263 092 bp (K7-325), respectively, and the other representing plasmids with sizes of 276 509 bp (pK6-320.1) and 246 731 bp (pK7-325). Strain K6-320.1 belonged to sequence type (ST)869, whereas strain K7-325 belonged to the pathogenic ST307. Both plasmids belonged to the IncFIB(K)/IncFII(K) group and carried several antibiotic resistance genes (ARGs), including tet(A) and blaCTX-M. Both plasmids (pK6-320.1 and pK7-325) were transferred to a GFP-tagged E. coli strain, leading to the acquisition of resistance against multiple classes of antibiotics. Several heavy-metal resistance genes (HMRGs) encoding resistance against silver (sil), copper (pco), and arsenic (ars) were also present on both plasmids. Our study demonstrates the presence of multidrug-resistant K. pneumoniae strains carrying conjugative plasmids encoding both ARGs and HMRGs that have potential for persistence in the environment and human microbiota.

Did the creeping vole sex chromosomes evolve through a cascade of adaptive responses to a selfish x chromosome?

The creeping vole Microtus oregoni exhibits remarkably transformed sex chromosome biology, with complete chromosome drive/drag, X-Y fusions, sex reversed X complements, biased X inactivation, and X chromosome degradation. Beginning with a selfish X chromosome, I propose a series of adaptations leading to this system, each compensating for deleterious consequences of the preceding adaptation: (1) YY embryonic inviability favored evolution of a selfish feminizing X chromosome; (2) the consequent Y chromosome transmission disadvantage favored X-Y fusion ("XP "); (3) Xist-based silencing of Y-derived XP genes favored a second X-Y fusion ("XM "); (4) X chromosome dosage-related costs in XP XM males favored the evolution of XM loss during spermatogenesis; (5) X chromosomal dosage-related costs in XM 0 females favored the evolution of XM drive during oogenesis; and (6) degradation of the non-recombining XP favored the evolution of biased X chromosome inactivation. I discuss recurrent rodent sex chromosome transformation, and selfish genes as a constructive force in evolution.

Automated Cytogenetic Classifier (ACC) to Support Accurate Clinical Decision Making in Acute Myeloid Leukaemia

Introduction Cytogenetic analysis is a recommended test in newly diagnosed AML but requires accurate interpretation of the International System for Human Cytogenomic Nomenclature (ISCN). Keeping up to date across recent published classifications - WHO, the European Leukaemia Net (ELN) and International Consensus Classification (ICC) - is critical for accurate prognostic assessment, but can be challenging for clinicians and laboratory scientists. Furthermore, eligibility to access state-funded treatments in the UK is dictated by certain cytogenetic results. We previously developed a webapp that integrates cytogenetic, molecular and clinical features to aid accurate classification, prognostication and provide treatment recommendations in AML (Coats et al, BJH 2021). To further improve the webapp we have developed a cytogenetic classifier to automatically interpret reports written in ISCN format. In this study we aimed to validate the automated cytogenetic classifier (ACC) in identifying abnormalities that feature in the new WHO, ELN and ICC classifications. Methodology Software for the ACC was written in Python to: 1. count the number and type of karyotypic abnormalities (subdivided into number of monosomies, complete chromosome gains and structural changes) 2. identify presence of clinically relevant, named, abnormalities that feature in international classifications (categorised as an individual abnormalityg. t(9;11)(p21.3;q23.3), or a generic abnormality e.g. any translocation involving 11q23.3). Cytogenetic reports of 50 consecutive AML cases with abnormal karyotypes analysed by Bristol Haematological Oncology Diagnostic Service (BHODS) were selected (training cohort). Each report was uploaded to the ACC and its output independently checked for accuracy. Following review of the software's performance, adjustments to the ACC were made. The updated ACC (Figure 1) analysed 98 additional reports from two other laboratories, UCLH and OUH (validation cohort). Results The ACC analysed 49/50 cytogenetic reports from the BHODS laboratory. One report could not be read by the ACC. The correct quantity of abnormalities was identified in 49/49 cases analysed. Where relevant named abnormalities occurred, these were detected in 79/80 instances. One case of acute promyelocytic leukaemia containing an idic(17p) was not identified as containing an abnormal 17p. Additionally, in one report with a complex karyotype, a monosomy 13 was erroneously detected by the ACC. Adjustments were made to the ACC such that it could analyse all 50 cases with 100% accuracy. The updated ACC then analysed 98 cases from the validation cohort. The correct quantity of abnormalities was detected in 97/ 98 cases. In one case with a highly complex karyotype containing a tetraploid clone, the ACC incorrectly counted 48, as opposed to 45, abnormalities. Clinically relevant, named, abnormalities were correctly identified in 178/179 instances. In one case, containing a 4-way translocation, the ACC did not identify the t(9;11) abnormality. Discussion Here, we show that the updated ACC is accurate in identifying the number and type of karyotypic abnormalities relevant to AML from different laboratories. Reduced accuracy in counting the number of abnormalities in highly complex karyotypes (i.e. >40 abnormalities) is not likely to impact clinical interpretation. Further validation work is needed to analyse 3 and 4-way translocations. User warnings have been added to the ACC for karyotypes containing multiple-translocations and tetraploid/triploid clones to improve clinical safety. The ACC can support clinical decision-making by streamlining cytogenetic interpretation from across multiple clinical classifications. It can be used as a stand-alone programme or integrated into more sophisticated decision support tools. The ACC has the potential to be updated to interpret cytogenetic reports from other classifications or haematological malignancies. The ACC is available as a webapp (Streamlit (ts482-cytogenetics-calculator-extractorstreamlit-w7dq9p.streamlit.app). Figure 1. Flowchart of the methodology of the updated automated cytogenetic classifier (ACC).

Complete trisomy 9 detected by noninvasive prenatal testing and confirmed by amniocentesis.

Complete chromosome 9 trisomy (T9) is a rare and fatal chromosomal disorder. We performed non-invasive prenatal testing (NIPT) in a patient with threatened abortion symptoms and found that the fetal was at risk for complete chromosome 9 trisomy. This shows that NIPT has certain accuracy in detecting trisomy of chromosome 9, which provide options for prenatal diagnosis of rare chromosomal abnormalities.

Telomere-to-telomere genome of the allotetraploid legume Sesbania cannabina reveals transposon-driven subgenome divergence and mechanisms of alkaline stress tolerance.

Alkaline soils pose an increasing problem for agriculture worldwide, but using stress-tolerant plants as green manure can improve marginal land. Here, we show that the legume Sesbania cannabina is very tolerant to alkaline conditions and, when used as a green manure, substantially improves alkaline soil. To understand genome evolution and the mechanisms of stress tolerance in this allotetraploid legume, we generated the first telomere-to-telomere genome assembly of S. cannabina spanning ∼2,087 Mb. The assembly included all centromeric regions, which contain centromeric satellite repeats, and complete chromosome ends with telomeric characteristics. Further genome analysis distinguished A and B subgenomes, which diverged approximately 7.9 million years ago. Comparative genomic analysis revealed that the chromosome homoeologs underwent large-scale inversion events (>10 Mb) and a significant, transposon-driven size expansion of the chromosome 5A homoeolog. We further identified four specific alkali-induced phosphate transporter genes in S. cannabina; these may function in alkali tolerance by relieving the deficiency in available phosphorus in alkaline soil. Our work highlights the significance of S. cannabina as a green tool to improve marginal lands and sheds light on subgenome evolution and adaptation to alkaline soils.

Improving the Completeness of Chromosome-Level Assembly by Recalling Sequences from Lost Contigs.

For a long time, the construction of complete reference genomes for complex eukaryotic genomes has been hindered by the limitations of sequencing technologies. Recently, the Pacific Biosciences (PacBio) HiFi data and Oxford Nanopore Technologies (ONT) Ultra-Long data, leveraging their respective advantages in accuracy and length, have provided an opportunity for generating complete chromosome sequences. Nevertheless, for the majority of genomes, the chromosome-level assemblies generated using existing methods still miss a high proportion of sequences due to losing small contigs in the step of assembly and scaffolding. To address this shortcoming, in this paper, we propose a novel method that is able to identify and fill the gaps in the chromosome-level assembly by recalling the sequences in the lost small contigs. Experimental results on both real and simulated datasets demonstrate that this method is able to improve the completeness of the chromosome-level assembly.

Screening and Analysis of Antifungal Strains Bacillus subtilis JF-4 and B. amylum JF-5 for the Biological Control of Fusarium Wilt of Banana.

This study aimed to identify the antagonistic bacteria from the rhizosphere of healthy bananas that can effectively suppress the Fusarium wilt of banana, and to further investigate the inhibitory mechanism. The primary and secondary screening techniques were implemented using the double-plate and fermentation antagonism methods. The strain was identified based on physiological and biochemical tests, 16S rRNA gene sequencing, and specific gene amplification. The effects of crude extract on the protein content, lipid peroxidation, and pectinase activity of mycelia were determined from the identified isolates. Two antagonistic bacteria, JF-4 and JF-5, were screened and initially identified as Bacillus subtilis (GenBank: OR125631) and B. amylum (GenBank: OR125632). The greenhouse experiment showed that the biological control efficiency of the two antagonists against the Fusarium wilt of banana was 48.3% and 40.3%, respectively. The catalase content produced by lipid peroxidation increased significantly after treatment with the crude extracts of JF-4 and JF-5 at concentrations of 0.69 μmol/L and 0.59 μmol/L, respectively. The protein and ergosterol content and pectinase activity decreased significantly. The two antagonistic bacteria might inhibit the growth of pathogens by enhancing lipid peroxidation and decreasing the synthesis of cell metabolites. Twenty compounds were identified by gas chromatography-mass spectrometry (GC-MS). B. subtilis JF-4 was further sequenced and assembled to obtain a complete circular chromosome genome of 681,804,824 bp. The genome consisted of a 4,310,825-bp-long scaffold. The findings of this study may help elucidate the mechanism behind this biocontrol isolate.

Prenatal diagnosis of paternal uniparental disomy for chromosome 2 in two fetuses with intrauterine growth restriction

Uniparental disomy (UPD) is when all or part of the homologous chromosomes are inherited from only one of the two parents. Currently, UPD has been reported to occur for almost all chromosomes. In this study, we report two cases of UPD for chromosome 2 (UPD2) encountered during prenatal diagnosis. The ultrasound findings of the fetuses from two unrelated families showed intrauterine growth restriction. The karyotype analyses were normal. The two fetuses both had complete paternal chromosome 2 uniparental disomy detected by whole-exome sequencing, but their clinical outcomes were significantly different, with fetal arrest in case 1 and birth in case 2. In this report, we analyzed and discussed the phenotypes of the fetuses in these two cases and reviewed the literature on UPD2.

Near chromosome-level genome assembly of the microsporidium Hamiltosporidium tvaerminnensis.

Microsporidia are intracellular parasitic fungi whose genomes rank among the smallest of all known eukaryotes. A number of outstanding questions remain concerning the evolution of their large-scale variation in genome architecture, responsible for genome size variation of more than an order of magnitude. This genome report presents the first near-chromosomal assembly of a large-genome microsporidium, Hamiltosporidium tvaerminnensis. Combined Oxford Nanopore, Pacific Biosciences (PacBio), and Illumina sequencing led to a genome assembly of 17 contigs, 11 of which represent complete chromosomes. Our assembly is 21.64 Mb in length, has an N50 of 1.44 Mb, and consists of 39.56% interspersed repeats. We introduce a novel approach in microsporidia, PacBio Iso-Seq, as part of a larger annotation pipeline for obtaining high-quality annotations of 3,573 protein-coding genes. Based on direct evidence from the full-length Iso-Seq transcripts, we present evidence for alternative polyadenylation and variation in splicing efficiency, which are potential regulation mechanisms for gene expression in microsporidia. The generated high-quality genome assembly is a necessary resource for comparative genomics that will help elucidate the evolution of genome architecture in response to intracellular parasitism.

Homoeologous Chromosome Pairing and Alien Introgression in Backcrossing Progenies Derived from Hybrids Solanum tuberosum (+) Mexican 2x (1 EBN) B-Genome Potato Species

We applied a genomic in situ hybridization (GISH) to analyze the genomic constitution of and meiotic pairing in interspecific somatic hybrids, and in a wide subset of backcrossing derivatives (BC1–BC5), from three interspecific combinations involving the cultivated potato, Solanum tuberosum (AAAA genome), and three diploid (1 EBN) wild Mexican potato species (genome BB)—S. tarnii, S. pinnatisectum, and S. bulbocastanum. The theoretically expected genomic composition was detected in the somatic hybrids (AAAABB) and in the BC1 progeny (AAAAB), whereas in the subsequent BC2–BC4 generations, the partial loss of alien chromosomes was observed and almost all the BC5 genotypes showed a complete chromosome elimination of wild species. GISH revealed a homoeologous pairing between the chromosomes of the A- and the B-genomes in each of the hybrid progenies. Using GISH, we selected introgression lines with a single chromosome of the wild species in a potato genome background, as well as introgression lines with intergenomic recombinant chromosomes. Moreover, via molecular screening, BC hybrids with diagnostic markers for the R-genes conferring resistance to late blight disease and to the quarantine pest of the potato–Columbia root-knot nematode—were selected. The potential application of the results obtained for the planning of introgressive schemes directed to the breeding of advanced lines with multiple disease and pest resistance is discussed.

Plassembler: an automated bacterial plasmid assembly tool.

With recent advances in sequencing technologies, it is now possible to obtain near-perfect complete bacterial chromosome assemblies cheaply and efficiently by combining a long-read-first assembly approach with short-read polishing. However, existing methods for assembling bacterial plasmids from long-read-first assemblies often misassemble or even miss bacterial plasmids entirely and accordingly require manual curation. Plassembler was developed to provide a tool that automatically assembles and outputs bacterial plasmids using a hybrid assembly approach. It achieves increased accuracy and computational efficiency compared to the existing gold standard tool Unicycler by removing chromosomal reads from the input read sets using a mapping approach. Plassembler is implemented in Python and is installable as a bioconda package using 'conda install -c bioconda plassembler'. The source code is available on GitHub at https://github.com/gbouras13/plassembler. The full benchmarking pipeline can be found at https://github.com/gbouras13/plassembler_simulation_benchmarking, while the benchmarking input FASTQ and output files can be found at https://doi.org/10.5281/zenodo.7996690.

A chromosome-scale genome assembly of the grape powdery mildew pathogen Erysiphe necator reveals its genomic architecture and previously unknown features of its biology.

Erysiphe necator is an obligate fungal pathogen that causes grape powdery mildew, globally the most important disease on grapevines. Previous attempts to obtain a quality genome assembly for this pathogen were hindered by its high repetitive DNA content. Here, chromatin conformation capture (Hi-C) with long-read PacBio sequencing was combined to obtain a chromosome-scale assembly and a high-quality annotation for E. necator isolate EnFRAME01. The resulting 81.1 Mb genome assembly is 98% complete and consists of 34 scaffolds, 11 of which represent complete chromosomes. All chromosomes contain large centromeric-like regions and lack synteny to the 11 chromosomes of the cereal PM pathogen Blumeria graminis. Further analysis of their composition showed that repeats and transposable elements (TEs) occupy 62.7% of their content. TEs were almost evenly interspersed outside centromeric and telomeric regions and massively overlapped with regions of annotated genes, suggesting that they could have a significant functional impact. Abundant gene duplicates were observed as well, particularly in genes encoding candidate secreted effector proteins. Moreover, younger in age gene duplicates exhibited more relaxed selection pressure and were more likely to be located physically close in the genome than older duplicates. A total of 122 genes with copy number variations among six isolates of E. necator were also identified and were enriched in genes that were duplicated in EnFRAME01, indicating they may reflect an adaptive variation. Taken together, our study illuminates higher-order genomic architectural features of E. necator and provides a valuable resource for studying genomic structural variations in this pathogen. IMPORTANCE Grape powdery mildew caused by the ascomycete fungus Erysiphe necator is economically the most important and recurrent disease in vineyards across the world. The obligate biotrophic nature of E. necator hinders the use of typical genetic methods to elucidate its pathogenicity and adaptation to adverse conditions, and thus comparative genomics has been a major method to study its genome biology. However, the current reference genome of E. necator isolate C-strain is highly fragmented with many non-coding regions left unassembled. This incompleteness prohibits in-depth comparative genomic analyses and the study of genomic structural variations (SVs) that are known to affect several aspects of microbial life, including fitness, virulence, and host adaptation. By obtaining a chromosome-scale genome assembly and a high-quality gene annotation for E. necator, we reveal the organization of its chromosomal content, unearth previously unknown features of its biology, and provide a reference for studying genomic SVs in this pathogen.

Whole and nearly complete mitochondrial genomes of an endemic and endangered neotropical rabbit (Romerolagus diazi) assembled using non-invasive eDNA metagenomics (field droppings)

ABSTRACT The volcano rabbit Romerolagus diazi is an endemic neotropical species restricted to the central part of the Trans-Mexican Volcanic Belt and considered ‘endangered’ by the International Union for Conservation of Nature Red List of Threatened Species. This study formally tested whether complete mitochondrial genomes of R. diazi can be assembled from non-enriched metagenomic libraries generated from eDNA (field droppings) using Illumina 300 bp pair-end reads. Using a direct assembly strategy, the mitogenome of R. diazi was assembled and circularised using the pipeline GetOrganelle with a coverage of 18.5× from one eDNA sample. A second strategy in which reads were first mapped to a reference mitochondrial genome and then assembled using the software SPAdes retrieved a second nearly complete mitogenome (2.7×) of R. diazi from a second eDNA sample. Complete or nearly complete mitochondrial genomes were not retrieved from two other eDNA samples. The organisation of the complete or nearly complete mitochondrial genomes recovered from eDNA samples is identical to that already described for the studied species. A phylomitogenomic analysis based on all protein coding genes (PCGs) demonstrated that the mitochondrial genomes assembled from eDNA can and do reliably identify the sequenced samples as belonging to R. diazi and distinguished the same samples from closely and distantly related cofamilial species. This study demonstrates the feasibility of retrieving fully accurate whole and nearly complete mitochondrial chromosomes of host species from eDNA without enrichment strategies and using straightforward bioinformatics pipelines. This new genomic resource developed for R. diazi represents a new tool to boost conservation strategies (bioprospecting and/or biomonitoring) and improve our understanding of the evolutionary history of this endemic and endangered species.