Large-scale analysis reveals that the genome features of simple sequence repeats are generally conserved at the family level in insects

Evolutionary and phylogenetic insights from the mitochondrial genomic analysis of Diceraeus melacanthus and D. furcatus (Hemiptera: Pentatomidae)

The mitochondrial genomes of D. melacanthus and D. furcatus were sequenced and used to investigate the phylogenetic relationships with 54 species of Pentatomidae. Their mitogenomes are 17,197 and 15,444 bp-long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22/21 transfer RNA genes, with conserved gene arrangement. Leu, Lys, and Ser were the most common amino acids in their PCGs. PCGs evolutionary analysis indicated their mitogenomes are under purifying selection, and the most conserved genes are from the cytochrome complex, reinforcing their suitability as markers for molecular taxonomy. We identified 490 mtSSRs in 56 Pentatomidae species, with large variation and a positive correlation between mtSSR number and genome size. Three mtSSRs were identified in each Diceraeus species. Only the mtSSR in the nad6 (D. melacanthus) and nad4 (D. furcatus) appear to have application as molecular markers for species characterization. Phylogenetic analysis confirmed the monophyly of Pentatomidae. However, our analysis challenged the monophyly of Pentatominae and Podopinae. We also detected unexpected relationships among some tribes and genera, highlighting the complexity of the internal taxonomic structure of Pentatomidae. Both Diceraeus species were grouped in the same clade with the remaining Carpocorini analyzed.

Development and validation of genome-wide SSR molecular markers of Tapes dorsatus

BACKGROUND

Tapes dorsatus is an economically important benthic animal in the Beibu Gulf of China. However, the deficiency of microsatellite markers has hindered the study of its genetics. The development of microsatellite markers will provide useful tools for genetic improvement, variety identification, phylogenetic analysis and resource conservation.

METHODS AND RESULTS

Within the genome sequence, 145,008 simple sequence repeats (SSRs) were identified, and 29,691 primer pairs were designed successfully. A total of 100 primer pairs were randomly synthesized for testing, and 93 primers yielded products. Sixty highly polymorphic primers were used to reveal the genetic diversity of 50 T. dorsatus individuals. The average number of alleles (Na) of the population was 10.40; the average number of effective alleles was 6.16, the average expected heterozygosity (He) was 0.82, and the average polymorphic information content was 0.80. The genetic structure of the population was detected, by which the population could be divided into three subpopulations.

CONCLUSION

We identified 145,008 SSRs in the genome of T. dorsatus and designed 29,691 primer pairs in this study. Of 100 synthesized primers, 60 were highly polymorphic and used to reveal the genetic diversity and structure of the population. The SSR markers identified here will provide useful tools and a foundation for genetic diversity, linkage mapping and molecular marker-aided breeding in T. dorsatus.

Differential Spreading of Microsatellites in Holocentric Chromosomes of Chagas Disease Vectors: Genomic and Evolutionary Implications

This study focused on analyzing the distribution of microsatellites in holocentric chromosomes of the Triatominae subfamily, insect vectors of Chagas disease. We employed a non-denaturing FISH technique to determine the chromosomal distribution of sixteen microsatellites across twenty-five triatomine species, involving five genera from the two principal tribes: Triatomini and Rhodniini. Three main hybridization patterns were identified: strong signals in specific chromosomal regions, dispersed signals dependent on microsatellite abundance and the absence of signals in certain chromosomal regions or entire chromosomes. Significant variations in hybridization patterns were observed between Rhodniini and Triatomini species. Rhodniini species displayed weak and scattered hybridization signals, indicating a low abundance of microsatellites in their genomes. In contrast, Triatomini species exhibited diverse and abundant hybridization patterns, suggesting that microsatellites are a significant repetitive component in their genomes. One particularly interesting finding was the high abundance of GATA repeats, and to a lesser extent AG repeats, in the Y chromosome of all analyzed Triatomini species. In contrast, the Y chromosome of Rhodniini species did not show enrichment in GATA and AG repeats. This suggests that the richness of GATA repeats on the Y chromosome likely represents an ancestral trait specific to the Triatomini tribe. Furthermore, this information can be used to elucidate the evolutionary relationships between Triatomini and other groups of reduviids, contributing to the understanding of the subfamily's origin. Overall, this study provides a comprehensive understanding of the composition and distribution of microsatellites within Triatominae genomes, shedding light on their significance in the evolutionary processes of these species.

The Chromosome-Level Genome Assembly of Bean Blossom Thrips (Megalurothrips usitatus) Reveals an Expansion of Protein Digestion-Related Genes in Adaption to High-Protein Host Plants

