Repetitive Genome Research Articles

Does time matter? Intraspecific diversity of ribosomal RNA genes in lineages of the allopolyploid model grass Brachypodium hybridum with different evolutionary ages

BackgroundPolyploidisation often results in genome rearrangements that may involve changes in both the single-copy sequences and the repetitive genome fraction. In this study, we performed a comprehensive comparative analysis of repetitive DNA, with a particular focus on ribosomal DNA (rDNA), in Brachypodium hybridum (2n = 4x = 30, subgenome composition DDSS), an allotetraploid resulting from a natural cross between two diploid species that resemble the modern B. distachyon (2n = 10; DD) and B. stacei (2n = 20; SS). Taking advantage of the recurrent origin of B. hybridum, we investigated two genotypes, Bhyb26 and ABR113, differing markedly in their evolutionary age (1.4 and 0.14 Mya, respectively) and which resulted from opposite cross directions. To identify the origin of rDNA loci we employed cytogenetic and molecular methods (FISH, gCAPS and Southern hybridisation), phylogenetic and genomic approaches.ResultsUnlike the general maintenance of doubled gene dosage in B. hybridum, the rRNA genes showed a remarkable tendency towards diploidisation at both locus and unit levels. While the partial elimination of 35S rDNA units occurred in the younger ABR113 lineage, unidirectional elimination of the entire locus was observed in the older Bhyb26 lineage. Additionally, a novel 5S rDNA family was amplified in Bhyb26 replacing the parental units. The 35S and 5S rDNA units were preferentially eliminated from the S- and D-subgenome, respectively. Thus, in the more ancient B. hybridum lineage, Bhyb26, 5S and 35S rRNA genes are likely expressed from different subgenomes, highlighting the complexity of polyploid regulatory networks.ConclusionComparative analyses between two B. hybridum lineages of distinct evolutionary ages revealed that although the recent lineage ABR113 exhibited an additive pattern of rDNA loci distribution, the ancient lineage Bhyb26 demonstrated a pronounced tendency toward diploidisation manifested by the reduction in the number of both 35S and 5S loci. In conclusion, the age of the allopolyploid appears to be a decisive factor in rDNA turnover in B. hybridum.

The Highly Repetitive Genome of Myxobolus rasmusseni, an Emerging Myxozoan Parasite of Fathead Minnows.

Myxozoans are a monophyletic taxon of approximately 2,400 described species of parasites from the phylum Cnidaria. The recent focus on their negative impacts on fisheries, on their evolution from free-living ancestors, and on their emergence into new fish host populations has stressed the critical need for genomic resources for this parasitic group. Here, we describe the genome assembly and annotation of Myxobolus rasmusseni, an emerging parasite of fathead minnows in Alberta, Canada. The assembly is 174.6 Mb in size, 68% of which is made up of repetitive elements, making it one of the most repetitive animal genomes sequenced to date. Through comparisons to other myxozoans, we show that widespread gene loss, a known phenomenon of this group of parasites, is consistent with closely related species. Additionally, we assembled the M. rasmusseni mitochondrial genome, which is nearly twice the size of the typical animal mitochondrial genome yet contains only five of the canonical mitochondrial protein-coding genes and open reading frames not found in other myxozoans. These results add to our understanding of the gene- and genome-level diversity observed in myxozoans.

Nanopore sequencing with unique molecular identifiers enables accurate mutation analysis and haplotyping in the complex lipoprotein(a) KIV-2 VNTR.

