De Novo Genome Assembly and Annotation of an Andean Caddisfly, Atopsyche davidsoni Sykora, 1991, a Model for Genome Research of High-Elevation Adaptations

Evolution of Opsin Genes in Caddisflies (Insecta: Trichoptera)

Insects have evolved complex and diverse visual systems in which light-sensing protein molecules called "opsins" couple with a chromophore to form photopigments. Insect photopigments group into three major gene families based on wavelength sensitivity: long wavelength (LW), short wavelength (SW), and ultraviolet wavelength (UV). In this study, we identified 123 opsin sequences from whole-genome assemblies across 25 caddisfly species (Insecta: Trichoptera). We discovered the LW opsins have the most diversity across species and form two separate clades in the opsin gene tree. Conversely, we observed a loss of the SW opsin in half of the trichopteran species in this study, which might be associated with the fact that caddisflies are active during low-light conditions. Lastly, we found a single copy of the UV opsin in all the species in this study, with one exception: Athripsodes cinereus has two copies of the UV opsin and resides within a clade of caddisflies with colorful wing patterns.

Comprehensive assessment of 11 de novo HiFi assemblers on complex eukaryotic genomes and metagenomes

Pacific Biosciences (PacBio) HiFi sequencing technology generates long reads (>10 kbp) with very high accuracy (<0.01% sequencing error). Although several de novo assembly tools are available for HiFi reads, there are no comprehensive studies on the evaluation of these assemblers. We evaluated the performance of 11 de novo HiFi assemblers on (1) real data for three eukaryotic genomes; (2) 34 synthetic data sets with different ploidy, sequencing coverage levels, heterozygosity rates, and sequencing error rates; (3) one real metagenomic data set; and (4) five synthetic metagenomic data sets with different composition abundance and heterozygosity rates. The 11 assemblers were evaluated using quality assessment tool (QUAST) and benchmarking universal single-copy ortholog (BUSCO). We also used several additional criteria, namely, completion rate, single-copy completion rate, duplicated completion rate, average proportion of largest category, average distance difference, quality value, run-time, and memory utilization. Results show that hifiasm and hifiasm-meta should be the first choice for assembling eukaryotic genomes and metagenomes with HiFi data. We performed a comprehensive benchmarking study of commonly used assemblers on complex eukaryotic genomes and metagenomes. Our study will help the research community to choose the most appropriate assembler for their data and identify possible improvements in assembly algorithms.

Characterization of the primary structure of the major silk gene, h-fibroin, across caddisfly (Trichoptera) suborders

Larvae of caddisflies (Trichoptera) produce silk to build various underwater structures allowing them to exploit a wide range of aquatic environments. The silk adheres to various substrates underwater and has high tensile strength, extensibility, and toughness and is of interest as a model for biomimetic adhesives. As a step toward understanding how the properties of underwater silk evolved in Trichoptera, we used genomic data to identify full-length sequences and characterize the primary structure of the major silk protein, h-fibroin, across the order. The h-fibroins have conserved termini and basic motif structure with high variation in repeating modules and variation in the percentage of amino acids, mainly proline. This finding might be linked to differences in mechanical properties related to the different silk usage and sets a starting point for future studies to screen and correlate amino acid motifs and other sequence features with quantifiable silk properties.

The impact of sequencing depth and relatedness of the reference genome in population genomic studies: A case study with two caddisfly species (Trichoptera, Rhyacophilidae, Himalopsyche)

Whole genome sequencing for generating SNP data is increasingly used in population genetic studies. However, obtaining genomes for massive numbers of samples is still not within the budgets of many researchers. It is thus imperative to select an appropriate reference genome and sequencing depth to ensure the accuracy of the results for a specific research question, while balancing cost and feasibility. To evaluate the effect of the choice of the reference genome and sequencing depth on downstream analyses, we used five confamilial reference genomes of variable relatedness and three levels of sequencing depth (3.5×, 7.5× and 12×) in a population genomic study on two caddisfly species: Himalopsyche digitata and H. tibetana. Using these 30 datasets (five reference genomes × three depths × two target species), we estimated population genetic indices (inbreeding coefficient, nucleotide diversity, pairwise F ST, and genome-wide distribution of F ST) based on variants and population structure (PCA and admixture) based on genotype likelihood estimates. The results showed that both distantly related reference genomes and lower sequencing depth lead to degradation of resolution. In addition, choosing a more closely related reference genome may significantly remedy the defects caused by low depth. Therefore, we conclude that population genetic studies would benefit from closely related reference genomes, especially as the costs of obtaining a high-quality reference genome continue to decrease. However, to determine a cost-efficient strategy for a specific population genomic study, a trade-off between reference genome relatedness and sequencing depth can be considered.