Megalurothrips usitatus (Bagnall) is a destructive pest of legumes, such as cowpea. The biology, population dynamics and control strategies of this pest have been well studied. However, the lack of a high-quality reference genome for M. usitatus has hindered the understanding of key biological questions, such as the mechanism of adaptation to feed preferentially on high-protein host plants and the resistance to proteinase inhibitors (PIs). In this study, we generated a high-resolution chromosome-level reference genome assembly (247.82 Mb, 16 chromosomes) of M. usitatus by combining Oxford Nanopore Technologies (ONT) and Hi-C sequencing. The genome assembly showed higher proportions of GC and repeat content compared to other Thripinae species. Genome annotation revealed 18,624 protein-coding genes, including 4613 paralogs that were preferentially located in TE-rich regions. GO and KEGG enrichment analyses of the paralogs revealed significant enrichment in digestion-related genes. Genome-wide identification uncovered 506 putative digestion-related enzymes; of those, proteases, especially their subgroup serine proteases (SPs), are significantly enriched in paralogs. We hypothesized that the diversity and expansion of the digestion-related genes, especially SPs, could be driven by mobile elements (TEs), which promote the adaptive evolution of M. usitatus to high-protein host plants with high serine protease inhibitors (SPIs). The current study provides a valuable genomic resource for understanding the genetic variation among different pest species adapting to different plant hosts.

Comparative genomics reveals insights into anuran genome size evolution

BACKGROUND

Amphibians, particularly anurans, display an enormous variation in genome size. Due to the unavailability of whole genome datasets in the past, the genomic elements and evolutionary causes of anuran genome size variation are poorly understood. To address this, we analyzed whole-genome sequences of 14 anuran species ranging in size from 1.1 to 6.8 Gb. By annotating multiple genomic elements, we investigated the genomic correlates of anuran genome size variation and further examined whether the genome size relates to habitat types.

RESULTS

Our results showed that intron expansions or contraction and Transposable Elements (TEs) diversity do not contribute significantly to genome size variation. However, the recent accumulation of transposable elements (TEs) and the lack of deletion of ancient TEs primarily accounted for the evolution of anuran genome sizes. Our study showed that the abundance and density of simple repeat sequences positively correlate with genome size. Ancestral state reconstruction revealed that genome size exhibits a taxon-specific pattern of evolution, with families Bufonidae and Pipidae experiencing extreme genome expansion and contraction events, respectively. Our result showed no relationship between genome size and habitat types, although large genome-sized species are predominantly found in humid habitats.

CONCLUSIONS

Overall, our study identified the genomic element and their evolutionary dynamics accounting for anuran genome size variation, thus paving a path to a greater understanding of the size evolution of the genome in amphibians.

Omics in the Red Palm Weevil Rhynchophorus ferrugineus (Olivier) (Coleoptera: Curculionidae): A Bridge to the Pest

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is the most devastating pest of palm trees worldwide. Mitigation of the economic and biodiversity impact it causes is an international priority that could be greatly aided by a better understanding of its biology and genetics. Despite its relevance, the biology of the RPW remains poorly understood, and research on management strategies often focuses on outdated empirical methods that produce sub-optimal results. With the development of omics approaches in genetic research, new avenues for pest control are becoming increasingly feasible. For example, genetic engineering approaches become available once a species's target genes are well characterized in terms of their sequence, but also population variability, epistatic interactions, and more. In the last few years alone, there have been major advances in omics studies of the RPW. Multiple draft genomes are currently available, along with short and long-read transcriptomes, and metagenomes, which have facilitated the identification of genes of interest to the RPW scientific community. This review describes omics approaches previously applied to RPW research, highlights findings that could be impactful for pest management, and emphasizes future opportunities and challenges in this area of research.

Comparison of microsatellite distribution in the genomes of Pteropus vampyrus and Miniopterus natalensis (Chiroptera)

BACKGROUND

Microsatellites are a ubiquitous occurrence in prokaryotic and eukaryotic genomes. Microsatellites have become one of the most popular classes of genetic markers due to their high reproducibility, multi-allelic nature, co-dominant mode of inheritance, abundance and wide genome coverage. We characterised microsatellites in the genomes and genes of two bat species, Pteropus vampyrus and Miniopterus natalensis. This characterisation was used for gene ontology analysis and the Kyoto Encyclopedia of Genes and Genomes pathway enrichment of coding sequences (CDS).

RESULTS

Compared to M. natalensis, the genome size of P. vampyrus is larger and contains more microsatellites, but the total diversity of both species is similar. Mononucleotide and dinucleotide repeats were the most diverse in the genome of the two species. In each bat species, the microsatellite bias was obvious. The microsatellites with the largest number of repeat motifs in P. vampyrus from mononucleotide to hexanucleotide were (A)_n, (AC)_n, (CAA)_n, (AAAC)_n, (AACAA)_n and (AAACAA)_n, with frequencies of 97.94%, 58.75%, 30.53%, 22.82%, 54.68% and 22.87%, respectively, while in M. natalensis were (A)_n, (AC)_n, (TAT)_n, (TTTA)_n, (AACAA)_n and (GAGAGG)_n, with of 92.00%, 34.08%, 40.36%, 21.83%, 25.42% and 12.79%, respectively. In both species, the diversity of microsatellites was highest in intergenic regions, followed by intronic, untranslated and exonic regions and lowest in coding regions. Location analysis indicated that microsatellites were mainly concentrated at both ends of the genes. Microsatellites in the CDS are thus subject to higher selective pressure. In the GO analysis, two unique GO terms were found only in P. vampyrus and M. natalensis, respectively. In KEGG enriched pathway, the biosynthesis of other secondary metabolites and metabolism of other amino acids in metabolism pathways were present only in M. natalensis. The combined biological process, cellular components and molecular function ontology are reflected in the GO analysis and six functional enrichments in KEGG annotation, suggesting advantageous mutations during species evolution.

