Developing conversed microsatellite markers and their implications in evolutionary analysis of the Bemisia tabaci complex

Disparities in Genetic Diversity Drive the Population Displacement of Two Invasive Cryptic Species of the Bemisia tabaci Complex in China

Within the whitefly Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) complex, two cryptic species, namely Middle East-Asia Minor 1 (MEAM1) and Mediterranean (MED), are important invasive pests affecting global agriculture and horticulture. They were introduced into China sequentially in the mid-1990s and around 2003, respectively. Subsequently, the latter invader MED has outcompeted the earlier invader MEAM1, becoming the dominant population in the field. Although extensive studies have explored the underlying mechanisms driving this shift, the contribution of population genetics remains notably underexplored. In this study, we analyzed the genetic diversity and structure of 22 MED and 8 MEAM1 populations from various regions of China using mitochondrial DNA sequencing and microsatellite genotyping. Our results indicate low and moderate levels of genetic differentiation among geographically separate populations of MED and MEAM1, respectively. Median-joining network analysis of mtCOI gene haplotypes revealed no clear geographic structuring for either, with common haplotypes observed across provinces, although MED had more haplotypes. Comparative analyses revealed that MED presented greater genetic diversity than MEAM1 on the basis of two markers. Furthermore, analysis of molecular variance supported these findings, suggesting that while some genetic variation exists between populations, a significant amount is also present within populations. These findings reveal the population genetics of the two invasive cryptic species of the B. tabaci complex in China and suggest that the disparities in genetic diversity drive the displacement of their populations in the field. This work also provides valuable information on the genetic factors influencing the population dynamics and dominance of these invasive whitefly species.

Development and characterization of microsatellite markers for population genetics of the cocoa pod borer Conopomorpha cramerella (Snellen) (Lepidoptera: Gracillaridae)

The cocoa pod borer (CPB) Conopomorpha cramerella (Snellen) (Lepidoptera: Gracillaridae) is one of the major constraints for cocoa production in South East Asia. In addition to cultural and chemical control methods, autocidal control tactics such as the Sterile Insect Technique (SIT) could be an efficient addition to the currently control strategy, however SIT implementation will depend on the population genetics of the targeted pest. The aim of the present work was to search for suitable microsatellite loci in the genome of CPB that is partially sequenced. Twelve microsatellites were initially selected and used to analyze moths collected from Indonesia, Malaysia, and the Philippines. A quality control verification process was carried out and seven microsatellites found to be suitable and efficient to distinguish differences between CPB populations from different locations. The selected microsatellites were also tested against a closely related species, i.e. the lychee fruit borer Conopomorpha sinensis (LFB) from Vietnam and eight loci were found to be suitable. The availability of these novel microsatellite loci will provide useful tools for the analysis of the population genetics and gene flow of these pests, to select suitable CPB strains to implement the SIT.

Development and characterization of microsatellite markers for the tsetse species Glossina brevipalpis and preliminary population genetics analyses

Tsetse flies, the vectors of African trypanosomes are of key medical and economic importance and one of the constraints for the development of Africa. Tsetse fly control is one of the most effective and sustainable strategies used for controlling the disease. Knowledge about population structure and level of gene flow between neighbouring populations of the target vector is of high importance to develop appropriate strategies for implementing effective management programmes. Microsatellites are commonly used to identify population structure and assess dispersal of the target populations and have been developed for several tsetse species but were lacking for Glossina brevipalpis. In this study, we screened the genome of G. brevipalpis to search for suitable microsatellite markers and nine were found to be efficient enough to distinguish between different tsetse populations. The availability of these novel microsatellite loci will help to better understand the population biology of G. brevipalpis and to assess the level of gene flow between different populations. Such information will help with the development of appropriate strategies to implement the sterile insect technique (SIT) in the framework of an area-wide integrated pest management (AW-IPM) approach to manage tsetse populations and ultimately address the trypanosomoses problem in these targeted areas.

Development of Novel Genomic Simple Sequence Repeat (g-SSR) Markers and Their Validation for Genetic Diversity Analyses in Kalmegh [Andrographis paniculata (Burm. F.) Nees]