Repetitive genome regions, such as variable number of tandem repeats (VNTR) or short tandem repeats (STR), are major constituents of the uncharted dark genome and evade conventional sequencing approaches. The protein-coding LPA kringle IV type-2 (KIV-2) VNTR (5.6kb per unit, 1-40 units per allele) is a medically highly relevant example with a particularly intricate structure, multiple haplotypes, intragenic homologies, and an intra-VNTR STR. It is the primary regulator of plasma lipoprotein(a) [Lp(a)] concentrations, an important cardiovascular risk factor. Lp(a) concentrations vary widely between individuals and ancestries. Multiple variants and functional haplotypes in the LPA gene and especially in the KIV-2 VNTR strongly contribute to this variance. We evaluated the performance of amplicon-based nanopore sequencing with unique molecular identifiers (UMI-ONT-Seq) for SNP detection, haplotype mapping, VNTR unit consensus sequence generation, and copy number estimation via coverage-corrected haplotypes quantification in the KIV-2 VNTR. We used 15 human samples and low-level mixtures (0.5 to 5%) of KIV-2 plasmids as a validation set. We then applied UMI-ONT-Seq to extract KIV-2 VNTR haplotypes in 48 multi-ancestry 1000 Genome samples and analyzed at scale a poorly characterized STR within the KIV-2 VNTR. UMI-ONT-Seq detected KIV-2 SNPs down to 1% variant level with high sensitivity, specificity, and precision (0.977 ± 0.018; 1.000 ± 0.0005; 0.993 ± 0.02) and accurately retrieved the full-length haplotype of each VNTR unit. Human variant levels were highly correlated with next-generation sequencing (R2 = 0.983) without bias across the whole variant level range. Six reads per UMI produced sequences of each KIV-2 unit with Q40 quality. The KIV-2 repeat number determined by coverage-corrected unique haplotype counting was in close agreement with droplet digital PCR (ddPCR), with 70% of the samples falling even within the narrow confidence interval of ddPCR. We then analyzed 62,679 intra-KIV-2 STR sequences and explored KIV-2 SNP haplotype patterns across five ancestries. UMI-ONT-Seq accurately retrieves the SNP haplotype and precisely quantifies the VNTR copy number of each repeat unit of the complex KIV-2 VNTR region across multiple ancestries. This study utilizes the KIV-2 VNTR, presenting a novel and potent tool for comprehensive characterization of medically relevant complex genome regions at scale.

Highly Repetitive Genome of Coniella granati (syn. Pilidiella granati), the Causal Agent of Pomegranate Fruit Rot, Encodes a Minimalistic Proteome with a Streamlined Arsenal of Effector Proteins.

This study describes the first genome sequence and analysis of Coniella granati, a fungal pathogen with a broad host range, which is responsible for postharvest crown rot, shoot blight, and canker diseases in pomegranates. C. granati is a geographically widespread pathogen which has been reported across Europe, Asia, the Americas, and Africa. Our analysis revealed a 46.8 Mb genome with features characteristic of hemibiotrophic fungi. Approximately one third of its genome was compartmentalised within 'AT-rich' regions exhibiting a low GC content (30 to 45%). These regions primarily comprised transposable elements that are repeated at a high frequency and interspersed throughout the genome. Transcriptome-supported gene annotation of the C. granati genome revealed a streamlined proteome, mirroring similar observations in other pathogens with a latent phase. The genome encoded a relatively compact set of 9568 protein-coding genes with a remarkable 95% having assigned functional annotations. Despite this streamlined nature, a set of 40 cysteine-rich candidate secreted effector-like proteins (CSEPs) was predicted as well as a gene cluster involved in the synthesis of a pomegranate-associated toxin. These potential virulence factors were predominantly located near repeat-rich and AT-rich regions, suggesting that the pathogen evades host defences through Repeat-Induced Point mutation (RIP)-mediated pseudogenisation. Furthermore, 23 of these CSEPs exhibited homology to known effector and pathogenicity genes found in other hemibiotrophic pathogens. The study establishes a foundational resource for the study of the genetic makeup of C. granati, paving the way for future research on its pathogenicity mechanisms and the development of targeted control strategies to safeguard pomegranate production.

Haematococcus lacustris genome assembly and annotation reveal diploid genetic traits and stress-induced gene expression patterns

The green alga Haematococcus lacustris (formerly Haematococcus pluvialis) is a primary source of astaxanthin, a ketocarotenoid with high antioxidant activity and several industrial applications. Here, the Haematococcus lacustris highly repetitive genome was reconstructed by exploiting next-generation sequencing integrated with Hi-C scaffolding, obtaining a 151 Mb genome assembly in 32 scaffolds at a near-chromosome level with high continuity. Surprisingly, the distribution of the single-nucleotide-polymorphisms identified demonstrates a diploid configuration for the Haematococcus genome, further validated by Sanger sequencing of heterozygous regions. Functional annotation and RNA-seq data enabled the identification of 13,946 nuclear genes, with >5000 genes not previously identified in this species, providing insights into the molecular basis for metabolic rearrangement in stressing conditions such as high light and/or nitrogen starvation, where astaxanthin biosynthesis is triggered. These data constitute a rich genetic resource for biotechnological manipulation of Haematococcus lacustris highlighting potential targets to improve astaxanthin and carotenoid productivity.

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution.

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