Genomic Insights into the Radiation-Resistant Capability of Sphingomonas qomolangmaensis S5-59T and Sphingomonas glaciei S8-45T, Two Novel Bacteria from the North Slope of Mount Everest

Mount Everest provides natural advantages to finding radiation-resistant extremophiles that are functionally mechanistic and possess commercial significance. (1) Background: Two bacterial strains, designated S5-59T and S8-45T, were isolated from moraine samples collected from the north slope of Mount Everest at altitudes of 5700m and 5100m above sea level. (2) Methods: The present study investigated the polyphasic features and genomic characteristics of S5-59^T and S8-45^T. (3) Results: The major fatty acids and the predominant respiratory menaquinone of S5-59^T and S8-45^T were summed as feature 3 (comprising C16:1 ω6c and/or C16:1 ω7c) and ubiquinone-10 (Q-10). Phylogenetic analyses based on 16S rRNA sequences and average nucleotide identity values among these two strains and their reference type strains were below the species demarcation thresholds of 98.65% and 95%. Strains S5-59^T and S8-45^T harbored great radiation resistance. The genomic analyses showed that DNA damage repair genes, such as mutL, mutS, radA, radC, recF, recN, etc., were present in the S5-59^T and S8-45^T strains. Additionally, strain S5-59^T possessed more genes related to DNA protection proteins. The pan-genome analysis and horizontal gene transfers revealed that strains of Sphingomonas had a consistently homologous genetic evolutionary radiation resistance. Moreover, enzymatic antioxidative proteins also served critical roles in converting ROS into harmless molecules that resulted in resistance to radiation. Further, pigments and carotenoids such as zeaxanthin and alkylresorcinols of the non-enzymatic antioxidative system were also predicted to protect them from radiation. (4) Conclusions: Type strains S5-59^T (=JCM 35564T =GDMCC 1.3193T) and S8-45^T (=JCM 34749T =GDMCC 1.2715T) represent two novel species of the genus Sphingomonas with the proposed name Sphingomonas qomolangmaensis sp. nov. and Sphingomonas glaciei sp. nov. The type strains, S5-59^T and S8-45^T, were assessed in a deeply genomic study of their radiation-resistant mechanisms and this thus resulted in a further understanding of their greater potential application for the development of anti-radiation protective drugs.

Transcriptome Analysis Provides Insights into Potentilla bifurca Adaptation to High Altitude

Potentilla bifurca is widely distributed in Eurasia, including the Tibetan Plateau. It is a valuable medicinal plant in the Tibetan traditional medicine system, especially for the treatment of diabetes. This study investigated the functional gene profile of Potentilla bifurca at different altitudes by RNA-sequencing technology, including de novo assembly of 222,619 unigenes from 405 million clean reads, 57.64% of which were annotated in Nr, GO, KEGG, Pfam, and Swiss-Prot databases. The most significantly differentially expressed top 50 genes in the high-altitude samples were derived from plants that responded to abiotic stress, such as peroxidase, superoxide dismutase protein, and the ubiquitin-conjugating enzyme. Pathway analysis revealed that a large number of DEGs encode key enzymes involved in secondary metabolites, including phenylpropane and flavonoids. In addition, a total of 298 potential genomic SSRs were identified in this study, which provides information on the development of functional molecular markers for genetic diversity assessment. In conclusion, this study provides the first comprehensive assessment of the Potentilla bifurca transcriptome. This provides new insights into coping mechanisms for non-model organisms surviving in harsh environments at high altitudes, as well as molecular evidence for the selection of superior medicinal plants.

Long-read HiFi sequencing correctly assembles repetitive heavy fibroin silk genes in new moth and caddisfly genomes

Insect silk is a versatile biomaterial. Lepidoptera and Trichoptera display some of the most diverse uses of silk, with varying strength, adhesive qualities, and elastic properties. Silk fibroin genes are long (>20 Kbp), with many repetitive motifs that make them challenging to sequence. Most research thus far has focused on conserved N- and C-terminal regions of fibroin genes because a full comparison of repetitive regions across taxa has not been possible. Using the PacBio Sequel II system and SMRT sequencing, we generated high fidelity (HiFi) long-read genomic and transcriptomic sequences for the Indianmeal moth (Plodia interpunctella) and genomic sequences for the caddisfly Eubasilissa regina. Both genomes were highly contiguous (N50  = 9.7 Mbp/32.4 Mbp, L50  = 13/11) and complete (BUSCO complete  = 99.3%/95.2%), with complete and contiguous recovery of silk heavy fibroin gene sequences. We show that HiFi long-read sequencing is helpful for understanding genes with long, repetitive regions.