Genomic resources for the Yellowfin tuna Thunnus albacares.

BACKGROUND: The Yellowfin tuna (Thunnus albacares) is a large tuna exploited by major fisheries in tropical and subtropical waters of all oceans except the Mediterranean Sea. Genomic studies of population structure, adaptive variation or of the genetic basis of phenotypic traits are needed to inform fisheries management but are currently limited by the lack of a reference genome for this species. Here we report a draft genome assembly and a linkage map for use in genomic studies of T. albacares. METHODS AND RESULTS: Illumina and PacBio SMRT sequencing were used in combination to generate a hybrid assembly that comprises 743,073,847 base pairs contained in 2,661 scaffolds. The assembly has a N50 of 351,587 and complete and partial BUSCO scores of 86.47% and 3.63%, respectively. Double-digest restriction associated DNA (ddRAD) was used to genotype the 2 parents and 164 of their F1 offspring resulting from a controlled breeding cross, retaining 19,469 biallelic single nucleotide polymorphism (SNP) loci. The SNP loci were used to construct a linkage map that features 24 linkage groups that represent the 24 chromosomes of yellowfin tuna. The male and female maps span 1,243.8 cM and 1,222.9 cM, respectively. The map was used to anchor the assembly in 24 super-scaffolds that contain 79% of the yellowfin tuna genome. Gene prediction identified 46,992 putative genes 20,203 of which could be annotated via gene ontology. CONCLUSIONS: The draft reference will be valuable to interpret studies of genome wide variation in T. albacares and other Scombroid species.

A consensus and saturated genetic map provides insight into genome anchoring, synteny of Solanaceae and leaf- and fruit-related QTLs in wolfberry (Lycium Linn.).

BACKGROUND: Lycium Linn. (Solanaceae) is a genus of economically important plants producing fruits and leaves with high nutritional value and medicinal benefits. However, genetic analysis of this plant and molecular breeding for quality improvement are limited by the lack of sufficient molecular markers. RESULTS: In this study, two parental strains, 'Ningqi No. 1' (Lycium barbarum L.) and 'Yunnan Gouqi' (Lycium yunnanense Kuang et A.M. Lu), and 200 F1 hybrid individuals were resequenced for genetic analysis. In total, 8,507 well-selected SNPs were developed, and a high-density genetic map (NY map) was constructed with a total genetic distance of 2,122.24 cM. A consensus genetic map was established by integrating the NY map and a previously published genetic map (NC map) containing 15,240 SNPs, with a total genetic distance of 3,058.19 cM and an average map distance of 0.21 cM. The 12 pseudochromosomes of the Lycium reference genome were anchored using this consensus genetic map, with an anchoring rate of 64.3%. Moreover, weak collinearities between the consensus map and the pepper, potato, and tomato genomes were observed. Twenty-five stable QTLs were identified for leaf- and fruit-related phenotypes, including fruit weight, fruit longitude, leaf length, the fruit index, and the leaf index; these stable QTLs were mapped to four different linkage groups, with LOD scores ranging from 2.51 to 19.37 and amounts of phenotypic variance explained from 6.2% to 51.9%. Finally, 82 out of 188 predicted genes underlying stable QTLs for fruit-related traits were differentially expressed according to RNA-seq analysis. CONCLUSIONS: A chromosome-level assembly can provide a foundation for further functional genomics research for wolfberry. The genomic regions of these stably expressed QTLs could be used as targets for further fine mapping and development of molecular markers for marker-assisted selection (MAS). The present study provided valuable information on saturated SNP markers and reliable QTLs for map-based cloning of functional genes related to yield and morphological traits in Lycium spp.

Inferring synteny between genome assemblies: a systematic evaluation.