Whole-genome assembly of a hybrid Trypanosoma cruzi strain assembled with Nanopore sequencing alone.

Trypanosoma cruzi is the causative agent of Chagas disease, which causes 10,000 deaths per year. Despite the high mortality associated with Chagas, relatively few parasite genomes have been assembled to date, with genome assemblies unavailable even for some commonly used laboratory strains. This is at least partially due to T. cruzi's highly complex and highly repetitive genome, which defies investigation using traditional short-read sequencing methods. In this study, we have generated a high-quality whole-genome assembly of the hybrid Tulahuen strain, a commercially available type VI strain, using long-read Nanopore sequencing without short-read scaffolding. The assembled genome contains 25% repeat regions, 17% variable multigene family members, and 27% transposable elements (TEs) and is of comparable quality with T. cruzi genome assemblies that utilized both long- and short-read data. Notably, we find that regions with TEs are significantly enriched for multicopy surface proteins, and that surface proteins are, on average, closer to TEs than to other coding regions. This finding suggests that mobile genetic elements such as transposons may drive recombination within surface protein gene families. This work demonstrates the feasibility of Nanopore sequencing to resolve complex regions of T. cruzi genomes, and with these resolved regions, provides support for a possible mechanism for genomic diversification.

A comprehensive benchmark of graph-based genetic variant genotyping algorithms on plant genomes for creating an accurate ensemble pipeline.

Although sequencing technologies have boosted the measurement of the genomic diversity of plant crops, it remains challenging to accurately genotype millions of genetic variants, especially structural variations, with only short reads. In recent years, many graph-based variation genotyping methods have been developed to address this issue and tested for human genomes. However, their performance in plant genomes remains largely elusive. Furthermore, pipelines integrating the advantages of current genotyping methods might be required, considering the different complexity of plant genomes. Here we comprehensively evaluate eight such genotypers in different scenarios in terms of variant type and size, sequencing parameters, genomic context, and complexity, as well as graph size, using both simulated and real data sets from representative plant genomes. Our evaluation reveals that there are still great challenges to applying existing methods to plants, such as excessive repeats and variants or high resource consumption. Therefore, we propose a pipeline called Ensemble Variant Genotyper (EVG) that can achieve better genotyping performance in almost all experimental scenarios and comparably higher genotyping recall and precision even using 5× reads. Furthermore, we demonstrate that EVG is more robust with an increasing number of graphed genomes, especially for insertions and deletions. Our study will provide new insights into the development and application of graph-based genotyping algorithms. We conclude that EVG provides an accurate, unbiased, and cost-effective way for genotyping both small and large variations and will be potentially used in population-scale genotyping for large, repetitive, and heterozygous plant genomes.

Cytogenomic Characterization of Transposable Elements and Satellite DNA in Passiflora L. Species.

The species Passiflora alata, P. cincinnata, and P. edulis have great economic value due to the use of their fruits for human consumption. In this study, we compared the repetitive genome fractions of these three species. The compositions of the repetitive DNA of these three species' genomes were analyzed using clustering and identification of the repetitive sequences with RepeatExplorer. It was found that repetitive DNA content represents 74.70%, 66.86%, and 62.24% of the genome of P. alata, P. edulis, and P. cincinnata, respectively. LTR Ty3/Gypsy retrotransposons represent the highest genome proportions in P. alata and P. edulis, while Ty1/Copia comprises the largest proportion of P. cincinnata genome. Chromosomal mapping by Fluorescent In Situ Hybridization (FISH) showed that LTR retrotransposons have a dispersed distribution along chromosomes. The subtelomeric region of chromosomes is where 145 bp satellite DNA is located, suggesting that these elements may play important roles in genome structure and organization in these species. In this work, we obtained the first global characterization of the composition of repetitive DNA in Passiflora, showing that an increase in genome size is related to an increase in repetitive DNA, which represents an important evolutionary route for these species.

Microevolution and Its Impact on Hypervirulence, Antimicrobial Resistance, and Vaccine Escape in Neisseria meningitidis.