Kalmegh (Andrographis paniculata (Burm. F.) Nees) is one of the most important medicinal plants and has been widely explored as traditional medicine. To exploit its natural genetic diversity and initiations of molecular breeding to develop novel cultivars or varieties, developments of genomic resources are essential. Four microsatellite-enriched genomic libraries-(CT)14, (GT)12, (AG)15 and (AAC)8-were constructed using the genomic DNA of A. paniculata. Initially, 183 recombinant colonies were screened for the presence of CT, GT, AG, and AAC microsatellite repeats, out of which 47 clones found positive for the desired simple sequence repeats (SSRs). It was found that few colonies had more than one desirable SSR. Thus, a sum of 67 SSRs were designed and synthesized for their validation among 42 A. paniculata accessions. Out of the 67 SSRs used for genotyping, only 41 were found to be polymorphic. The developed set of g-SSR markers showed substantial genetic variability among the selected A. paniculata accessions, with an average polymorphic information content (PIC) value of 0.32. Neighbor-joining tree analysis, population structure analysis, analysis of molecular variance (AMOVA), and principal coordinate analysis (PCoA) illustrated the considerable genetic diversity among them. The novel g-SSR markers developed in the present study could be important genomic resources for future applications in A. paniculata.

Pattern recognition receptors in grass carp Ctenopharyngodon idella: II. Organization and expression analysis of NOD-like receptors

Nucleotide-binding domain and leucine-rich repeat containing receptors (NLRs) are a pivotal intracellular pattern recognition receptor family. However, studies on NLR genes in important economic fish grass carp (Ctenopharyngodon idella) are sporadic. The accumulations of genomic resource and transcriptomic sequences make it feasible to conduct a systematic analysis of these genes. In this study, we systematically conducted the genome-wide study of C. idella NLR genes and characterized their phylogeny, gene structure, conserved domain, evolutionary mechanism, and expression profiles post viral or bacterial challenge. A total of 65 NLR genes were identified and clustered into five subfamilies based on structural and phylogenetic features, including eight NODs (NLR-A), five NLRP-like receptors (NLR-B), forty-seven teleost-specific NLRs (NLR-C), two members with a B30.2 domain at the C-terminal (NLR-B30.2), and three additional NLRs (other NLRs). Gene structure analysis showed that NLRs were significantly different, with exon numbers from 3 to 31. Conserved domain analysis showed that most members of C. idella NLRs had additional domains besides the typical NLR domains. Gene duplication analysis indicated that the evolution of the NLR gene family was mainly related to segment duplication. mRNA expression analysis indicated that many members were differently expressed in multiple tissues post grass carp reovirus (GCRV) or Aeromonas hydrophila infection. The expression was particularly enhanced in liver post GCRV infection, and obviously lower post A. hydrophila infection than that post GCRV infection in spleen. These results provide systematic basic data for further functional studies of NLR, and insight into the immune responses of piscine fish NLRs to pathogen infections.

Insight into the microbial world of Bemisia tabaci cryptic species complex and its relationships with its host

The 37 currently recognized Bemisia tabaci cryptic species are economically important species and contain both primary and secondary endosymbionts, but their diversity has never been mapped systematically across the group. To achieve this, PacBio sequencing of full-length bacterial 16S rRNA gene amplicons was carried out on 21 globally collected species in the B. tabaci complex, and two samples from B. afer were used here as outgroups. The microbial diversity was first explored across the major lineages of the whole group and 15 new putative bacterial sequences were observed. Extensive comparison of our results with previous endosymbiont diversity surveys which used PCR or multiplex 454 pyrosequencing platforms showed that the bacterial diversity was underestimated. To validate these new putative bacteria, one of them (Halomonas) was first confirmed to be present in MED B. tabaci using Hiseq2500 and FISH technologies. These results confirmed PacBio is a reliable and informative venue to reveal the bacterial diversity of insects. In addition, many new secondary endosymbiotic strains of Rickettsia and Arsenophonus were found, increasing the known diversity in these groups. For the previously described primary endosymbionts, one Portiera Operational Taxonomic Units (OTU) was shared by all B. tabaci species. The congruence of the B. tabaci-host and Portiera phylogenetic trees provides strong support for the hypothesis that primary endosymbionts co-speciated with their hosts. Likewise, a comparison of bacterial alpha diversities, Principal Coordinate Analysis, indistinct endosymbiotic communities harbored by different species and the co-divergence analyses suggest a lack of association between overall microbial diversity with cryptic species, further indicate that the secondary endosymbiont-mediated speciation is unlikely to have occurred in the B. tabaci species group.