BACKGROUND: Genome assemblies across all domains of life are being produced routinely. Initial analysis of a new genome usually includes annotation and comparative genomics. Synteny provides a framework in which conservation of homologous genes and gene order is identified between genomes of different species. The availability of human and mouse genomes paved the way for algorithm development in large-scale synteny mapping, which eventually became an integral part of comparative genomics. Synteny analysis is regularly performed on assembled sequences that are fragmented, neglecting the fact that most methods were developed using complete genomes. It is unknown to what extent draft assemblies lead to errors in such analysis. RESULTS: We fragmented genome assemblies of model nematodes to various extents and conducted synteny identification and downstream analysis. We first show that synteny between species can be underestimated up to 40% and find disagreements between popular tools that infer synteny blocks. This inconsistency and further demonstration of erroneous gene ontology enrichment tests raise questions about the robustness of previous synteny analysis when gold standard genome sequences remain limited. In addition, assembly scaffolding using a reference guided approach with a closely related species may result in chimeric scaffolds with inflated assembly metrics if a true evolutionary relationship was overlooked. Annotation quality, however, has minimal effect on synteny if the assembled genome is highly contiguous. CONCLUSIONS: Our results show that a minimum N50 of 1 Mb is required for robust downstream synteny analysis, which emphasizes the importance of gold standard genomes to the science community, and should be achieved given the current progress in sequencing technology.

Chromosome-Level Genome Assembly of Papilio elwesi Leech, 1889 (Lepidoptera: Papilionidae).

A rarely seen butterfly species, the large swallowtail butterfly Papilio elwesi Leech, 1889 (Lepidoptera: Papilionidae), endemic to the Chinese mainland, has been declared a state-protected animal in China since 2000, but its genome is not yet available. To obtain high-quality genome assembly and annotation, we sequenced the genome and transcriptome of P. elwesi using the PacBio and PromethION platforms, respectively. The final assembled genome was 358.51 Mb, of which 97.59% was anchored to chromosomes (30 autosomes and 1 Z sex chromosome), with a contig/scaffold N50 length of 6.79/12.32 Mb and 99.0% (n = 1367) BUSCO completeness. The genome annotation pointed to 36.82% (131.99 Mb) repetitive elements and 1296 non-coding RNAs in the genome, along with 13,681 protein-coding genes that cover 98.6% (1348) of the BUSCO genes. Among the 11,499 identified gene families, 104 underwent significantly rapid expansions or contractions, and these rapidly expanding families play roles in detoxification and metabolism. Additionally, strong synteny exists between the chromosomes of P. elwesi and P. machaon. The chromosome-level genome of P. elwesi could serve as an important genomic resource for furthering our understanding of butterfly evolution and for more in-depth genomic analyses.

The Evolution of Mitochondrial Genomes between Two Cymbidium Sister Species: Dozens of Circular Chromosomes and the Maintenance and Deterioration of Genome Synteny.

Plant mitochondrial genomes (mitogenomes) exhibit fluid genome architectures, which could lead to the rapid erosion of genome synteny over a short evolutionary time scale. Among the species-rich orchid family, the leafy Cymbidium lancifolium and leafless Cymbidium macrorhizon are sister species with remarkable differences in morphology and nutritional physiology. Although our understanding of the evolution of mitochondria is incomplete, these sister taxa are ideal for examining this subject. In this study, the complete mitogenomes of C. lancifolium and C. macrorhizon, totaling 704,244 bp and 650,751 bp, respectively, were assembled. In the 2 mitogenomes, 38 protein-coding genes, 18 cis- and 6 trans-spliced introns, and approximately 611 Kb of homologous sequences are identical; overall, they have 99.4% genome-wide similarity. Slight variations in the mitogenomes of C. lancifolium and C. macrorhizon in repeat content (21.0 Kb and 21.6 Kb, respectively) and mitochondrial DNA of plastid origin (MIPT; 38.2 Kb and 37.5 Kb, respectively) were observed. The mitogenome architectures of C. lancifolium and C. macrorhizon are complex and comprise 23 and 22 mini-circular chromosomes, respectively. Pairwise comparisons indicate that the two mitogenomes are largely syntenic, and the disparity in chromosome numbers is likely due to repeat-mediated rearrangements among different chromosomes. Notably, approximately 93.2 Kb C. lancifolium mitochondrial sequences lack any homology in the C. macrorhizon mitogenome, indicating frequent DNA gains and losses, which accounts mainly for the size variation. Our findings provide unique insights into mitogenome evolution in leafy and leafless plants of sister species and shed light on mitogenome dynamics during the transition from mixotrophy to mycoheterotrophy.

