Sequence Analysis of the Complete Mitochondrial Genome of a Medicinal Plant, Vitex rotundifolia Linnaeus f. (Lamiales: Lamiaceae).

In the present study, we depicted the complete mitochondrial genome of a valuable medicinal plant, Vitex rotundifolia. The mitochondrial genome of V. rotundifolia, mapped as a circular molecule, spanned 380,980 bp in length and had a GC content of 45.54%. The complete genome contained 38 protein-coding genes, 19 transfer RNAs (tRNAs), and 3 ribosomal RNAs (rRNAs). We found that there were only 38.73% (147.54 kb), 36.28% (138.23 kb), and 52.22% (198.96 kb) of the homologous sequences in the mitochondrial genome of V. rotundifolia, as compared with the mitochondrial genomes of Scutellaria tsinyunensis, Boea hygrometrica, and Erythranthe lutea, respectively. A multipartite structure mediated by the homologous recombinations of the three direct repeats was found in the V. rotundifolia mitochondrial genome. The phylogenetic tree was built based on 10 species of Lamiales, using the maximum likelihood method. Moreover, this phylogenetic analysis is the first to present the evolutionary relationship of V. rotundifolia with the other species in Lamiales, based on the complete mitochondrial genome.

Clinical Pharmacogenetics Implementation Consortium Guideline for the Use of Aminoglycosides Based on MT-RNR1 Genotype.