CONCLUSIONS

Our study gives a comparative characterization of the genomes of microsatellites composition in the two bat species. And also allow further study on the effect of microsatellites on gene function as well as provide an insight into the molecular basis for species adaptation to new and changing environments.

Development of Nuclear DNA Markers for Applications in Genetic Diversity Study of Oil Palm-Pollinating Weevil Populations

The oil palm-pollinating weevil (Elaeidobius kamerunicus Faust) was introduced from Cameroon, West Africa, to Malaysia in 1981, and subsequently, to other oil palm-growing countries as well. This study aims to develop a set of robust E. kamerunicus-specific nuclear DNA markers to directly assess the genetic diversity of the weevil populations. A total of 19,148 SNP and 223,200 SSR were discovered from 48 weevils representing three origins (Peninsular Malaysia, Sabah, and Riau) using RAD tag sequencing. Subsequent filtering steps further reduced these to 1000 SNP and 120 SSR. The selected 220 SNP exhibited a polymorphism information content (PIC) of 0.2387 (±0.1280), and 8 SSR had the PIC of 0.5084 (±0.1928). These markers were found to show sufficient polymorphism, making it possible to assign 180 weevils into three major clusters from Ghana, Cameroon, and Southeast Asia (mainly in Malaysia and Indonesia). These DNA markers successfully confirmed the Cameroon origin of the Southeast Asian cluster. However, the presence of null alleles in the SSR markers, due to limited flexibility of the probe design on the short RAD tags, led to an underestimation of heterozygosity within the populations. Hence, the developed SNP markers turned out to be more efficient than the SSR markers in the genetic diversity assessment of the E. kamerunicus populations. The genetic information provides useful insight into developing guidelines for the genetic monitoring and conservation planning of E. kamerunicus.

Mutation and selection processes regulating short tandem repeats give rise to genetic and phenotypic diversity across species

Short tandem repeats (STRs) are units of 1-6 bp that repeat in a tandem fashion in DNA. Along with single nucleotide polymorphisms and large structural variations, they are among the major genomic variants underlying genetic, and likely phenotypic, divergence. STRs experience mutation rates that are orders of magnitude higher than other well-studied genotypic variants. Frequent copy number changes result in a wide range of alleles, and provide unique opportunities for modulating complex phenotypes through variation in repeat length. While classical studies have identified key roles of individual STR loci, the advent of improved sequencing technology, high-quality genome assemblies for diverse species, and bioinformatics methods for genome-wide STR analysis now enable more systematic study of STR variation across wide evolutionary ranges. In this review, we explore mutation and selection processes that affect STR copy number evolution, and how these processes give rise to varying STR patterns both within and across species. Finally, we review recent examples of functional and adaptive changes linked to STRs.

Genetic structure of Spirometra mansoni (Cestoda: Diphyllobothriidae) populations in China revealed by a Target SSR-seq method

BACKGROUND

In China, the plerocercoid of the cestode Spirometra mansoni is the main causative agent of human and animal sparganosis. However, the population genetic structure of this parasite remains unclear. In this study, we genotyped S. mansoni isolates with the aim to improve current knowledge on the evolution and population diversity of this cestode.

METHODS

We first screened 34 perfect simple sequence repeats (SSRs) using all available omic data and then constructed target sequencing technology (Target SSR-seq) based on the Illumina NovaSeq platform. Next, a series of STRUCTURE. clustering, principal component, analysis of molecular variance and TreeMix analyses were performed on 362 worm samples isolated from 12 different hosts in 16 geographical populations of China to identify the genetic structure.

RESULTS

A total of 170 alleles were detected. The whole population could be organized and was found to be derived from the admixture of two ancestral clusters. TreeMix analysis hinted that possible gene flow occurred from Guizhou (GZ) to Sichuan (SC), SC to Jaingxi (JX), SC to Hubei (HB), GZ to Yunnan (YN) and GZ to Jiangsu (JS). Both neighbor-joining clustering and principal coordinate analysis showed that isolates from intermediate hosts tend to cluster together, while parasites from definitive hosts revealed greater genetic differences. Generally, a S. mansoni population was observed to harbor high genetic diversity, moderate genetic differentiation and a little genetic exchange among geographical populations.

CONCLUSIONS

A Target SSR-seq genotyping method was successfully developed, and an in-depth view of genetic diversity and genetic relationship will have important implications for the prevention and control of sparganosis.