High-quality genomes reveal significant genetic divergence and cryptic speciation in the model organism Folsomia candida (collembola).

The collembolan Folsomia candida Willem, 1902, is widely distributed throughout the world and has been frequently used as a test organism in soil ecology and ecotoxicology studies. However, it is questioned as an ideal "standard" because of differences in reproductive modes and cryptic genetic diversity between strains from various geographical origins. In this study, we obtained two high-quality chromosome-level genomes of F. candida, for a parthenogenetic strain (named FCDK, 219.08 Mb, 25,139 protein-coding genes) and a sexual strain (named FCSH, 153.09 Mb, 21,609 protein-coding genes), reannotated the genome of the parthenogenetic strain reported by Faddeeva-Vakhrusheva et al. in 2017 (named FCBL, 221.7 Mb, 25,980 protein-coding genes) and conducted comparative genomic analyses of the three strains. High genome similarities between FCDK and FCBL based on synteny, genome architecture, mitochondrial and nuclear gene sequences suggest that they are conspecific. The seven chromosomes of FCDK are each 25%-54% larger than the corresponding chromosomes of FCSH, showing obvious repetitive element expansions and large-scale inversions and translocations but no whole-genome duplication. The strain-specific genes, expanded gene families and genes in nonsyntenic chromosomal regions identified in FCDK are highly related to the broader environmental adaptation of parthenogenetic strains. In addition, FCDK has fewer strain-specific microRNAs than FCSH, and their mitochondrial and nuclear genes have diverged greatly. In conclusion, FCDK/FCBL and FCSH have accumulated independent genetic changes and evolved into distinct species after 10 million years ago. Our work provides important genomic resources for studying the mechanisms of rapidly cryptic speciation and soil arthropod adaptation to soil ecosystems.

Comparative genomics of the ADA clade within the Nostocales.

The ADA clade of Nostocales cyanobacteria, a group that is prominent in current harmful algal bloom events, now includes over 40 genome sequences with the recent addition of sixteen novel sequenced genomes (Dreher et al., Harmful Algae, 2021). Fourteen genomes are complete (closed), enabling highly detailed assessments of gene content and genome architecture. ADA genomes contain 5 rRNA operons, genes expected to support a photoautotrophic and diazotrophic lifestyle, and a varied array of genes for the synthesis of bioactive secondary metabolites. Genes for the production of the taste-and-odor compound geosmin and the four major classes of cyanotoxins - anatoxin-a, cylindrospermopsin, microcystin and saxitoxin - are represented in members of the ADA clade. Notably, the gene array for the synthesis of cylindrospermopsin by Dolichospermum sp. DET69 was located on a plasmid, raising the possibility of facile horizontal transmission. However, genes supporting independent conjugative transfer of this plasmid are lacking. Further, analysis of genomic loci containing this and other cyanotoxin gene arrays shows evidence that these arrays have long-term stability and do not appear to be genomic islands easily capable of horizontal transmission to other cells. There is considerable diversity in the gene complements of individual ADA genomes, including the variable presence of physiologically important genes: genomes in three species-level subclades lack the gas vesicle genes that facilitate a planktonic lifestyle, and, surprisingly, the genome of Cuspidothrix issatschenkoi CHARLIE-1, a reported diazotroph, lacks the genes for nitrogen fixation. Notably, phylogenetically related genomes possess limited synteny, indicating a prominent role for chromosome rearrangements during ADA strain evolution. The genomes contain abundant insertion sequences and repetitive transposase genes, which could be the main drivers of genome rearrangement through active transposition and homologous recombination. No prophages were found, and no evidence of viral infection was observed in the bloom population samples from which the genomes discussed here were derived. Phages thus seem to have a limited influence on ADA evolution.