Neisseria meningitidis is commensal of the human pharynx and occasionally invades the host, causing the life-threatening illness invasive meningococcal disease. The meningococcus is a highly diverse and adaptable organism thanks to natural competence, a propensity for recombination, and a highly repetitive genome. These mechanisms together result in a high level of antigenic variation to invade diverse human hosts and evade their innate and adaptive immune responses. This review explores the ways in which this diversity contributes to the evolutionary history and population structure of the meningococcus, with a particular focus on microevolution. It examines studies on meningococcal microevolution in the context of within-host evolution and persistent carriage; microevolution in the context of meningococcal outbreaks and epidemics; and the potential of microevolution to contribute to antimicrobial resistance and vaccine escape. A persistent theme is the idea that the process of microevolution contributes to the development of new hyperinvasive meningococcal variants. As such, microevolution in this species has significant potential to drive future public health threats in the form of hypervirulent, antibiotic-resistant, vaccine-escape variants. The implications of this on current vaccination strategies are explored.

Whole-genome resequencing facilitates the development of a 50K single nucleotide polymorphism genotyping array for Scots pine (Pinus sylvestris L.) and its transferability to other pine species.

Scots pine (Pinus sylvestris L.) is one of the most widespread and economically important conifer species in the world. Applications like genomic selection and association studies, which could help accelerate breeding cycles, are challenging in Scots pine because of its large and repetitive genome. For this reason, genotyping tools for conifer species, and in particular for Scots pine, are commonly based on transcribed regions of the genome. In this article, we present the Axiom Psyl50K array, the first single nucleotide polymorphism (SNP) genotyping array for Scots pine based on whole-genome resequencing, that represents both genic and intergenic regions. This array was designed following a two-step procedure: first, 192 trees were sequenced, and a 430K SNP screening array was constructed. Then, 480 samples, including haploid megagametophytes, full-sib family trios, breeding population, and range-wide individuals from across Eurasia were genotyped with the screening array. The best 50K SNPs were selected based on quality, replicability, distribution across the draft genome assembly, balance between genic and intergenic regions, and genotype-environment and genotype-phenotype associations. Of the final 49 877 probes tiled in the array, 20 372 (40.84%) occur inside gene models, while the rest lie in intergenic regions. We also show that the Psyl50K array can yield enough high-confidence SNPs for genetic studies in pine species from North America and Eurasia. This new genotyping tool will be a valuable resource for high-throughput fundamental and applied research of Scots pine and other pine species.

Orientia tsutsugamushi: comprehensive analysis of the mobilome of a highly fragmented and repetitive genome reveals the capacity for ongoing lateral gene transfer in an obligate intracellular bacterium.

Obligate intracellular bacteria, or those only capable of growth inside other living cells, have limited opportunities for horizontal gene transfer with other microbes due to their isolated replicative niche. The human pathogen Ot, an obligate intracellular bacterium causing scrub typhus, encodes an unusually high copy number of a ~40 gene mobile genetic element that typically facilitates genetic transfer across microbes. This proliferated element is heavily degraded in Ot and previously assumed to be inactive. Here, we conducted a detailed analysis of this element in eight Ot strains and discovered two strains with at least one intact copy. This implies that the element is still capable of moving across Ot populations and suggests that the genome of this bacterium may be even more dynamic than previously appreciated. Our work raises questions about intracellular microbial evolution and sounds an alarm for gene-based efforts focused on diagnosing and combatting scrub typhus.

The highly improved genome of Ixodes scapularis with X and Y pseudochromosomes.

Ixodes scapularis, the black-legged tick, is the principal vector of the Lyme disease spirochete, Borrelia burgdorferi, and is responsible for most of the ∼470,000 estimated Lyme disease cases annually in the USA. Ixodes scapularis can transmit six additional pathogens of human health significance. Because of its medical importance, I. scapularis was the first tick genome to be sequenced and annotated. However, the first assembly, I. scapularis Wikel (IscaW), was highly fragmented because of the technical challenges posed by the long, repetitive genome sequences characteristic of arthropod genomes and the lack of long-read sequencing techniques. Although I. scapularis has emerged as a model for tick research because of the availability of new tools such as embryo injection and CRISPR-Cas9-mediated gene editing yet the lack of chromosome-scale scaffolds has slowed progress in tick biology and the development of tools for their control. Here we combine diverse technologies to produce the I. scapularis Gulia-Nuss (IscGN) genome assembly and gene set. We used DNA from eggs and male and female adult ticks and took advantage of Hi-C, PacBio HiFi sequencing, and Illumina short-read sequencing technologies to produce a chromosome-level assembly. In this work, we present the predicted pseudochromosomes consisting of 13 autosomes and the sex pseudochromosomes: X and Y, and a markedly improved genome annotation compared with the existing assemblies and annotations.