Microsatellite Variation in the Most Devastating Beetle Pests (Coleoptera: Curculionidae) of Agricultural and Forest Crops

Weevils, classified in the family Curculionidae (true weevils), constitute a group of phytophagous insects of which many species are considered significant pests of crops. Within this family, the red palm weevil (RPW), Rhynchophorus ferrugineus, has an integral role in destroying crops and has invaded all countries of the Middle East and many in North Africa, Southern Europe, Southeast Asia, Oceania, and the Caribbean Islands. Simple sequence repeats (SSRs), also termed microsatellites, have become the DNA marker technology most applied to study population structure, evolution, and genetic diversity. Although these markers have been widely examined in many mammalian and plant species, and draft genome assemblies are available for many species of true weevils, very little is yet known about SSRs in weevil genomes. Here we carried out a comparative analysis examining and comparing the relative abundance, relative density, and GC content of SSRs in previously sequenced draft genomes of nine true weevils, with an emphasis on R. ferrugineus. We also used Illumina paired-end sequencing to generate draft sequence for adult female RPW and characterized it in terms of perfect SSRs with 1–6 bp nucleotide motifs. Among weevil genomes, mono- to trinucleotide SSRs were the most frequent, and mono-, di-, and hexanucleotide SSRs exhibited the highest GC content. In these draft genomes, SSR number and genome size were significantly correlated. This work will aid our understanding of the genome architecture and evolution of Curculionidae weevils and facilitate exploring SSR molecular marker development in these species.

Development of SSR markers and genetic diversity analysis based on RAD-seq technology among Chinese populations of Daphnia magna

BACKGROUND

Daphnia magna belongs to the Cladocera order and plays an important role in the aquatic ecosystem. With the intensification of water pollution, the wild population of D. magna has declined rapidly in recent years, and insufficient molecular markers have limited effective research and conservation of this species.

METHODS AND RESULTS

26 novel microsatellite (SSR) markers were developed in an artificially domesticated D. magna and 12 wild D. magna populations using restriction site-associated DNA sequencing (RAD-seq). The results showed that the observed heterozygosity (Ho) and expected heterozygosity (He) ranged from 0.083 to 0.999 and 0.085 to 0.862, respectively. The PIC ranged from 0.368 to 0.805. These results indicate that the developed SSR marker is highly polymorphic. Nei's genetic identity (H) ranged from 0.0926 to 0.3462. Shannon's Information index (I) ranged from 0.1333 to 0.4799. Genetic distance and Nei's genetic identity analysis, NJ tree diagram analysis, and PCoA analysis were conducted on populations of D. magna from different regions. The results show that the D. magna genetic relationship between Liaoning and Shanxi, Hunan and Anhui, and Beijing and Hainan are relatively close, while the genetic structure of D. magna in Guangdong, Jiangsu, and Sichuan is quite different from other sampling sites. An analysis of population genetic structure divided the D. magna samples into two major groups.

CONCLUSIONS

These results indicate that the genetic structure of D. magna differs considerably in different regions. Our research results and the newly developed polymorphic SSR markers for D. magna are of great significance in terms of the genetic breeding of D. magna, identification of wild and artificially domesticated populations and conservation genetics research.

Genome-Wide Survey and Development of the First Microsatellite Markers Database (AnCorDB) in Anemone coronaria L

Anemone coronaria L. (2n = 2x = 16) is a perennial, allogamous, highly heterozygous plant marketed as a cut flower or in gardens. Due to its large genome size, limited efforts have been made in order to develop species-specific molecular markers. We obtained the first draft genome of the species by Illumina sequencing an androgenetic haploid plant of the commercial line “MISTRAL^® Magenta”. The genome assembly was obtained by applying the MEGAHIT pipeline and consisted of 2 × 10⁶ scaffolds. The SciRoKo SSR (Simple Sequence Repeats)-search module identified 401.822 perfect and 188.987 imperfect microsatellites motifs. Following, we developed a user-friendly “Anemone coronaria Microsatellite DataBase” (AnCorDB), which incorporates the Primer3 script, making it possible to design couples of primers for downstream application of the identified SSR markers. Eight genotypes belonging to eight cultivars were used to validate 62 SSRs and a subset of markers was applied for fingerprinting each cultivar, as well as to assess their intra-cultivar variability. The newly developed microsatellite markers will find application in Breeding Rights disputes, developing genetic maps, marker assisted breeding (MAS) strategies, as well as phylogenetic studies.

An empirical analysis of mtSSRs: could microsatellite distribution patterns explain the evolution of mitogenomes in plants?