Chromosome-level genome assembly of Japanese chestnut (Castanea crenata Sieb. et Zucc.) reveals conserved chromosomal segments in woody rosids.

Japanese chestnut (Castanea crenata Sieb. et Zucc.), unlike other Castanea species, is resistant to most diseases and wasps. However, genomic data of Japanese chestnut that could be used to determine its biotic stress resistance mechanisms have not been reported to date. In this study, we employed long-read sequencing and genetic mapping to generate genome sequences of Japanese chestnut at the chromosome level. Long reads (47.7 Gb; 71.6× genome coverage) were assembled into 781 contigs, with a total length of 721.2 Mb and a contig N50 length of 1.6 Mb. Genome sequences were anchored to the chestnut genetic map, comprising 14,973 single nucleotide polymorphisms (SNPs) and covering 1,807.8 cM map distance, to establish a chromosome-level genome assembly (683.8 Mb), with 69,980 potential protein-encoding genes and 425.5 Mb repetitive sequences. Furthermore, comparative genome structure analysis revealed that Japanese chestnut shares conserved chromosomal segments with woody plants, but not with herbaceous plants, of rosids. Overall, the genome sequence data of Japanese chestnut generated in this study is expected to enhance not only its genetics and genomics but also the evolutionary genomics of woody rosids.

Chromosome-level genome reference and genome editing of the tea geometrid.

The tea geometrid is a destructive insect pest on tea plants, which seriously affects tea production in terms of both yield and quality and causes severe economic losses. The tea geometrid also provides an important study system to address the ecological adaptive mechanisms underlying its unique host plant adaptation and protective resemblance. In this study, we fully sequenced and de novo assembled the reference genome of the tea geometrid, Ectropis grisescens, using long sequencing reads. We presented a highly continuous, near-complete genome reference (787.4 Mb; scaffold N50: 26.9 Mb), along with the annotation of 18,746 protein-coding genes and 53.3% repeat contents. Importantly, we successfully placed 97.8% of the assembly in 31 chromosomes based on Hi-C interactions and characterized the sex chromosome based on sex-biased sequencing coverage. Multiple quality-control assays and chromosome-scale synteny with the model species all supported the high quality of the presented genome reference. We focused biological annotations on gene families related to the host plant adaptation and camouflage in the tea geometrid and performed comparisons with other representative lepidopteran species. Important findings include the E. grisescens-specific expansion of CYP6 P450 genes that might be involved in metabolism of tea defensive chemicals and unexpected massive expansion of gustatory receptor gene families that suggests potential polyphagy for this tea pest. Furthermore, we developed an efficient genome editing system based on CRISPR/Cas9 technology and successfully implement mutagenesis of a Hox gene in the tea geometrid. Our study provides key genomic resources both for exploring unique mechanisms underlying the ecological adaptation of tea geometrids and for developing environment-friendly strategies for tea pest management.

A chromosome-level genome assembly of rice leaffolder, Cnaphalocrocis medinalis.

The rice leaffolder, Cnaphalocrocis medinalis Guenée (Crambidae, Lepidoptera), is an important agricultural pest that causes serious losses to rice production in rice-growing regions with high humidity and temperature. However, a lack of genomic resources limits in-depth understanding of its biological characteristics and ecological adaptation. Here, we sequenced the genome of rice leaffolder using the Illumina and PacBio platforms, yielding a genome assembly of 528.3 Mb with a contig N50 of 524.6 kb. A high percentage (96.4%) of Benchmarking Universal Single-Copy Orthologs (BUSCOs) were successfully detected, suggesting high-level completeness of the genome assembly. In total, 39.5% of the genome consists of repeat sequences and 15,045 protein-coding genes were annotated. Comparative phylogenomic analysis showed that some gene families associated with hormone biosynthesis expanded in rice leaffolder. Next, we used the Hi-C technique to produce a chromosome-level genome assembly with a scaffold N50 of 16.1 Mb by anchoring 3,248 scaffolds to 31 chromosomes. The rice leaffolder genome showed high chromosomal synteny with the genome of four other lepidopteran insects. By comparing coverage ratios from the genome resequencing of male and female pupae, we identified near intact Z and W chromosomes. The W chromosome is estimated as 20.75 Mb, which is the most complete known W chromosome in Lepidoptera. The protein-coding genes on the W chromosome were significantly enriched in metabolic pathways. In all, the high-quality genome assembly and the near-intact W chromosome of rice leaffolder should be a useful resource for the fields of insect migration, chromosome evolution and pest control.