Aminoglycosides are widely used antibiotics with notable side effects, such as nephrotoxicity, vestibulotoxicity, and sensorineural hearing loss (cochleotoxicity). MT-RNR1 is a gene that encodes the 12s rRNA subunit and is the mitochondrial homologue of the prokaryotic 16s rRNA. Some MT-RNR1 variants (i.e., m.1095T>C; m.1494C>T; m.1555A>G) more closely resemble the bacterial 16s rRNA subunit and result in increased risk of aminoglycoside-induced hearing loss. Use of aminoglycosides should be avoided in individuals with an MT-RNR1 variant associated with an increased risk of aminoglycoside-induced hearing loss unless the high risk of permanent hearing loss is outweighed by the severity of infection and safe or effective alternative therapies are not available. We summarize evidence from the literature supporting this association and provide therapeutic recommendations for the use of aminoglycosides based on MT-RNR1 genotype (updates at https://cpicpgx.org/guidelines/ and www.pharmgkb.org).

FISH-Based Karyotype Analyses of Four Dracaena Species.

The genus Dracaena is the main source of dragon's blood, which is a plant resin and has been used as traditional medicine since ancient times in different civilizations. However, the chromosome numbers and karyotypes present in this genus remain poorly understood. In this study, fluorescence in situ hybridization (FISH) using oligonucleotide probes for ribosomal DNAs (5S and 45S rDNA) and telomeric repeats (TTTAGGG)3 was applied to analyze 4 related species: Dracaena terniflora Roxb., Dracaena cambodiana Pierre ex Gagnep., Aizong (Dracaena sp.), and Dracaena cochinchinensis (Lour.) S.C. Chen. In all 4 species, both 5S and 45S rDNA showed hybridization signals in the paracentromeric region of a pair of chromosomes; the sizes of the 45S rDNA signals were larger than those of the 5S rDNA. Importantly, the telomeric repeat signals were located in the telomeric regions of almost all chromosomes. The results indicated that the chromosome number of all 4 Dracaena species is 2n = 40, and the lengths of the mitotic metaphase chromosomes range from 0.99 to 2.98 µm. Our results provide useful cytogenetic information, which will be beneficial to future studies in genome structure of the genus Dracaena.

Phylogenetic and morphological significance of an overlooked flying squirrel (Pteromyini, Rodentia) from the eastern Himalayas with the description of a new genus.

The flying squirrels (Pteromyini, Rodentia) are the most diverse and widely distributed group of gliding mammals. Taxonomic boundaries and relationships within flying squirrels remain an area of active research in mammalogy. The discovery of new specimens of Pteromys ( Hylopetes) leonardi Thomas, 1921 previously considered a synonym of Hylopetes alboniger, in Yunnan Province, China allowed a morphological and genetic reassessment of the status of this taxon. Phylogenetic reconstruction was implemented using sequences of two mitochondrial (12S ribosomal DNA and 16S ribosomal DNA) and one nuclear (interphotoreceptor retinoid-binding protein) gene fragments. Morphological assessments involved examinations of features preserved on skins, skulls, and penises of museum specimens, supplemented with principal component analysis of craniometric data. Together these assessments revealed that this taxon should be recognized not only as a distinct species, and should also be placed within a new genus, described here as Priapomys.

Estimating included animal species in mixed crude drugs derived from animals using massively parallel sequencing.

We developed a method that can detect each animal species of origin for crude drugs derived from multiple animal species based on massively parallel sequencing analysis of mitochondrial genes. The crude drugs derived from animals investigated in this study were Cervi Parvum Cornu and Trogopterorum feces, which are derived from a mix of different animal species, two chopped cicada sloughs, and two commercial Kampo drugs. The mitochondrial 12S rRNA, 16S rRNA, and cytochrome oxidase subunit I gene regions were amplified and sequenced using MiSeq. The ratios of haplotype to total number of sequences reads were calculated after sequence extraction and trimming. Haplotypes that exceeded the threshold were defined as positive haplotypes, which were compared with all available sequences using BLAST. In the Cervi Parvum Cornu and Trogopterorum feces samples, the haplotype ratios corresponded roughly to the mixture ratios, although there was a slight difference from mixture ratios depending on the gene examined. This method could also roughly estimate the compositions of chopped cicada sloughs and Kampo drugs. This analysis, whereby the sequences of several genes are elucidated, is better for identifying the included animal species. This method should be useful for quality control of crude drugs and Kampo drugs.

Cytomolecular analysis of mutants, breeding lines, and varieties of camelina (Camelina sativa L. Crantz).

Camelina sativa L. Crantz (Brassicaceae family), known as camelina, has gained new attention as a re-emerging oil seed crop. With a unique seed oil profile, with the majority of the fatty acids consisting of linolenic (C18:3), oleic (C18:1), linoleic (C18:2), and eicosenoic (C20:1), camelina oil is reported to be useful as a food oil and biofuel. However, there are still many unknown factors about the structure and genetic variability of this crop. Chromosomal localization of ribosomal DNA was performed using fluorescence in situ hybridization (FISH) with 5S rDNA and 25S rDNA sequences as molecular probes on mitotic chromosomes of enzymatically digested root-tip meristematic cells. Here, we present for the first time a comparative analysis of selected genotypes (cultivars, breeding lines and mutants) of C. sativa with the use of cytogenetic techniques. The main aim of the study was to determine the intraspecific and interspecific polymorphisms in the structure of chromosomes of selected accessions using conserved 5S and 25S rDNA repetitive sequences as molecular probes. The results were compared with C. microcarpa (closely related to C. sativa) rDNA gene loci distribution. The presence of minor rDNA sites was discussed and compared with other Brassicaceae species. In addition, demonstration karyograms of C. sativa and C. microcarpa mapped with rDNA probes were prepared based on the cv. "Przybrodzka" and GE2011-02 genotype, respectively. The use of 5S and 25S rDNA probes provided an insight on the genome structure of C. sativa at the cytogenetic level and can help to understand the genome organization of this crop. The putative role of cytogenetic markers in phylogenetic analyses of camelina was discussed, as well.

Candida metrosideri pro tempore sp. nov. and Candida ohialehuae pro tempore sp. nov., two antifungal-resistant yeasts associated with Metrosideros polymorpha flowers in Hawaii.

Flowers produce an array of nutrient-rich exudates in which microbes can thrive, making them hotspots for microbial abundance and diversity. During a diversity study of yeasts inhabiting the flowers of Metrosideros polymorpha (Myrtaceae) in the Hawai'i Volcanoes National Park (HI, USA), five isolates were found to represent two novel species. Morphological and physiological characterization, and sequence analysis of the small subunit ribosomal RNA (rRNA) genes, the D1/D2 domains of the large subunit rRNA genes, the internal transcribed spacer (ITS) regions, and the genes encoding the largest and second largest subunits of the RNA polymerase II (RPB1 and RPB2, respectively), classified both species in the family Metschnikowiaceae, and we propose the names Candida metrosideri pro tempore sp. nov. (JK22T = CBS 16091 = MUCL 57821) and Candida ohialehuae pro tempore sp. nov. (JK58.2T = CBS 16092 = MUCL 57822) for such new taxa. Both novel Candida species form a well-supported subclade in the Metschnikowiaceae containing species associated with insects, flowers, and a few species of clinical importance. The ascosporic state of the novel species was not observed. The two novel yeast species showed elevated minimum inhibitory concentrations to the antifungal drug amphotericin B (>4 µg/mL). The ecology and phylogenetic relationships of C. metrosideri and C. ohialehuae are also discussed.

The complete chloroplast genome of Gleditsia sinensis and Gleditsia japonica: genome organization, comparative analysis, and development of taxon specific DNA mini-barcodes.

Chloroplast genomes have been widely considered an informative and valuable resource for molecular marker development and phylogenetic reconstruction in plant species. This study evaluated the complete chloroplast genomes of the traditional Chinese medicine Gleditsia sinensis and G. japonica, an adulterant of the former. The complete chloroplast genomes of G. sinensis and G. japonica were found to be of sizes 163,175 bp and 162,391 bp, respectively. A total of 111 genes were identified in each chloroplast genome, including 77 coding sequences, 30 tRNA, and 4 rRNA genes. Comparative analysis demonstrated that the chloroplast genomes of these two species were highly conserved in genome size, GC contents, and gene organization. Additionally, nucleotide diversity analysis of the two chloroplast genomes revealed that the two short regions of ycf1b were highly diverse, and could be treated as mini-barcode candidate regions. The mini-barcode of primers ZJ818F-1038R was proven to precisely discriminate between these two species and reflect their biomass ratio accurately. Overall, the findings of our study will shed light on the genetic evolution and guide species identification of G. sinensis and G. japonica.

The complete chloroplast genome and characteristics analysis of Callistemon rigidus R.Br.

Callistemon rigidus R.Br. one of the traditional Chinese medicinal plants, is acrid-flavored and mild-natured, with the prominent effects reducing swelling, resolving phlegm, and dispelling rheumatism. Clinically, it has been commonly used to treat cold, cough and asthma, pain and swelling from impact injuries, eczema, rheumatic arthralgia. The chloroplast genome study on Callistemon rigidus R.Br. is a few seen. This study demonstrates the data collected from the assembly and annotation of the chloroplast (cp) genome of Callistemon rigidus R.Br., followed by furthers comparative analysis with the cp genomes of closely related species. C. rigidus R.Br. showed a cp genome in the size of 158, 961 bp long with 36.78% GC content, among which a pair of inverted repeats (IRs) of 26, 671 bp separated a large single-copy (LSC) region of 87, 162 bp and a small single-copy (SSC) region of 18, 457 bp. Altogether 131 genes were hosted, including 37 transfer RNAs, 8 ribosomal RNAs, and 86 protein-coding genes. 284 simple sequence repeats (SSRs) were also marked out. A comparative analysis of the genome structure and the sequence data of closely related species unveiled the conserved gene order in the IR and LSC/SSC regions, a quite constructive finding for future phylogenetic research. Overall, this study providing C. rigidus R.Br. genomic resources could positively contribute to the evolutionary study and the phylogenetic reconstruction of Myrtaceae.

Design and validation of a real-time PCR technique for assessing the level of inclusion of fungus- and yeast-based additives in feeds.

A reliable method for quantification of non-viable microbe-based nutritional and zootechnical additives introduced into feed is essential in order to ensure regulatory compliance, feed safety and product authenticity in industrial applications. In the present work, we developed a novel real-time quantitative polymerase chain reaction (qPCR) -based analysis protocol for monitoring two microbial additives in feed matrices. To evaluate the applicability of the method, pelleted wheat- and maize-based broiler chicken diets containing a non-viable phytase-producing strain of Aspergillus niger produced in solid state fermentation (150 or 300 g/t) and a non-viable selenium-enriched Saccharomyces cerevisiae (100 or 200 g/t) as model feed ingredients, were manufactured and subjected to analysis. Power analysis of the qPCR results indicated that 2 to 6 replicate feed samples were required to distinguish the product doses applied, which confirms that the microbial DNA was efficiently recovered and that potential PCR inhibitors present in the feed material were successfully removed in DNA extraction. The analysis concept described here was shown to be an accurate and sensitive tool for monitoring the inclusion levels of non-viable, unculturable microbial supplements in animal diets.

The origin identification method for crude drugs derived from arthropods and annelids using molecular biological techniques.

We evaluated whether the origins of crude drugs derived from arthropods and annelids could be identified using molecular biological techniques. DNA was extracted from 20 crude drugs prepared from different animals using a commercial kit with added phenol treatment. The target regions used to identify origin were the mitochondrial 16S ribosomal RNA (rRNA), 12S rRNA, and cytochrome oxidase subunit I (COI) gene regions. Extracted DNA was amplified by polymerase chain reaction, and then sequenced by the Sanger method. The aligned sequences were compared with all available sequences using BLAST to estimate the origins of the crude drugs. The origin of crude drugs used in this study could be estimated using this method. The COI region was the best for identifying origin among three regions examined, based on the success rate of PCR amplification and analysis. Moreover, the 12S rRNA region was also useful for origin identification, with the exception of the earthworm. However, the origin of some crude drugs could not be strictly identified due to matches to various species in all three regions. One likely cause was that the species of origin of a crude drug has not been registered in DNA databases. We found that even the same crude drug from the same pharmaceutical company had different origins by production lot or import source country. Therefore, this method is useful not only for DNA-based origin identification but also quality control of production lots.

Gene Losses and Variations in Chloroplast Genome of Parasitic Plant Macrosolen and Phylogenetic Relationships within Santalales.

Macrosolen plants are parasitic shrubs, several of which are important medicinal plants, that are used as folk medicine in some provinces of China. However, reports on Macrosolen are limited. In this study, the complete chloroplast genome sequences of Macrosolen cochinchinensis, Macrosolen tricolor and Macrosolen bibracteolatus are reported. The chloroplast genomes were sequenced by Illumina HiSeq X. The length of the chloroplast genomes ranged from 129,570 bp (M. cochinchinensis) to 126,621 bp (M. tricolor), with a total of 113 genes, including 35 tRNA, eight rRNA, 68 protein-coding genes, and two pseudogenes (ycf1 and rpl2). The simple sequence repeats are mainly comprised of A/T mononucleotide repeats. Comparative genome analyses of the three species detected the most divergent regions in the non-coding spacers. Phylogenetic analyses using maximum parsimony and maximum likelihood strongly supported the idea that Loranthaceae and Viscaceae are monophyletic clades. The data obtained in this study are beneficial for further investigations of Macrosolen in respect to evolution and molecular identification.

Occurrence of oleaginous yeast from mangrove forest in Thailand.

A total of 191 yeasts were isolated from 197 samples collected from eight estuarine mangrove forests along four different coastlines of Thailand (Andaman Sea and the East, North and West coasts of the Gulf of Thailand). Of these, 178 isolates were identified as 32 species in 16 genera of Ascomycota, 12 species in nine genera of Basidiomycota, and 13 isolates as potential new species, respectively. Mangroves located along the Andaman Sea coastline had a higher yeast diversity at the species and genera levels than those along the Gulf of Thailand. Kluyveromyces siamensis was the most frequently isolated species, whilst Candida tropicalis was the only species isolated at all eight sites. Screening isolated yeast strains belonging to genera previously reported as oleaginous yeast plus the 13 potential new species, revealed two oleaginous strains, Rhodotorula sphaerocarpa 11-14.4 and Saitozyma podzolica 11-11.3.1. Both of these strains were isolated from the same mangrove forest on the Andaman Sea coastline. They could accumulate lipid when suspended in glucose solution without any supplementation, while the fatty acid composition and oil profile of Rh. sphaerocarpa 11-14.4 and Sait. podzolica 11-11.3.1 were similar to vegetable oil and cocoa butter, respectively.

Dynamic Chloroplast Genome Rearrangement and DNA Barcoding for Three Apiaceae Species Known as the Medicinal Herb "Bang-Poong".

Three Apiaceae species Ledebouriella seseloides, Peucedanum japonicum, and Glehnia littoralis are used as Asian herbal medicines, with the confusingly similar common name "Bang-poong". We characterized the complete chloroplast (cp) genomes and 45S nuclear ribosomal DNA (45S nrDNA) sequences of two accessions for each species. The complete cp genomes of G. littoralis, L. seseloides, and P. japonicum were 147,467, 147,830, and 164,633 bp, respectively. Compared to the other species, the P. japonicum cp genome had a huge inverted repeat expansion and a segmental inversion. The 45S nrDNA cistron sequences of the three species were almost identical in size and structure. Despite the structural variation in the P. japonicum cp genome, phylogenetic analysis revealed that G. littoralis diverged 5-6 million years ago (Mya), while P. japonicum diverged from L. seseloides only 2-3 Mya. Abundant copy number variations including tandem repeats, insertion/deletions, and single nucleotide polymorphisms, were found at the interspecies level. Intraspecies-level polymorphism was also found for L. seseloides and G. littoralis. We developed nine PCR barcode markers to authenticate all three species. This study characterizes the genomic differences between L. seseloides, P. japonicum, and G. littoralis; provides a method of species identification; and sheds light on the evolutionary history of these three species.

The plastid genome and its implications in barcoding specific-chemotypes of the medicinal herb Pogostemon cablin in China.

Pogostemon cablin (Blanco) Benth. (Patchouli) is not only an important essential oil plant, but also a valuable medicinal plant in China. P. cablin in China can be divided into three cultivars (Shipai, Gaoyao, and Hainan) and two chemotypes (pogostone-type and patchoulol-type). The pogostone-type and patchoulol-type are, respectively, used for medicinals and perfumes. In this study, we sequenced and characterized the plastid genomes for all three Chinese cultivars and aimed to develop a chemotype-specific barcode for future quality control. The plastid genomes of P. cablin cultivars ranged from 152,461 to 152,462 bp in length and comprise 114 genes including 80 protein coding genes, 30 tRNA genes, and four rRNA genes. Phylogenetic analyses suggested that P. cablin cultivars clustered with the other two Pogostemon species with strong support. Although extremely conserved in P. cablin plastid genomes, 58 cpSSRs were filtered out among the three cultivars. One single variable locus, cpSSR, was discovered. The cpSSR genotypes successfully matched the chemotypes of Chinese patchouli, which was further supported by PCR-based Sanger sequences in more Chinese patchouli samples. The barcode developed in this study is thought to be a simple and reliable quality control method for Chinese P. cablin on the market.

C9orf72 Poly(PR) Dipeptide Repeats Disturb Biomolecular Phase Separation and Disrupt Nucleolar Function.

Repeat expansion in the C9orf72 gene is the most common cause of the neurodegenerative disorder amyotrophic lateral sclerosis (C9-ALS) and is linked to the unconventional translation of five dipeptide-repeat polypeptides (DPRs). The two enriched in arginine, poly(GR) and poly(PR), infiltrate liquid-like nucleoli, co-localize with the nucleolar protein nucleophosmin (NPM1), and alter the phase separation behavior of NPM1 in vitro. Here, we show that poly(PR) DPRs bind tightly to a long acidic tract within the intrinsically disordered region of NPM1, altering its phase separation with nucleolar partners to the extreme of forming large, soluble complexes that cause droplet dissolution in vitro. In cells, poly(PR) DPRs disperse NPM1 from nucleoli and entrap rRNA in static condensates in a DPR-length-dependent manner. We propose that R-rich DPR toxicity involves disrupting the role of phase separation by NPM1 in organizing ribosomal proteins and RNAs within the nucleolus.

Complete Chloroplast Genomes and Comparative Analysis of Sequences Evolution among Seven Aristolochia (Aristolochiaceae) Medicinal Species.

Aristolochiaceae, comprising about 600 species, is a unique plant family containing aristolochic acids (AAs). In this study, we sequenced seven species of Aristolochia, and retrieved eleven chloroplast (cp) genomes published for comparative genomics analysis and phylogenetic constructions. The results show that the cp genomes had a typical quadripartite structure with conserved genome arrangement and moderate divergence. The cp genomes range from 159,308 bp to 160,520 bp in length and have a similar GC content of 38.5%⁻38.9%. A total number of 113 genes were identified, including 79 protein-coding genes, 30 tRNAs and four rRNAs. Although genomic structure and size were highly conserved, the IR-SC boundary regions were variable between these seven cp genomes. The trnH-GUG genes, are one of major differences between the plastomes of the two subgenera Siphisia and Aristolochia. We analyzed the features of nucleotide substitutions, distribution of repeat sequences and simple sequences repeats (SSRs), positive selections in the cp genomes, and identified 16 hotspot regions for genomes divergence that could be utilized as potential markers for phylogeny reconstruction. Phylogenetic relationships of the family Aristolochiaceae inferred from the 18 cp genome sequences were consistent and robust, using maximum parsimony (MP), maximum likelihood (ML), and Bayesian analysis (BI) methods.

Identification and expression analysis of a microRNA cluster derived from pre-ribosomal RNA in Papaver somniferum L. and Papaver bracteatum L.

Opium poppy (Papaver somniferum L.) is one of the ancient medical crops, which produces several important alkaloids such as morphine, noscapine, sanguinarine and codeine. MicroRNAs are endogenous non-coding RNAs that play important regulatory roles in plant diverse biological processes. Many plant miRNAs are encoded as single transcriptional units, in contrast to animal miRNAs, which are often clustered. Herein, using computational approaches, a total of 22 miRNA precursors were identified, which five of them were located as a clustered in pre-ribosomal RNA. Afterward, the transcript level of the precursor and the mature of clustered miRNAs in two species of the Papaveraceae family, i.e. P. somniferum L. and P. bracteatum L, were quantified by RT-PCR. With respect to obtained results, these clustered miRNAs were expressed differentially in different tissues of these species. Moreover, using target prediction and Gene Ontology (GO)-based on functional classification indicated that these miRNAs might play crucial roles in various biological processes as well as metabolic pathways. In this study, we discovered the clustered miRNA derived from pre-rRNA, which may shed some light on the importance of miRNAs in the plant kingdom.

[Rapid identification of Gekko gecko by loop-mediated isothermal amplification based on 12S rRNA sequence].

Gekko gecko (Tokay Gecko) is a valuable traditional Chinese medicine. In this study, the loop-mediated isothermal amplification (LAMP) technique was introduced for visual rapid identification of G. gecko from adulterants. A total of sixty-five 12S rRNA sequences of fourteen species of G. gecko and its adulterants were obtained. The results showed that G. gecko could be identified from its adulterants through BLAST analysis based on 12S rRNA regions. The 12S rRNA sequences of ten batches of G. gecko were conserved. There were only two haplotypes and three variation sites in the available regions for primers design. Six specific LAMP primers were successfully designed online based on 12S rRNA sequences. The visual rapid detection of G. gecko could be achieved with the optimized conditions (64 °C for 1 h and 80 °C for 5 min). And the required minimal template concentration was 5 µg·L⁻¹ while conventional PCR with 0.5 mg·L⁻¹. Consequently, the LAMP method established from this study was rapid, specific, highly sensitive, and simple. It could be applied to detect G. gecko from its adulterants efficiently.

The Complete Chloroplast Genome Sequences of the Medicinal Plant Forsythia suspensa (Oleaceae).

Forsythia suspensa is an important medicinal plant and traditionally applied for the treatment of inflammation, pyrexia, gonorrhea, diabetes, and so on. However, there is limited sequence and genomic information available for F. suspensa. Here, we produced the complete chloroplast genomes of F. suspensa using Illumina sequencing technology. F. suspensa is the first sequenced member within the genus Forsythia (Oleaceae). The gene order and organization of the chloroplast genome of F. suspensa are similar to other Oleaceae chloroplast genomes. The F. suspensa chloroplast genome is 156,404 bp in length, exhibits a conserved quadripartite structure with a large single-copy (LSC; 87,159 bp) region, and a small single-copy (SSC; 17,811 bp) region interspersed between inverted repeat (IRa/b; 25,717 bp) regions. A total of 114 unique genes were annotated, including 80 protein-coding genes, 30 tRNA, and four rRNA. The low GC content (37.8%) and codon usage bias for A- or T-ending codons may largely affect gene codon usage. Sequence analysis identified a total of 26 forward repeats, 23 palindrome repeats with lengths >30 bp (identity > 90%), and 54 simple sequence repeats (SSRs) with an average rate of 0.35 SSRs/kb. We predicted 52 RNA editing sites in the chloroplast of F. suspensa, all for C-to-U transitions. IR expansion or contraction and the divergent regions were analyzed among several species including the reported F. suspensa in this study. Phylogenetic analysis based on whole-plastome revealed that F. suspensa, as a member of the Oleaceae family, diverged relatively early from Lamiales. This study will contribute to strengthening medicinal resource conservation, molecular phylogenetic, and genetic engineering research investigations of this species.