Microsatellites (SSRs) are tandem repeat sequences in eukaryote genomes, including plant cytoplasmic genomes. The mitochondrial genome (mtDNA) has been shown to vary in size, number, and distribution of SSRs among different plant groups. Thus, SSRs contribute with genomic diversity in mtDNAs. However, the abundance, distribution, and evolutionary significance of SSRs in mtDNA from a wide range of algae and plants have not been explored. In this study, the mtDNAs of 204 plant and algal species were investigated related to the presence of SSRs. The number of SSRs was positively correlated with genome size. Its distribution is dependent on plant and algal groups analyzed, although the cluster analysis indicates the conservation of some common motifs in algal and terrestrial plants that reflect common ancestry of groups. Many SSRs in coding and non-coding regions can be useful for molecular markers. Moreover, mitochondrial SSRs are highly abundant, representing an important source for natural or induced genetic variation, i.e., for biotechnological approaches that can modulate mtDNA gene regulation. Thus, this comparative study increases the understanding of the plant and algal SSR evolution and brings perspectives for further studies.

Insights into the origin of the invasive populations of Trioza erytreae in Europe using microsatellite markers and mtDNA barcoding approaches

The African citrus psyllid Trioza erytreae is one of the major threats to citrus industry as the vector of the incurable disease known as huanglongbing (HLB) or citrus greening. The psyllid invaded the northwest of the Iberian Peninsula 6 years ago. The invasion alarmed citrus growers in the Mediterranean basin, the largest citrus producing area in Europe, which is still free of HLB. Before our study, no research had been carried out on the genetic diversity of T. erytreae populations that have invaded the Iberian Peninsula and the archipelagos of the Macaronesia (Madeira and the Canary Islands). In this study, combining microsatellites markers and mtDNA barcoding analysis, we characterize the genetic diversity, structure and maternal relationship of these new invasive populations of T. erytreae and those from Africa. Our results suggest that the outbreaks of T. erytreae in the Iberian Peninsula may have derived from the Canary Islands. The populations of T. erytreae that invaded Macaronesia and the Iberian Peninsula are likely to have originated from southern Africa. We anticipate our results to be a starting point for tracking the spread of this invasive pest outside of Africa and to be important for optimizing contingency and eradication plans in newly invaded and free areas.

The Draft Genome of Yellow Stem Borer, an Agriculturally Important Pest, Provides Molecular Insights into Its Biology, Development and Specificity Towards Rice for Infestation

Simple Summary

Yellow stem borer (YSB), is the most destructive and widely occurring pest that attacks rice throughout the growing season. Rice (Oryza sativa L.) is a major staple cereal worldwide, providing essential caloric requirements for more than half of the world’s population. Annual losses to rice borers are approximately 5–10%, but losses in individual fields may reach up to 50–60%. The use of traditional pest management strategies in controlling YSB is somewhat challenging due to its unique internal feeding habit. Genome sequence information of economically important crop pests is important for designing or developing pest-resistant rice varieties. In an approach to achieve this, we present our first-ever study on the draft genome sequence of YSB. The information provided from our current study might be useful in developing genome-based approaches for the management of pest species.

Abstract

Yellow stem borer (YSB), Scirpophaga incertulas (Walker) (Lepidoptera: Crambidae), a major monophagous insect pest of rice, causes significant yield losses. The rice–YSB interaction is very dynamic, making it difficult for management. The development of resistant lines has been unsuccessful as there are no effective resistant sources in the germplasm. Genome information is necessary for a better understanding of interaction with rice in terms of its recognition, response, and infestation mechanism. The draft genome of YSB is predicted to have 46,057 genes with an estimated size of 308 Mb, being correlated with the flow cytometry analysis. The existence of complex metabolic mechanisms and genes related to specific behavior was identified, being conditioned by a higher level of regulation. We deciphered the possible visual, olfactory, and gustatory mechanisms responsible for its evolution as a monophagous pest. Comparative genomic analysis revealed that YSB is unique in the way it has evolved. The obvious presence of high-immunity-related genes, well-developed RNAi machinery, and diverse effectors provides a means for developing genomic tools for its management. The identified 21,696 SSR markers can be utilized for diversity analysis of populations across the rice-growing regions. We present the first draft genome of YSB. The information emanated paves a way for biologists to design novel pest management strategies as well as for the industry to design new classes of safer and specific insecticide molecules.

Comparison of the Microsatellite Distribution Patterns in the Genomes of Euarchontoglires at the Taxonomic Level

Microsatellite or simple sequence repeat (SSR) instability within genes can induce genetic variation. The SSR signatures remain largely unknown in different clades within Euarchontoglires, one of the most successful mammalian radiations. Here, we conducted a genome-wide characterization of microsatellite distribution patterns at different taxonomic levels in 153 Euarchontoglires genomes. Our results showed that the abundance and density of the SSRs were significantly positively correlated with primate genome size, but no significant relationship with the genome size of rodents was found. Furthermore, a higher level of complexity for perfect SSR (P-SSR) attributes was observed in rodents than in primates. The most frequent type of P-SSR was the mononucleotide P-SSR in the genomes of primates, tree shrews, and colugos, while mononucleotide or dinucleotide motif types were dominant in the genomes of rodents and lagomorphs. Furthermore, (A)n was the most abundant motif in primate genomes, but (A)n, (AC)n, or (AG)n was the most abundant motif in rodent genomes which even varied within the same genus. The GC content and the repeat copy numbers of P-SSRs varied in different species when compared at different taxonomic levels, reflecting underlying differences in SSR mutation processes. Notably, the CDSs containing P-SSRs were categorized by functions and pathways using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes annotations, highlighting their roles in transcription regulation. Generally, this work will aid future studies of the functional roles of the taxonomic features of microsatellites during the evolution of mammals in Euarchontoglires.