Microsatellite Borders and Micro-sequence Conservation in Juglans

Walnuts (Juglans spp.) are economically important nut and timber species with a worldwide distribution. Using the published Persian walnut genome as a reference for the assembly of short reads from six Juglans species and several interspecific hybrids, we identified simple sequence repeats in 12 Juglans nuclear and organellar genomes. The genome-wide distribution and polymorphisms of nuclear and organellar microsatellites (SSRs) for most Juglans genomes have not been previously studied. We compared the frequency of nuclear SSR motifs and their lengths across Juglans, and identified section-specific chloroplast SSR motifs. Primer pairs were designed for more than 60,000 SSR-containing sequences based on alignment against assembled scaffold sequences. Of the >60,000 loci, 39,000 were validated by e-PCR using unique primer pairs. We identified primers containing 100% sequence identity in multiple species. Across species, sequence identity in the SSR-flanking regions was generally low. Although SSRs are common and highly dispersed in the genome, their flanking sequences are conserved at about 90 to 95% identity within Juglans and within species. In a few rare cases, flanking sequences are identical across species of Juglans. This comprehensive report of nuclear and organellar SSRs in Juglans and the generation of validated SSR primers will be a useful resource for future genetic analyses, walnut breeding programs, high-level taxonomic evaluations, and genomic studies in Juglandaceae.

Genome-wide identification and expression analyses of WRKY transcription factor family members from chickpea (Cicer arietinum L.) reveal their role in abiotic stress-responses

BACKGROUND

WRKY proteins play a vital role in the regulation of several imperative plant metabolic processes and pathways, especially under biotic and abiotic stresses. Although WRKY genes have been characterized in various major crop plants, their identification and characterization in pulse legumes is still in its infancy. Chickpea (Cicer arietinum L.) is the most important pulse legume grown in arid and semi-arid tropics.

OBJECTIVE

In silico identification and characterization of WRKY transcription factor-encoding genes in chickpea genome.

METHODS

For this purpose, a systematic genome-wide analysis was carried out to identify the non-redundant WRKY transcription factors in the chickpea genome.

RESULTS

We have computationally identified 70 WRKY-encoding non-redundant genes which were randomly distributed on all the chickpea chromosomes except chromosome 8. The evolutionary phylogenetic analysis classified the WRKY proteins into three major groups (I, II and III) and seven sub-groups (IN, IC, IIa, IIb, IIc, IId and IIe). The gene structure analysis revealed the presence of 2-7 introns among the family members. Along with the presence of absolutely conserved signatory WRKY domain, 19 different domains were also found to be conserved in a group-specific manner. Insights of gene duplication analysis revealed the predominant role of segmental duplications for the expansion of WRKY genes in chickpea. Purifying selection seems to be operated during the evolution and expansion of paralogous WRKY genes. The transcriptome data-based in silico expression analysis revealed the differential expression of CarWRKY genes in root and shoot tissues under salt, drought, and cold stress conditions. Moreover, some of these genes showed identical expression pattern under these stresses, revealing the possibility of involvement of these genes in conserved abiotic stress-response pathways.

CONCLUSION

This genome-wide computational analysis will serve as a base to accelerate the functional characterization of WRKY TFs especially under biotic and abiotic stresses.

The Complete Plastome Sequences of Eleven Capsicum Genotypes: Insights into DNA Variation and Molecular Evolution