High levels of intra-strain structural variation in Drosophila simulans X pericentric heterochromatin.

Large genome structural variations can impact genome regulation and integrity. Repeat-rich regions like pericentric heterochromatin are vulnerable to structural rearrangements although we know little about how often these rearrangements occur over evolutionary time. Repetitive genome regions are particularly difficult to study with genomic approaches, as they are missing from most genome assemblies. However, cytogenetic approaches offer a direct way to detect large rearrangements involving pericentric heterochromatin. Here, we use a cytogenetic approach to reveal large structural rearrangements associated with the X pericentromeric region of Drosophila simulans. These rearrangements involve large blocks of satellite DNA-the 500-bp and Rsp-like satellites-which colocalize in the X pericentromeric heterochromatin. We find that this region is polymorphic not only among different strains, but between isolates of the same strain from different labs, and even within individual isolates. On the one hand, our observations raise questions regarding the potential impact of such variation at the phenotypic level and our ability to control for such genetic variability. On the other hand, this highlights the very rapid turnover of the pericentric heterochromatin most likely associated with genomic instability of the X pericentromere. It represents a unique opportunity to study the dynamics of pericentric heterochromatin, the evolution of associated satellites on a very short time scale, and to better understand how structural variation arises.

An inducible CRISPR-Kill system for temporally controlled cell type-specific cell ablation in Arabidopsis thaliana.

The application of the CRISPR/Cas system as a biotechnological tool for genome editing has revolutionized plant biology. Recently, the repertoire was expanded by CRISPR-Kill, enabling CRISPR/Cas-mediated tissue engineering through genome elimination by tissue-specific expression. Using the Cas9 nuclease from Staphylococcus aureus (SaCas9), CRISPR-Kill relies on the induction of multiple double-strand breaks (DSBs) in conserved repetitive genome regions, such as the rDNA, causing cell death of the targeted cells. Here, we show that in addition to spatial control by tissue-specific expression, temporal control of CRISPR-mediated cell death is feasible in Arabidopsis thaliana. We established a chemically inducible tissue-specific CRISPR-Kill system that allows the simultaneous detection of targeted cells by fluorescence markers. As proof of concept, we were able to eliminate lateral roots and ablate root stem cells. Moreover, using a multi-tissue promoter, we induced targeted cell death at defined time points in different organs at select developmental stages. Thus, using this system makes it possible to gain new insights into the developmental plasticity of certain cell types. In addition to enabling tissue engineering in plants, our system provides an invaluable tool to study the response of developing plant tissue to cell elimination through positional signaling and cell-to-cell communication.

Diversity analyses in two ornamental and large-genome Ranunculaceae species based on a low-cost Klenow NGS-based protocol

Persian buttercup ( Ranunculus asiaticus L.) and poppy anemone ( Anemone coronaria L.) are ornamental, outcrossing, perennial species belonging to the Ranunculaceae family, characterized by large and highly repetitive genomes. We applied K-seq protocol in both species to generate high-throughput sequencing data and produce a large number of genetic polymorphisms. The technique entails the application of Klenow polymerase-based PCR using short primers designed by analyzing k-mer sets in the genome sequence. To date the genome sequence of both species has not been released, thus we designed primer sets based on the reference the genome sequence of the related species Aquilegia oxysepala var. kansuensis (Brühl). A whole of 11,542 SNPs were selected for assessing genetic diversity of eighteen commercial varieties of R. asiaticus , while 1,752 SNPs for assessing genetic diversity in six cultivars of A. coronaria . UPGMA dendrograms were constructed and in R. asiaticus integrated in with PCA analysis. This study reports the first molecular fingerprinting within Persian buttercup, while the results obtained in poppy anemone were compared with a previously published SSR-based fingerprinting, proving K-seq to be an efficient protocol for the genotyping of complex genetic backgrounds.

Low coverage sequencing provides insights into the key features of the nuclear and mitochondrial genomes of the Alligator Snapping Turtle Macrochelys temminckii