The Developmental Transcriptome of Bagworm, Metisa plana (Lepidoptera: Psychidae) and Insights into Chitin Biosynthesis Genes

Bagworm, Metisa plana (Lepidoptera: Psychidae) is a ubiquitous insect pest in the oil palm plantations. M. plana infestation could reduce the oil palm productivity by 40% if it remains untreated over two consecutive years. Despite the urgency to tackle this issue, the genome and transcriptome of M. plana have not yet been fully elucidated. Here, we report a comprehensive transcriptome dataset from four different developmental stages of M. plana, comprising of egg, third instar larva, pupa and female adult. The de novo transcriptome assembly of the raw data had produced a total of 193,686 transcripts, which were then annotated against UniProt, NCBI non-redundant (NR) database, Gene Ontology, Cluster of Orthologous Group, and Kyoto Encyclopedia of Genes and Genomes databases. From this, 46,534 transcripts were annotated and mapped to 146 known metabolic or signalling KEGG pathways. The paper further identified 41 differentially expressed transcripts encoding seven genes in the chitin biosynthesis pathways, and their expressions across each developmental stage were further analysed. The genetic diversity of M. plana was profiled whereby there were 21,516 microsatellite sequences and 379,895 SNPs loci found in the transcriptome of M. plana. These datasets add valuable transcriptomic resources for further study of developmental gene expression, transcriptional regulations and functional gene activities involved in the development of M. plana. Identification of regulatory genes in the chitin biosynthesis pathway may also help in developing an RNAi-mediated pest control management by targeting certain pathways, and functional studies of the genes in M. plana.

Genome-enabled insights into the biology of thrips as crop pests

BACKGROUND

The western flower thrips, Frankliniella occidentalis (Pergande), is a globally invasive pest and plant virus vector on a wide array of food, fiber, and ornamental crops. The underlying genetic mechanisms of the processes governing thrips pest and vector biology, feeding behaviors, ecology, and insecticide resistance are largely unknown. To address this gap, we present the F. occidentalis draft genome assembly and official gene set.

RESULTS

We report on the first genome sequence for any member of the insect order Thysanoptera. Benchmarking Universal Single-Copy Ortholog (BUSCO) assessments of the genome assembly (size = 415.8 Mb, scaffold N50 = 948.9 kb) revealed a relatively complete and well-annotated assembly in comparison to other insect genomes. The genome is unusually GC-rich (50%) compared to other insect genomes to date. The official gene set (OGS v1.0) contains 16,859 genes, of which ~ 10% were manually verified and corrected by our consortium. We focused on manual annotation, phylogenetic, and expression evidence analyses for gene sets centered on primary themes in the life histories and activities of plant-colonizing insects. Highlights include the following: (1) divergent clades and large expansions in genes associated with environmental sensing (chemosensory receptors) and detoxification (CYP4, CYP6, and CCE enzymes) of substances encountered in agricultural environments; (2) a comprehensive set of salivary gland genes supported by enriched expression; (3) apparent absence of members of the IMD innate immune defense pathway; and (4) developmental- and sex-specific expression analyses of genes associated with progression from larvae to adulthood through neometaboly, a distinct form of maturation differing from either incomplete or complete metamorphosis in the Insecta.

CONCLUSIONS

Analysis of the F. occidentalis genome offers insights into the polyphagous behavior of this insect pest that finds, colonizes, and survives on a widely diverse array of plants. The genomic resources presented here enable a more complete analysis of insect evolution and biology, providing a missing taxon for contemporary insect genomics-based analyses. Our study also offers a genomic benchmark for molecular and evolutionary investigations of other Thysanoptera species.

Genetic Diversity and Population Structure of Cannabis Based on the Genome-Wide Development of Simple Sequence Repeat Markers

Cannabis has been used as a source of nutrition, medicine, and fiber. However, lack of genomic simple sequence repeat (SSR) markers had limited the genetic research on Cannabis species. In the present study, 92,409 motifs were identified, and 63,707 complementary SSR primer pairs were developed. The most abundant SSR motifs had six repeat units (36.60%). The most abundant type of motif was dinucleotides (70.90%), followed by trinucleotides, tetranucleotides, and pentanucleotides. We randomly selected 80 pairs of genomic SSR markers, of which 69 (86.25%) were amplified successfully; 59 (73.75%) of these were polymorphic. Genetic diversity and population structure were estimated using the 59 (72 loci) validated polymorphic SSRs and three phenotypic markers. Three hundred ten alleles were identified, and the major allele frequency ranged from 0.26 to 0.85 (average: 0.56), Nei’s genetic diversity ranged from 0.28 to 0.82 (average: 0.56), and the expected heterozygosity ranged from 0.28 to 0.81 (average: 0.56). The polymorphism information content ranged from 0.25 to 0.79 (average: 0.50), the observed number of alleles ranged from 2 to 8 (average: 4.13), and the effective number of alleles ranged from 0.28 to 0.81 (average: 0.5). The Cannabis population did not show mutation-drift equilibrium following analysis via the infinite allele model. A cluster analysis was performed using the unweighted pair group method using arithmetic means based on genetic distances. Population structure analysis was used to divide the germplasms into two subgroups. These results provide guidance for the molecular breeding and further investigation of Cannabis.