Members of the genus Capsicum are of great economic importance, including both wild forms and cultivars of peppers and chilies. The high number of potentially informative characteristics that can be identified through next-generation sequencing technologies gave a huge boost to evolutionary and comparative genomic research in higher plants. Here, we determined the complete nucleotide sequences of the plastomes of eight Capsicum species (eleven genotypes), representing the three main taxonomic groups in the genus and estimated molecular diversity. Comparative analyses highlighted a wide spectrum of variation, ranging from point mutations to small/medium size insertions/deletions (InDels), with accD, ndhB, rpl20, ycf1, and ycf2 being the most variable genes. The global pattern of sequence variation is consistent with the phylogenetic signal. Maximum-likelihood tree estimation revealed that Capsicum chacoense is sister to the baccatum complex. Divergence and positive selection analyses unveiled that protein-coding genes were generally well conserved, but we identified 25 positive signatures distributed in six genes involved in different essential plastid functions, suggesting positive selection during evolution of Capsicum plastomes. Finally, the identified sequence variation allowed us to develop simple PCR-based markers useful in future work to discriminate species belonging to different Capsicum complexes.

Cassava whitefly, Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) in East African farming landscapes: a review of the factors determining abundance

Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is a pest species complex that causes widespread damage to cassava, a staple food crop for millions of households in East Africa. Species in the complex cause direct feeding damage to cassava and are the vectors of multiple plant viruses. Whilst significant work has gone into developing virus-resistant cassava cultivars, there has been little research effort aimed at understanding the ecology of these insect vectors. Here we assess critically the knowledge base relating to factors that may lead to high population densities of sub-Saharan African (SSA) B. tabaci species in cassava production landscapes of East Africa. We focus first on empirical studies that have examined biotic or abiotic factors that may lead to high populations. We then identify knowledge gaps that need to be filled to deliver sustainable management solutions. We found that whilst many hypotheses have been put forward to explain the increases in abundance witnessed since the early 1990s, there are little published data and these tend to have been collected in a piecemeal manner. The most critical knowledge gaps identified were: (i) understanding how cassava cultivars and alternative host plants impact population dynamics and natural enemies; (ii) the impact of natural enemies in terms of reducing the frequency of outbreaks and (iii) the use and management of insecticides to delay the development of resistance. In addition, there are several fundamental methodologies that need to be developed and deployed in East Africa to address some of the more challenging knowledge gaps.

Tandem Action of Natural and Chemical Stressors in Stream Ecosystems: Insights from a Population Genetic Perspective

Agricultural and urban land use has dramatically increased over the last century and one consequence is the release of anthropogenic chemicals into aquatic ecosystems. One of the rarely studied consequences is the effect of land use change on internal concentrations of organic micropollutants (OMPs) in aquatic invertebrates and its effects on their genotype diversity. Here, we applied population genetic and internal concentrations of OMPs analyses to determine evolutionary implications of chemical pollution on Gammarus pulex populations from a natural and two agricultural streams. Along 14 consecutive months sampled, 26 different OMPs were quantified in G. pulex extracts with the highest number, concentration, and toxic pressure in the anthropogenically stressed stream ecosystems. Our results indicate distinct internal OMP profiles and changes in both genetic variation and genetic structure in streams affected by anthropogenic activity. Genetic variation was attributed to chemical pollution whereas changes in the genetic structure were attributed to environmental disturbances, such as changes in discharge in the impacted stream ecosystems, which worked both independently and in tandem. Finally, we conclude that human-impacted streams are subjected to severe alterations in their population genetic patterns compared to nonimpacted stream ecosystems.

Genome sequencing of the sweetpotato whitefly Bemisia tabaci MED/Q

The sweetpotato whitefly Bemisia tabaci is a highly destructive agricultural and ornamental crop pest. It damages host plants through both phloem feeding and vectoring plant pathogens. Introductions of B. tabaci are difficult to quarantine and eradicate because of its high reproductive rates, broad host plant range, and insecticide resistance. A total of 791 Gb of raw DNA sequence from whole genome shotgun sequencing, and 13 BAC pooling libraries were generated by Illumina sequencing using different combinations of mate-pair and pair-end libraries. Assembly gave a final genome with a scaffold N50 of 437 kb, and a total length of 658 Mb. Annotation of repetitive elements and coding regions resulted in 265.0 Mb TEs (40.3%) and 20 786 protein-coding genes with putative gene family expansions, respectively. Phylogenetic analysis based on orthologs across 14 arthropod taxa suggested that MED/Q is clustered into a hemipteran clade containing A. pisum and is a sister lineage to a clade containing both R. prolixus and N. lugens. Genome completeness, as estimated using the CEGMA and Benchmarking Universal Single-Copy Orthologs pipelines, reached 96% and 79%. These MED/Q genomic resources lay a foundation for future 'pan-genomic' comparisons of invasive vs. noninvasive, invasive vs. invasive, and native vs. exotic Bemisia, which, in return, will open up new avenues of investigation into whitefly biology, evolution, and management.

