Developing population identification tool based on polymorphism of rDNA for traditional Chinese medicine: Artemisia annua L

BACKGROUND

Artemisia annua, a well-known traditional Chinese medicine, is the main source for production of artemisinin, an anti-malaria drug. A. annua is distributed globally, with great diversity of morphological characteristics and artemisinin contents. Diverse traits among A. annua populations impeded the stable production of artemisinin, which needs an efficient tool to identify strains and assess population genetic homogeneity.

PURPOSE

In this study, ribosomal DNA (rDNA), were characterized for A. annua for strains identification and population genetic homogeneity assessment.

METHODS

The ribosomal RNA (rRNA) genes were identified using cmscan and assembled using rDNA unit of LQ-9 as a reference. rDNA among Asteraceae species were compared performing with 45S rDNA. The rDNA copy number was calculated based on sequencing depth. The polymorphisms of rDNA sequences were identified with bam-readcount, and confirmed by Sanger sequencing and restriction enzyme experiment. The ITS2 amplicon sequencing was used to verify the stability of ITS2 haplotype analysis.

RESULTS

Different from other Asteraceae species, 45S and 5S linked-type rDNA was only found in Artemisia genus. Rich polymorphisms of copy number and sequence of rDNA were identified in A. annua population. The haplotype composition of internal transcribed spacer 2 (ITS2) region which had moderate sequence polymorphism and relative short size was significantly different among A. annua strains. A population discrimination method was developed based on ITS2 haplotype analysis with high-throughput sequencing.

CONCLUSION

This study provides comprehensive characteristics of rDNA and suggests that ITS2 haplotype analysis is ideal tool for A. annua strain identification and population genetic homogeneity assessment.

Macroevolutionary patterns in intragenomic rDNA variability among planktonic foraminifera

Ribosomal intragenomic variability in prokaryotes and eukaryotes is a genomic feature commonly studied for its inflationary impact on molecular diversity assessments. However, the evolutionary mechanisms and distribution of this phenomenon within a microbial group are rarely explored. Here, we investigate the intragenomic variability in 33 species of planktonic foraminifera, calcifying marine protists, by inspecting 2,403 partial SSU sequences obtained from single-cell clone libraries. Our analyses show that polymorphisms are common among planktonic foraminifera species, but the number of polymorphic sites significantly differs among clades. With our molecular simulations, we could assess that most of these mutations are located in paired regions that do not affect the secondary structure of the SSU fragment. Finally, by mapping the number of polymorphic sites on the phylogeny of the clades, we were able to discuss the evolution and potential sources of intragenomic variability in planktonic foraminifera, linking this trait to the distinctive nuclear and genomic dynamics of this microbial group.

Allelic diversity and forensic estimations of the Beijing Hans: Comparative data on sequence-based and length-based STRs

Short tandem repeat (STR) profiling is routinely used in forensic genetics. At present, STR analysis is mainly performed by capillary electrophoresis (CE). However, due to limitations associated with the CE method, STR genotyping has been limited to length polymorphisms only. Because next generation sequencing (NGS) is capable of providing full resolution STR data at the sequence variation level, the individual identification capability of forensic STR loci could be significantly improved. Here we present sequence-based STR data for the Beijing Han population in which 291 individuals were screened for 23 commonly used forensic STRs using the SeqTypeR24 CASE kit on an Ion PGM platform. In total, 234 length-based alleles and 356 sequence-based alleles, which included 22 novel core repeat sequences, were observed. The sequence-based matching probability and power of discrimination were superior to the length-based numbers for 16 loci bearing micro-variant alleles. Combined matching probability reached 8.2 × 10^-29 for 23 STR loci at the sequence level. This was two orders of magnitude higher than the parameters at length level and provides a data base for sequence-based STR casework applications.

Non-synonymous amino acid alterations in PfEBA-175 modulate the merozoite ligand's ability to interact with host's Glycophorin A receptor

The pathological outcome of malaria due to Plasmodium falciparum infection depends largely on erythrocyte invasion by blood-stage merozoites which employ a cascade of interactions occurring between parasite ligands and RBC receptors. In a previous study exploring the genetic diversity of region-II of PfEBA-175, a ligand that plays a crucial part in parasite's RBC entry through Glycophorin A (GPA) receptor, we demonstrated that F2 domain of region-II underwent positive selection in Indian P. falciparum population through the accumulation of non-synonymous polymorphisms. Here, we examine the functional impact of two highly prevalent non-synonymous alterations in F2, namely Q584E & E592A, using a battery of molecular, biophysical and in-silico techniques. Application of circular dichroism, FTIR, fluorescence spectroscopy reveals that secondary and three-dimensional folding of recombinant-F2 protein carrying 584E and 592A residues (F2-Mut) differs significantly from that carrying 584Q and 592E (F2-3D7). A comparison of spectroscopic and thermodynamic parameters shows that F2-Mut is capable of forming a complex with GPA with higher efficiency compared to F2-3D7. In silico docking predicts both artemisinin and artesunate possess the capacity of slipping into the GPA binding crevices of PfEBA-175 and disrupt PfEBA-GPA association. However, the estimated affinity of artesunate towards PfEBA-175 with 584E and 592A residues is higher than that of artemisinin. Thermodynamic parameters computed using isotherms are concordant with this in-silico prediction. Together, our data suggest that the presence of amino acid alterations in F2 provide structural and functional stability favoring PfEBA-GPA interaction and artesunate can efficiently disrupt the interaction between GPA and PfEBA-175 even carrying altered amino acid residues. The present study alerts the malaria research community by presenting evidence that the parasite is gaining evolutionary fitness by cultivating genetic alterations in many of its proteins.