Genome Wide Identification and Comparative Analysis of the Serpin Gene Family in Brachypodium and Barley.

Serpins (serine protease inhibitors) constitute one of the largest and most widely distributed superfamilies of protease inhibitors and have been identified in nearly all organisms. To gain significant insights, a comprehensive in silico analysis of the serpin gene family was carried out in the model plant for temperate grasses Brachypodium distachyon and barley Hordeum vulgare using bioinformatic tools at the genome level for the first time. We identified a total of 27 BdSRPs and 25 HvSRP genes in Brachypodium and barley, respectively, showing an unexpectedly high gene number in these model plants. Gene structure, conserved motifs and phylogenetic comparisons of serpin genes supported the role of duplication events in the expansion and evolution of serpin gene family. Further, purifying selection pressure was found to be a main driving force in the evolution of serpin genes. Genome synteny analysis indicated that BdSRP genes were present in syntenic regions of barley, rice, sorghum and maize, suggesting that they evolved before the divergence of these species from common ancestor. The distinct expression pattern in specific tissues further suggested a specialization of functions during development and in plant defense. These results suggest that the LR serpins (serpins with Leu-Arg residues at P2-P1') identified here can be utilized as candidates for exploitation in disease resistance, pest control and preventing stress-induced cell death. Additionally, serpins were identified that could lead to further research aimed at validating and functionally characterizing the role of potential serpin genes from other plants.

Developing a High-Quality Linkage Map for the Atlantic Killifish Fundulus heteroclitus.

Killifish (Fundulus heteroclitus) are widely distributed among different aquatic environments where they demonstrate an impressive range of highly-plastic and locally adaptive phenotypes. High-throughput sequencing has begun to unravel the mechanisms and evolutionary history of these interesting features by establishing relationships in the genotype-phenotype map. However, some genotype-phenotype analyses require a higher order of contiguity than what initial scaffolded (fragmented genome assembly where contigs have been assemble into scaffolds) genome assemblies can provide. Here, we used 5,685 high-quality RAD-Seq markers from a single mapping family to order 84% of the scaffolded genome assembly to 24 chromosomes. This serves to: 1) expand the killifish genomic toolkit, 2) estimate genome-wide recombination rates, and 3) compare genome synteny to humans and other fishes. After initially building our map, we found that the selection of thresholds for sequence data filtration highly impacted scaffold placement in the map. We outline each step of the approach that dramatically improved our map to help guide others toward more effective linkage mapping for genome assembly. Our final map supports strong conservation of genomic synteny among closely related fish species and reveals previously described chromosomal rearrangements between more distantly related clades. However, we also commonly found minor scaffold misorientations in F. heteroclitus and in other assemblies, suggesting that further mapping (such as optical mapping) is necessary for finer scale resolution of genome structure. Lastly, we discuss the problems that would be expected from misoriented/unplaced scaffolds and stress the importance of a quality mapped genome as a key feature for further investigating population and comparative genomic questions with F. heteroclitus and other taxa.

Draft Sequencing of the Heterozygous Diploid Genome of Satsuma (Citrus unshiu Marc.) Using a Hybrid Assembly Approach.