Transcriptome analysis and development of simple sequence repeat (SSR) markers in Zingiber striolatum Diels

Illumina-based paired-end sequencing technology was used for the high-throughput transcriptome sequencing of combined Zingiber striolatum Diels tissues (i.e., root, stem, leaf, flower, and fruit tissues). More than 130 million sequencing reads were generated, and a de novo assembly yielded 287,959 contigs and 112,107 unigenes with an average length of 1029 and 28,891 bp, respectively. Similarity searches with known sequences led to the identification of 51,804 (46.21%) genes. Of the annotated unigenes, 6867 and 51,987 were assigned to Gene Ontology and Clusters of Orthologous Groups categories, respectively. Additionally, 8384 simple sequence repeats (SSRs) were identified as potential molecular markers in the unigenes. Thirty pairs of polymerase chain reaction primers were designed and used to validate the unigenes and assess the associated genomic polymorphism. The PCR amplification products for 25 primer pairs were of the expected size. These primers may represent usable molecular markers. The thousands of SSR markers identified in the present study may be useful for analyses of genetic diversity, genetic linkage mapping, and the identification and improvement of varieties during the breeding of Z. striolatum Diels. The unigene sequences and SSR markers described herein may serve as valuable resources for future investigations of Z. striolatum Diels.

Blood transcriptomics of captive forest musk deer (Moschus berezovskii) and possible associations with the immune response to abscesses

Forest musk deer (Moschus berezovskii; FMD) are both economically valuable and highly endangered. A problem for FMD captive breeding programs has been the susceptibility of FMD to abscesses. To investigate the mechanisms of abscess development in FMD, the blood transcriptomes of three purulent and three healthy individuals were generated. A total of ~39.68 Gb bases were generated using Illumina HiSeq 4000 sequencing technology and 77,752 unigenes were identified after assembling. All the unigenes were annotated, with 63,531 (81.71%) mapping to at least one database. Based on these functional annotations, 45,798 coding sequences (CDS) were detected, along with 12,697 simple sequence repeats (SSRs) and 65,536 single nucleotide polymorphisms (SNPs). A total of 113 unigenes were found to be differentially expressed between healthy and purulent individuals. Functional annotation indicated that most of these differentially expressed genes were involved in the regulation of immune system processes, particularly those associated with parasitic and bacterial infection pathways.

Development of genome- and transcriptome-derived microsatellites in related species of snapping shrimps with highly duplicated genomes

Molecular markers are powerful tools for studying patterns of relatedness and parentage within populations and for making inferences about social evolution. However, the development of molecular markers for simultaneous study of multiple species presents challenges, particularly when species exhibit genome duplication or polyploidy. We developed microsatellite markers for Synalpheus shrimp, a genus in which species exhibit not only great variation in social organization, but also interspecific variation in genome size and partial genome duplication. From the four primary clades within Synalpheus, we identified microsatellites in the genomes of four species and in the consensus transcriptome of two species. Ultimately, we designed and tested primers for 143 microsatellite markers across 25 species. Although the majority of markers were disomic, many markers were polysomic for certain species. Surprisingly, we found no relationship between genome size and the number of polysomic markers. As expected, markers developed for a given species amplified better for closely related species than for more distant relatives. Finally, the markers developed from the transcriptome were more likely to work successfully and to be disomic than those developed from the genome, suggesting that consensus transcriptomes are likely to be conserved across species. Our findings suggest that the transcriptome, particularly consensus sequences from multiple species, can be a valuable source of molecular markers for taxa with complex, duplicated genomes.

Genome microsatellite diversity within the Apicomplexa phylum

The Apicomplexa phylum groups include unicellular and obligate intracellular protozoan parasites with an apical complex used for attachment and invasion to host cells. In this study, we analyze single sequence repeats (SSRs) in the whole genome of 20 apicomplexan organisms that represent four different lineages within the phylum. Only perfect SSRs with at least 12 nucleotides and composed of 2-6 mers were included. To better understand the association of SSR types with the genomic regions, the SSRs were classified accordingly with the genomic location into exon, intron and intergenic categories. Our results showed heterogeneous SSRs density within the studied genomes. However, the most frequent SSRs types were di- and tri-nucleotide repeats. The former was associated with intergenic regions, while the latter was associated with exon regions.