An assessment on the intrapopulational and intraindividual genetic diversity in LSU rDNA in the harmful algal blooms-forming dinoflagellate Margalefidinium (= Cochlodinium) fulvescens based on clonal cultures and bloom samples from Jiaozhou Bay, China

Large subunit ribosomal DNA (LSU rDNA) sequences have been increasingly used to infer the phylogeny and species identity of organisms, a few previous studies, however, have observed high intraspecific and even intraindividual variability in LSU rDNA in some dinoflagellate species due to, assumably, large copy numbers of rDNA in dinoflagellates. Since the copy number of LSU rDNA varies tremendously among dinoflagellate species, the intraspecific and intraindividual diversity for a species of particular interest thus needs to be investigated individually. As a toxic and HABs-forming dinoflagellate, Margalefidinium (= Cochlodinium) fulvescens has been observed to approach blooming density in Jiaozhou Bay, China since 2015 after numerous blooms having been reported from other countries. In trying to identify the source of this newly observed HABs-forming species in China by sequencing the LSU rDNA for both field samples and clonal cultures, we noticed and thus further investigated high intrapopulational and intraindividual genetic diversities of the dinoflagellate. The D1-D6 region of the LSU rDNA (1,435 bases) was amplified from 7 field samples (pooled cells) and 11 clonal cultures, cloned, sequenced, and analyzed phylogenetically for 2,341 sequences obtained. All the numbers of sequences obtained from each clonal culture were far less than the estimated rDNA copy number in M. fulvescens. In the clone library, only one unique sequence was contained in all samples as the most dominant sequence. We found high intrapopulational and intraindividual genetic diversity in M. fulvescens as reflected in the number of polymorphic sites and unique sequences in the clone library for different field samples and clonal cultures in comparison to other species. The mean number of nucleotide differences of each sequence from different field samples and clonal cultures were 6.43 and 4.42 bases, respectively, with the highest being 132 bases, nearly 10%. The sequences with highest variability may be easily annotated as different species if they were obtained from environmental genomic studies because sequence-based species identification in meta-barcoding studies often use "97% identity" threshold. Based on that the mean and overall intrapopulational genetic diversity calculated for 7 field samples was equivalent to the mean and overall intraindividual variability for 11 clonal cultures in indices of genetic diversity, together with the result of AMOVA analysis, we infer that the variability within individual cells (i.e. variability among LSU rDNA polymorphic copies) caused both the intraindividual and intrapopulational genetic diversities observed in the M. fulvescens population, and a higher interpopulational diversity may exist among different geographic populations. The results provide an insightful basis for such a comprehensive interpopulational comparison and important implications for identifying species and establishing new taxa based on the similarity comparison to reference sequences deposited in databases.

Characterization of sequence variation at 30 autosomal STRs in Chinese Han and Tibetan populations

Massively parallel sequencing (MPS) technologies have the ability to reveal sequence variations within STR alleles as well as their nominal allele lengths, which have traditionally been detected by CE instruments. Recently, Thermo Fisher Scientific has updated the MPS-STR panel, named the Precision ID GlobalFiler next-generation sequencing (NGS) STR Panel version 2, with primers redesigned to add two pentanucleotide tandem repeat loci and profile interpretation supported by the Converge software. Using the Ion Chef System, the Ion S5XL System, and the Converge software, genetic variations were characterized within STR repeat and flanking regions of 30 autosomal STR markers in 115 unrelated individuals from two Chinese population groups (58 Tibetans and 57 Hans). Nineteen STRs demonstrated a relative increase in diversity with the variant sequence alleles compared with those of traditional nominal length alleles. In total, 390 alleles were identified by their sequences compared with 258 alleles that were identified by length. Of these 92 sequence variants found within the STR repeat regions, 40 variants were located in STR flanking regions. Additionally, the agreement of the results with CE data was evaluated, as was the ability of this new MPS panel to analyze case-type (11 samples) and artificially degraded samples (seven samples in triplicate). The results generated from this study illustrate that extensive sequence variation exists in commonly used STR markers in the selected population samples and indicate that this NGS STR panel has the potential to be used as an effective tool for human forensics.

Complete genome of DENV2 isolated from mosquitoes in Mexico