The alligator snapping turtle Macrochelys temminckii is a culturally, ecologically, and evolutionary relevant species of conservation concern. In this study, we conducted a genome survey of M. temminckii. Using a low-coverage short read sequencing strategy, this study estimated the genome size, repetitive genome content, annotated and quantified repetitive elements, assembled the 45S rRNA DNA operon, and characterized in detail the mitochondrial genome of M. temminckii. Using a k-mer strategy, the estimated haploid genome size varied between 3.77 and 3.19 Gbp, which is within the range previously reported for other representatives of the family Chelydridae. Repetitive genome content estimates using different k-mers (21 to 51) indicated that more than 75 % of the genome of M. temminckii comprised repetitive elements. Taking into account only annotated repetitive elements, the most common repetitive elements were classified as Class I - Long Interspersed Nuclear Element (LINE) which were more abundant than Class I – Penelope and Class I - Long Terminal Repeat (LTR) Ty3 mobile elements. Less abundant repeat element families in the nuclear genome of M. temminckii included Class I - DIRS mobile elements and Satellite DNA. The nuclear ribosomal operon was partially assembled into three contigs, one encoding the complete ssrRNA gene, a second encoding the complete 5.8S rRNA gene, and a third comprising the full lsrRNA gene. The AT-rich complete mitochondrial genome was 16,570 bp long. These new genomic resources are of utmost importance to aid in the development of conservation plans for this iconic freshwater turtle.

CRISPR-Cas9-based repeat depletion for high-throughput genotyping of complex plant genomes.

High-throughput genotyping enables the large-scale analysis of genetic diversity in population genomics and genome-wide association studies that combine the genotypic and phenotypic characterization of large collections of accessions. Sequencing-based approaches for genotyping are progressively replacing traditional genotyping methods because of the lower ascertainment bias. However, genome-wide genotyping based on sequencing becomes expensive in species with large genomes and a high proportion of repetitive DNA. Here we describe the use of CRISPR-Cas9 technology to deplete repetitive elements in the 3.76-Gb genome of lentil (Lens culinaris), 84% consisting of repeats, thus concentrating the sequencing data on coding and regulatory regions (single-copy regions). We designed a custom set of 566,766 gRNAs targeting 2.9 Gbp of repeats and excluding repetitive regions overlapping annotated genes and putative regulatory elements based on ATAC-seq data. The novel depletion method removed ∼40% of reads mapping to repeats, increasing those mapping to single-copy regions by ∼2.6-fold. When analyzing 25 million fragments, this repeat-to-single-copy shift in the sequencing data increased the number of genotyped bases of ∼10-fold compared to nondepleted libraries. In the same condition, we were also able to identify ∼12-fold more genetic variants in the single-copy regions and increased the genotyping accuracy by rescuing thousands of heterozygous variants that otherwise would be missed because of low coverage. The method performed similarly regardless of the multiplexing level, type of library or genotypes, including different cultivars and a closely related species (L. orientalis). Our results showed that CRISPR-Cas9-driven repeat depletion focuses sequencing data on single-copy regions, thus improving high-density and genome-wide genotyping in large and repetitive genomes.

The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

A practical comparison of the next-generation sequencing platform and assemblers using yeast genome.

Assembling fragmented whole-genomic information from the sequencing data is an inevitable process for further genome-wide research. However, it is intricate to select the appropriate assembly pipeline for unknown species because of the species-specific genomic properties. Therefore, our study focused on relatively more static proclivities of sequencing platforms and assembly algorithms than the fickle genome sequences. A total of 212 draft and polished de novo assemblies were constructed under the different sequencing platforms and assembly algorithms with the repetitive yeast genome. Our comprehensive data indicated that sequencing reads from Oxford Nanopore with R7.3 flow cells generated more continuous assemblies than those derived from the PacBio Sequel, although the homopolymer-based assembly errors and chimeric contigs exist. In addition, the comparison between two second-generation sequencing platforms showed that Illumina NovaSeq 6000 provides more accurate and continuous assembly in the second-generation-sequencing-first pipeline, but MGI DNBSEQ-T7 provides a cheap and accurate read in the polishing process. Furthermore, our insight into the relationship among the computational time, read length, and coverage depth provided clues to the optimal pipelines of yeast assembly.