Transcriptome-Derived Tetranucleotide Microsatellites and Their Associated Genes from the Giant Panda (Ailuropoda melanoleuca)

Recently, an increasing number of microsatellites or simple sequence repeats (SSRs) have been found and characterized from transcriptomes. Such SSRs can be employed as putative functional markers to easily tag corresponding genes, which play an important role in biomedical studies and genetic analysis. However, the transcriptome-derived SSRs for giant panda (Ailuropoda melanoleuca) are not yet available. In this work, we identified and characterized 20 tetranucleotide microsatellite loci from a transcript database generated from the blood of giant panda. Furthermore, we assigned their predicted transcriptome locations: 16 loci were assigned to untranslated regions (UTRs) and 4 loci were assigned to coding regions (CDSs). Gene identities of 14 transcripts contained corresponding microsatellites were determined, which provide useful information to study the potential contribution of SSRs to gene regulation in giant panda. The polymorphic information content (PIC) values ranged from 0.293 to 0.789 with an average of 0.603 for the 16 UTRs-derived SSRs. Interestingly, 4 CDS-derived microsatellites developed in our study were also polymorphic, and the instability of these 4 CDS-derived SSRs was further validated by re-genotyping and sequencing. The genes containing these 4 CDS-derived SSRs were embedded with various types of repeat motifs. The interaction of all the length-changing SSRs might provide a way against coding region frameshift caused by microsatellite instability. We hope these newly gene-associated biomarkers will pave the way for genetic and biomedical studies for giant panda in the future. In sum, this set of transcriptome-derived markers complements the genetic resources available for giant panda.

Distribution and population genetic variation of cryptic species of the Alpine mayfly Baetis alpinus (Ephemeroptera: Baetidae) in the Central Alps

BACKGROUND

Many species contain evolutionarily distinct groups that are genetically highly differentiated but morphologically difficult to distinguish (i.e., cryptic species). The presence of cryptic species poses significant challenges for the accurate assessment of biodiversity and, if unrecognized, may lead to erroneous inferences in many fields of biological research and conservation.

RESULTS

We tested for cryptic genetic variation within the broadly distributed alpine mayfly Baetis alpinus across several major European drainages in the central Alps. Bayesian clustering and multivariate analyses of nuclear microsatellite loci, combined with phylogenetic analyses of mitochondrial DNA, were used to assess population genetic structure and diversity. We identified two genetically highly differentiated lineages (A and B) that had no obvious differences in regional distribution patterns, and occurred in local sympatry. Furthermore, the two lineages differed in relative abundance, overall levels of genetic diversity as well as patterns of population structure: lineage A was abundant, widely distributed and had a higher level of genetic variation, whereas lineage B was less abundant, more prevalent in spring-fed tributaries than glacier-fed streams and restricted to high elevations. Subsequent morphological analyses revealed that traits previously acknowledged as intraspecific variation of B. alpinus in fact segregated these two lineages.

CONCLUSIONS

Taken together, our findings indicate that even common and apparently ecologically well-studied species may consist of reproductively isolated units, with distinct evolutionary histories and likely different ecology and evolutionary potential. These findings emphasize the need to investigate hidden diversity even in well-known species to allow for appropriate assessment of biological diversity and conservation measures.

A comprehensive characterization of simple sequence repeats in pepper genomes provides valuable resources for marker development in Capsicum

The sequences of the full set of pepper genomes including nuclear, mitochondrial and chloroplast are now available for use. However, the overall of simple sequence repeats (SSR) distribution in these genomes and their practical implications for molecular marker development in Capsicum have not yet been described. Here, an average of 868,047.50, 45.50 and 30.00 SSR loci were identified in the nuclear, mitochondrial and chloroplast genomes of pepper, respectively. Subsequently, systematic comparisons of various species, genome types, motif lengths, repeat numbers and classified types were executed and discussed. In addition, a local database composed of 113,500 in silico unique SSR primer pairs was built using a homemade bioinformatics workflow. As a pilot study, 65 polymorphic markers were validated among a wide collection of 21 Capsicum genotypes with allele number and polymorphic information content value per marker raging from 2 to 6 and 0.05 to 0.64, respectively. Finally, a comparison of the clustering results with those of a previous study indicated the usability of the newly developed SSR markers. In summary, this first report on the comprehensive characterization of SSR motifs in pepper genomes and the very large set of SSR primer pairs will benefit various genetic studies in Capsicum.