Dengue virus is the most prevalent arbovirus in Mexico, and although the diversity of this virus has been studied, the vast majority of sequences have been derived from viruses isolated from the human host. In this work, we aimed to sequence and to analyze DENVs derived from wild mosquitoes captured in Acapulco Guerrero, Mexico. We succeeded in determining three full genome sequences of such viruses and were able to compare them with other reported sequences from human and mosquito-derived DENVs. We found 15 nonsynonymous and 88 synonymous substitutions that were present more frequently in mosquito viruses than what would be expected by chance, although the limited number of genomes reported so far puts a constraint on the conclusions that can be derived from these analyses. Also, given the high depth of coverage attained in one of the genomes a variant analysis was carried out, finding 68 polymorphic sites in this genome. Interestingly, six of them corresponded to SNV that were detected as potentially differential between mosquitoes and humans, indicating that a that at least some positions may be maintained as polymorphic, which may facilitate host transmission.

Novel Y-chromosome Short Tandem Repeat Variants Detected Through the Use of Massively Parallel Sequencing

Massively parallel sequencing (MPS) technology is capable of determining the sizes of short tandem repeat (STR) alleles as well as their individual nucleotide sequences. Thus, single nucleotide polymorphisms (SNPs) within the repeat regions of STRs and variations in the pattern of repeat units in a given repeat motif can be used to differentiate alleles of the same length. In this study, MPS was used to sequence 28 forensically-relevant Y-chromosome STRs in a set of 41 DNA samples from the 3 major U.S. population groups (African Americans, Caucasians, and Hispanics). The resulting sequence data, which were analyzed with STRait Razor v2.0, revealed 37 unique allele sequence variants that have not been previously reported. Of these, 19 sequences were variations of documented sequences resulting from the presence of intra-repeat SNPs or alternative repeat unit patterns. Despite a limited sampling, two of the most frequently-observed variants were found only in African American samples. The remaining 18 variants represented allele sequences for which there were no published data with which to compare. These findings illustrate the great potential of MPS with regard to increasing the resolving power of STR typing and emphasize the need for sample population characterization of STR alleles.

A high-throughput pipeline for detecting locus-specific polymorphism in hexaploid wheat (Triticum aestivum L.)

BACKGROUND

Bread wheat (Triticum aestivum L., 2n = 6x = 42) is an allohexaploid with a huge genome. Due to the presence of extensive homoeologs and paralogs, generating locus-specific sequences can be challenging, especially when a large number of sequences are required. Traditional methods of generating locus-specific sequences are rather strenuous and time-consuming if large numbers of sequences are to be handled.

RESULTS

To improve the efficiency of isolating sequences for targeted loci, a time-saving and high-throughput pipeline integrating orthologous sequence alignment, genomic sequence retrieving, and multiple sequence alignment was developed. This pipeline was successfully employed in retrieving and aligning homoeologous sequences and 83% of the primers designed based on the pipeline successfully amplified fragments from the targeted subgenomes.

CONCLUSIONS

The high-throughput pipeline developed in this study makes it feasible to efficiently identify locus-specific sequences for large numbers of sequences. It could find applications in all research projects where locus-specific sequences are required. In addition to generating locus-specific markers, the pipeline was also used in our laboratory to identify differentially expressed genes among the three subgenomes of bread wheat. Importantly, the pipeline is not only valuable for research in wheat but should also be applicable to other allopolyploid species.

Effect of temperature on symptom expression and sequence polymorphism of grapevine yellow speckle viroid 1 in grapevine

Symptom expression of yellow speckle disease was studied in a row of 32 individual Vitis vinifera cv. Chardonnay vines in the warmer summer of 2009/2010 as compared with the cooler summer of 2011/2012 in South Australia. RT-PCR analysis showed that all these vines were positive for hop stunt viroid, grapevine yellow speckle viroid 1 (GYSVd-1) and grapevine rupestris stem pitting-associated virus. Four vines named Vines 1, 8, 11 and 15 were selected for further analyses. Vines 1 and 8 had never shown yellow speckle (YS) symptoms, Vine 11 had always been symptomatic, and Vine 15 showed YS symptoms only in the summer of 2009/2010, but not in 2011/2012. Analysis of partial nucleotide (nt) sequence of GYSVd-1 from these vines showed two major sequence polymorphisms in the pathogenicity domain coinciding with the YS symptoms and the prevailing temperature in each season. One group designated UA group had a uridine (U) at position 309 and an adenine (A) at position 311, while another group designated AU group had the other way around: an "A" at position 309 and a "U" at position 311. The AU group had never been reported before. The AU group was a minor variant in the GYSVd-1 population and not present in symptomatic Vine 11. In contrast, the UA group was dominant and present in all the vine samples. Surprisingly, all the asymptomatic vines, but not symptomatic vines, had the signature of the AU group. Whether the AU group is associated with the YS symptom expression is interesting. Our result provides a new insight into the sequence variability of viroid-inducing symptoms during two significantly different growing seasons.