Comparison of microsatellite distribution patterns in twenty-nine beetle genomes

Simple sequence repeats (SSRs) represent an important source of genetic variation that provides a basis for adaptation to different environments in organisms. In this study, we examined the distribution patterns of SSRs in twenty-nine beetle genomes and carried out Gene Ontology (GO) analysis of CDSs embedded with perfect SSRs (P-SSRs). The results demonstrated that imperfect SSRs (I-SSRs) represented the most abundant SSR category in beetle genomes and in different genomic regions (CDS, exon, and intron regions). The numbers of P-SSRs, I-SSRs, compound SSRs, and variable number tandem repeats were positively correlated with beetle genome size, whereas neither the frequency nor the density of the SSRs was correlated with genome size. Moreover, our results demonstrated that common genomic features of P-SSRs within the same suborder or family of Coleoptera were rare. Mono-, di-, tri-, or tetranucleotide SSRs were the most abundant P-SSR categories in beetle genomes. The preferred predominant repeat motif among the mononucleotide P-SSRs was (A)n, but the most frequent repeat motifs for other length classes varied differentially among these genomes. Furthermore, the P-SSR type with the highest GC content differed in the beetle genomes and in different genomic regions. CV (coefficient of variability) analysis demonstrated that the repeat copy numbers of P-SSRs presented relatively higher variation in introns than in CDSs and exons. The GO terms of CDSs containing P-SSRs for molecular functions were mainly enriched in "binding" and "transcription". Our findings will be useful for studying the functional roles of microsatellite heterogeneity in beetle adaptation.

A Novel Software and Method for the Efficient Development of Polymorphic SSR Loci Based on Transcriptome Data

Traditional methods for developing polymorphic microsatellite loci without reference sequences are time-consuming and labor-intensive, and the polymorphisms of simple sequence repeat (SSR) loci developed from expressed sequence tag (EST) databases are generally poor. To address this issue, in this study, we developed a new software (PSSRdt) and established an effective method for directly obtaining polymorphism details of SSR loci by analyzing diverse transcriptome data. The new method includes three steps, raw data processing, PSSRdt application, and loci extraction and verification. To test the practicality of the method, we successfully obtained 1940 potential polymorphic SSRs from the transcript dataset combined with 44 pea aphid transcriptomes. Fifty-two SSR loci obtained by the new method were selected for validating the polymorphic characteristics by genotyping in pea aphid individuals. The results showed that over 92% of SSR loci were polymorphic and 73.1% of loci were highly polymorphic. Our new software and method provide an innovative approach to microsatellite development based on RNA-seq data, and open a new path for the rapid mining of numerous loci with polymorphism to add to the body of research on microsatellites.

Genome-wide mining of microsatellites in king cobra (Ophiophagus hannah) and cross-species development of tetranucleotide SSR markers in Chinese cobra (Naja atra)

The complete genome sequence provides the opportunity for genome-wide and coding region analysis of SSRs in the king cobra and for cross-species identification of microsatellite markers in the Chinese cobra. In the Ophiophagus hannah genome, tetranucleotide repeats (38.03%) were the most abundant category, followed by dinucleotides (23.03%), pentanucleotides (13.07%), mononucleotides (11.78%), trinucleotides (11.49%) and hexanucleotides (2.6%). Twenty predominant motifs in the O. hannah genome were (A)n (C)n, (AC)n, (AG)n, (AT)n, (AGG)n, (AAT)n, (AAG)n, (AAC)n, (ATG)n, (ATAG)n, (AAGG)n, (ATCT)n, (CCTT)n, (ATTT)n, (AAAT)n, (AATAG)n, (ATTCT)n, (ATATGT)n, (AGATAT)n. In total, 4344 SSRs were found in coding sequences (CDSs). Tetranucleotides (52.79%) were the most abundant microsatellite type in CDS, followed by trinucleotides (28.50%), dinucleotides (11.02%), pentanucleotides (4.42%), mononucleotides (1.77%), and hexanucleotides (1.50%). A total of 984 CDSs containing microsatellites were assigned 11152 Gene Ontology (GO) functional terms. Gene Ontology (GO) analysis demonstrated that cellular process, cell and binding were the most frequent GO terms in biological process, cellular component and molecular function, respectively. Thirty-two novel highly polymorphic (PIC > 0.5) SSR markers for Naja atra were developed from cross-species amplification based on the tetranucleotide microsatellite sequences in the king cobra genome. The number of alleles (N_A) per locus had between 3 and 11 alleles with an average of 6.5, the polymorphism information content (PIC) value ranged from 0.521 to 0.858 (average = 0.707), the observed heterozygosity (Ho) of 32 microsatellite loci ranged from 0.292 to 0.875 (mean = 0.678), the expected heterozygosity (H_E) ranged from 0.561 to 0.889 (average = 0.761), and 3 microsatellite loci exhibited statistically significant departure from Hardy-Weinberg equilibrium (HWE) after Bonferroni correction (p < 0.003).