Satsuma (Citrus unshiu Marc.) is one of the most abundantly produced mandarin varieties of citrus, known for its seedless fruit production and as a breeding parent of citrus. De novo assembly of the heterozygous diploid genome of Satsuma ("Miyagawa Wase") was conducted by a hybrid assembly approach using short-read sequences, three mate-pair libraries, and a long-read sequence of PacBio by the PLATANUS assembler. The assembled sequence, with a total size of 359.7 Mb at the N50 length of 386,404 bp, consisted of 20,876 scaffolds. Pseudomolecules of Satsuma constructed by aligning the scaffolds to three genetic maps showed genome-wide synteny to the genomes of Clementine, pummelo, and sweet orange. Gene prediction by modeling with MAKER-P proposed 29,024 genes and 37,970 mRNA; additionally, gene prediction analysis found candidates for novel genes in several biosynthesis pathways for gibberellin and violaxanthin catabolism. BUSCO scores for the assembled scaffold and predicted transcripts, and another analysis by BAC end sequence mapping indicated the assembled genome consistency was close to those of the haploid Clementine, pummel, and sweet orange genomes. The number of repeat elements and long terminal repeat retrotransposon were comparable to those of the seven citrus genomes; this suggested no significant failure in the assembly at the repeat region. A resequencing application using the assembled sequence confirmed that both kunenbo-A and Satsuma are offsprings of Kishu, and Satsuma is a back-crossed offspring of Kishu. These results illustrated the performance of the hybrid assembly approach and its ability to construct an accurate heterozygous diploid genome.

Genome-wide SNP identification for the construction of a high-resolution genetic map of Japanese flounder (Paralichthys olivaceus): applications to QTL mapping of Vibrio anguillarum disease resistance and comparative genomic analysis.

High-resolution genetic maps are essential for fine mapping of complex traits, genome assembly, and comparative genomic analysis. Single-nucleotide polymorphisms (SNPs) are the primary molecular markers used for genetic map construction. In this study, we identified 13,362 SNPs evenly distributed across the Japanese flounder (Paralichthys olivaceus) genome. Of these SNPs, 12,712 high-confidence SNPs were subjected to high-throughput genotyping and assigned to 24 consensus linkage groups (LGs). The total length of the genetic linkage map was 3,497.29 cM with an average distance of 0.47 cM between loci, thereby representing the densest genetic map currently reported for Japanese flounder. Nine positive quantitative trait loci (QTLs) forming two main clusters for Vibrio anguillarum disease resistance were detected. All QTLs could explain 5.1-8.38% of the total phenotypic variation. Synteny analysis of the QTL regions on the genome assembly revealed 12 immune-related genes, among them 4 genes strongly associated with V. anguillarum disease resistance. In addition, 246 genome assembly scaffolds with an average size of 21.79 Mb were anchored onto the LGs; these scaffolds, comprising 522.99 Mb, represented 95.78% of assembled genomic sequences. The mapped assembly scaffolds in Japanese flounder were used for genome synteny analyses against zebrafish (Danio rerio) and medaka (Oryzias latipes). Flounder and medaka were found to possess almost one-to-one synteny, whereas flounder and zebrafish exhibited a multi-syntenic correspondence. The newly developed high-resolution genetic map, which will facilitate QTL mapping, scaffold assembly, and genome synteny analysis of Japanese flounder, marks a milestone in the ongoing genome project for this species.

Genome variation in the hyperthermophilic archaeon Aeropyrum.

Aeropyrum spp are aerobic, heterotrophic, and hyperthermophilic marine archaea. There are two closely related Aeropyrum species, Aeropyrum camini and Aeropyrum pernix, which are isolated from geographically distinct locations. Recently, we compared their genome sequences to determine their genomic variation. They possess highly conserved small genomes, reflecting their close relationship. The entire genome similarity may result from their survival strategies in adapting to extreme environmental conditions. Meanwhile, synteny disruptions were observed in some regions including clustered regularly interspaced short palindromic repeats elements. Further, the largest portion of their non-orthologous genes were genes in the two proviral regions of A. pernix (Aeropyrum pernix spindle-shaped virus 1 and Aeropyrum pernix ovoid virus 1) or ORFans considered to be derived from viruses. Our data shows that genomic diversification of Aeropyrum spp may be substantially induced by viruses. This suggests that Aeropyrum spp may have a large pan-genome that can be extended by viruses, while each of the species shares a highly conserved small genome specializing for extreme environments.