Distribution, function and evolution characterization of microsatellite in Sargassum thunbergii (Fucales, Phaeophyta) transcriptome and their application in marker development

Using transcriptome data to mine microsatellite and develop markers has growingly become prevalent. However, characterizing the possible function of microsatellite is relatively rare. In this study, we explored microsatellites in the transcriptome of the brown alga Sargassum thunbergii and characterized the frequencies, distribution, function and evolution, and developed primers to validate these microsatellites. Our results showed that Tri-nucleotide is the most abundant, followed by di- and mono-nucleotide. The length of microsatellite was significantly affected by the repeat motif size. The density of microsatellite in the CDS region is significantly lower than that in the UTR region. The annotation of the transcripts containing microsatellite showed that 573 transcripts have GO terms and can be categorized into 42 groups. Pathways enrichment showed that microsatellites were significantly overrepresented in the genes involved in pathways such as Ubiquitin mediated proteolysis, RNA degradation, Spliceosome, etc. Primers flanking 961 microsatellite loci were designed, and among the 30 pairs of primer selected randomly for availability test, 23 were proved to be efficient. These findings provided new insight into the function and evolution of microsatellite in transcriptome, and the identified microsatellite loci within the annotated gene will be useful for developing functional markers in S. thunbergii.

Genetic differentiation of the regional Plutella xylostella populations across the Taiwan Strait based on identification of microsatellite markers

Movement of individuals through events, such as storms or crop transportation, may affect survival and distribution of insect pests, as well as population genetic structure at a regional scale. Understanding what factors contribute to gene flow in pest populations remains very important for sustainable pest management. The diamondback moth (Plutella xylostella) is an insect pest well known for its capacity of moving over short to long distances. Here, we used newly isolated microsatellite markers to analyze the genetic structure of nine populations across the Taiwan Strait of China (Taiwan and Fujian). A total of 12,152 simple sequence repeats (SSRs) were initially identified from the P. xylostella transcriptome (~94 Mb), with an average of 129 SSRs per Mb. Nine SSRs were validated to be polymorphic markers, and eight were used for this population genetic study. Our results showed that the P. xylostella populations could be divided into distinct two clusters, which is likely due to the year-round airflows in this region. A pattern of isolation by distance among the local populations within Fujian was found, and may be related to vegetable transportation. Considering the complexity of the P. xylostella population genetic structure from local and regional to global levels, we propose that developing ecologically sound strategies for managing this pest will require knowledge of the link between behavioral and population ecology and its genetic structure.

Sequence variation of Bemisia tabaci Chemosensory Protein 2 in cryptic species B and Q: New DNA markers for whitefly recognition

Bemisia tabaci Gennadius biotypes B and Q are two of the most important worldwide agricultural insect pests. Genomic sequences of Type-2 B. tabaci chemosensory protein (BtabCSP2) were cloned and sequenced in B and Q biotypes, revealing key biotype-specific variations in the intron sequence. A Q260 sequence was found specifically in Q-BtabCSP2 and Cucumis melo LN692399, suggesting ancestral horizontal transfer of gene between the insect and the plant through bacteria. A cleaved amplified polymorphic sequences (CAPS) method was then developed to differentiate B and Q based on the sequence variation in exon of BtabCSP2 gene. The performances of CSP2-based CAPS for whitefly recognition were assessed using B. tabaci field collections from Shandong Province (P.R. China). Our SacII based CAPS method led to the same result compared to mitochondrial cytochrome oxidase-based CAPS method in the field collections. We therefore propose an explanation for CSP origin and a new rapid simple molecular method based on genomic DNA and chemosensory gene to differentiate accurately the B and Q whiteflies of the Bemisia complex around the world.