Mapping and identification of potential target genes from short-RNA seq for the control of Pieris rapae larvae

Pieris rapae is a serious pest of brassicas worldwide. We performed de novo assembly of P. rapae transcriptome by next-generation sequencing and assembled approximately 65,727,422 clean paired-end reads into 32,118 unigenes, of which 13,585 were mapped to 255 pathways in the KEGG database. A total of 6173 novel transcripts were identified from reads directly mapped to P. rapae genome. Additionally, 1490 SSRs, 301,377 SNPs, and 29,284 InDels were identified as potential molecular markers to explore polymorphism within P. rapae populations. We screened and mapped 36 transcripts related to OBP, CSP, SNMP, PBAN, and OR. We analyzed the expression profiles of 7 selected genes involved in pheromone transport and degradation by quantitative real-time PCR; these genes are sex-specific and differentially expressed in the developmental stages. Overall, the comprehensive transcriptome resources described in this study could help understand and identify molecular targets particularly reproduction-related genes for developing effective P. rapae management tools.

Patterns of microsatellite distribution across eukaryotic genomes

Background

Microsatellites, or Simple Sequence Repeats (SSRs), are short tandem repeats of 1–6 nt motifs present in all genomes. Emerging evidence points to their role in cellular processes and gene regulation. Despite the huge resource of genomic information currently available, SSRs have been studied in a limited context and compared across relatively few species.

Results

We have identified ~ 685 million eukaryotic microsatellites and analyzed their genomic trends across 15 taxonomic subgroups from protists to mammals. The distribution of SSRs reveals taxon-specific variations in their exonic, intronic and intergenic densities. Our analysis reveals the differences among non-related species and novel patterns uniquely demarcating closely related species. We document several repeats common across subgroups as well as rare SSRs that are excluded almost throughout evolution. We further identify species-specific signatures in pathogens like Leishmania as well as in cereal crops, Drosophila, birds and primates. We also find that distinct SSRs preferentially exist as long repeating units in different subgroups; most unicellular organisms show no length preference for any SSR class, while many SSR motifs accumulate as long repeats in complex organisms, especially in mammals.

Conclusions

We present a comprehensive analysis of SSRs across taxa at an unprecedented scale. Our analysis indicates that the SSR composition of organisms with heterogeneous cell types is highly constrained, while simpler organisms such as protists, green algae and fungi show greater diversity in motif abundance, density and GC content. The microsatellite dataset generated in this work provides a large number of candidates for functional analysis and for studying their roles across the evolutionary landscape.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5516-5) contains supplementary material, which is available to authorized users.

Target SSR-Seq: A Novel SSR Genotyping Technology Associate With Perfect SSRs in Genetic Analysis of Cucumber Varieties

Simple sequence repeats (SSR) - also known as microsatellites - have been used extensively in genetic analysis, fine mapping, quantitative trait locus (QTL) mapping, as well as marker-assisted selection (MAS) breeding and other techniques. Despite a plethora of studies reporting that perfect SSRs with stable motifs and flanking sequences are more efficient for genetic research, the lack of a high throughput technology for SSR genotyping has limited their use as genetic targets in many crops. In this study, we developed a technology called Target SSR-seq that combined the multiplexed amplification of perfect SSRs with high throughput sequencing. This method can genotype plenty of SSR loci in hundreds of samples with highly accurate results, due to the substantial coverage afforded by high throughput sequencing. We also detected 844 perfect SSRs based on 182 resequencing datasets in cucumber, of which 91 SSRs were selected for Target SSR-seq. Finally, 122 SSRs, including 31 SSRs for varieties identification, were used to genotype 382 key cucumber varieties readily available in Chinese markets using our Target SSR-seq method. Libraries of PCR products were constructed and then sequenced on the Illumina HiSeq X Ten platform. Bioinformatics analysis revealed that 111 filtered SSRs were accurately genotyped with an average coverage of 1289× at an extremely low cost; furthermore, 398 alleles were observed in 382 cucumber cultivars. Genetic analysis identified four populations: northern China type, southern China type, European type, and Xishuangbanna type. Moreover, we acquired a set of 16 core SSRs for the identification of 382 cucumber varieties, of which 42 were isolated as backbone cucumber varieties. This study demonstrated that Target SSR-seq is a novel and efficient method for genetic research.