Comparative genomics and phylogenetic analysis of chloroplast genomes of Asian Caryodaphnopsis taxa (Lauraceae)

The genus Caryodaphnopsis, a member of the Lauraceae family, is characterized by seeds that are rich in oil, as well as highly exploitable fruits and wood. The Asian taxa within this genus exhibit complex morphological variations, posing challenges to their accurate classification and impeding their effective use and development as a resource. In this study, we sequenced the chloroplast genomes of 31 individuals representing nine Asian taxa within the Caryodaphnopsis genus. Our primary objectives were to reveal structural variations in these chloroplast genomes through comparative analyses and to infer the species' phylogenetic relationships. Our findings revealed that all chloroplast genomes had a tetrad structure, ranged in length from 148,828 to 154,946 bp, and harbored 128-131 genes. Notably, contraction of the IR region led to the absence of some genes in eight taxa. A comprehensive analysis identified 1267 long repetitive sequences and 2176 SSRs, 286 SNPs, and 135 indels across the 31 chloroplast genomes. The Ka/Ks ratio analysis indicated potential positive selection on the matK, rpl22, and rpoC2 genes. Furthermore, we identified six variable regions as promising barcode regions. Phylogenetic analysis grouped the nine Asian taxa into six branches, with C. henryi forming the basal group from which three distinct complexes emerged. This study contributes significantly to the current understanding of the evolutionary dynamics and phylogenetic relationships within the genus Caryodaphnopsis. Furthermore, the identified molecular markers hold potential for molecular barcoding applications in population genetics, providing valuable tools for future research and conservation efforts within this diverse genus.

Chloroplast genome data of five Amygdalus species: Clarifying genome structure and phylogenetic relationships

Amygdalus species have considerable ecological and economic value, however, the phylogenetic relationships among Amygdalus remain controversy. In this study, we sequenced and assembled the chloroplast (cp) genomes of five Amygdalus species: Prunus communis, P. mongolica, P. pedunculata, P. triloba, and P. mira. We then conducted comparative genomic analyses and constructed their phylogenetic relationships. The genome length ranged from 157,870 to 158,451 bp, and 131 genes were annotated (86 protein-coding genes, 37 tRNAs, and 8 rRNAs). Additionally, 49-57 simple sequence repeats were detected, with most in the large single-copy region and with AT base preferences. Comparative genomic analyses revealed high similarities in structure, order, and gene content. However, we identified four highly divergent sequences: trnR-UCU-atpA, nbdhC-trnV-UAC, ycf4-cemA, and rpl32-trnL-UAG. The phylogenomic relationship analysis suggested that the Amygdalus species were grouped together, in which P. pedunculata, P. triloba, and Prunus tangutica were categorized into a branch, P. mongolica and Prunus davidiana were clustered a branch. This study provides an improved understanding of the genetic relationships among the Amygdalus and provides a basis for the development and utilization of Amygdalus resources.

Characterization of Firmiana danxiaensis plastomes and comparative analysis of Firmiana: insight into its phylogeny and evolution

BACKGROUND

Firmiana danxiaensis is a critically endangered and ecologically important tree currently only found in four locations in Danxia or Karst habitats in northern Guangdong Province, China. The specialized habitat preference makes it an ideal model species for study of adaptive evolution. Meanwhile, the phylogenetic relationships of F. danxiaensis in four locations under two landforms are unclear. Therefore, we sequenced its complete chloroplast (cp.) genomes and conducted comprehensive interspecific and intrageneric plastome studies.

RESULTS

The F. danxiaensis plastomes in four locations showed a typical quadripartite and circular structure that ranged from 160,832 to 161,206 bp in size, with 112 unique genes encoded. Comparative genomics showed that the plastomes of F. danxiaensis were relatively conserved with high similarity of genome organization, gene number, GC content and SSRs. While the genomes revealed higher biased codon preferences in Karst habitat than those in Danxia habitats. Eighteen and 11 divergent hotpots were identified at interspecific and intrageneric levels for species identification and further phylogenetic studies. Seven genes (clpP, accD, ccsA, ndhH, rpl20, rpoC2, and rps4) were under positive selection and may be related to adaptation. Phylogenetic analysis revealed that F. danxiaensis is sister to F. major and F. simplex. However, the interspecific relationships are not consistent with the habitat types.

CONCLUSIONS

The characteristics and interspecific relationship of F. danxiaensis plastomes provide new insights into further integration of geographical factors, environmental factors, and genetic variations on the genomic study of F. danxiaensis. Together, our study will contribute to the study of species identification, population genetics, and conservation biology of F. danxiaensis.

Genome-wide identification and functional analysis of the cellulose synthase-like gene superfamily in common oat (Avena sativa L.)

Hemicelluloses constitute approximately one-third of the plant cell wall and can be used as a dietary fiber and food additive, and as raw materials for biofuels. Although genes involved in hemicelluloses synthesis have been investigated in some model plants, no comprehensive analysis has been conducted in common oat at present. In this study, we identified and systematically analyzed the cellulose synthase-like gene (Csl) family members in common oat and investigated them using various bioinformatics tools. The results showed that there are 76 members of the oat Csl gene family distributed on 17 chromosomes, and phylogenetic analysis indicated that the 76 Csl genes belong to the CslA, CslC, CslD, CslE, CslF, CslH, and CslJ subfamilies. A total of 14 classes of cis-acting elements were identified in the promoter regions, including hormone response, light response, cell development, and defense stress elements. The collinearity analysis identified 28 pairs of segmentally duplicated genes, most of which were found on chromosomes 2D and 6A. Expression pattern analysis showed that oat Csl genes display strong tissue-specific expression; of the 76 Csl genes, 33 were significantly up-regulated in stems and 30 were up-regulated in immature seeds. The expression of most members of the AsCsl gene family is repressed by abiotic stress, while the expression of some members is up-regulated by light. Immunoelectron microscopy shows that the product of AsCsl61, a member of CslF subfamily, mediates (1,3; 1,4)-β-D-glucan synthesis in transgenic Arabidopsis. These findings provide a fundamental understanding of the structural, functional, and evolutionary features of the oat Csl genes and may contribute to our general understanding of hemicellulose biosynthesis. Moreover, this information will be helpful in designing experiments for genetic manipulation of mixed-linkage glucan (MLG) synthesis with the goal of quality improvement in oat.

Comparative analysis of codon usage patterns in the chloroplast genomes of nine forage legumes

Leguminosae is one of the three largest families of angiosperms after Compositae and Orchidaceae. It is widely distributed and grows in a variety of environments, including plains, mountains, deserts, forests, grasslands, and even waters where almost all legumes can be found. It is one of the most important sources of starch, protein and oil in the food of mankind and also an important source of high-quality forage material for animals, which has important economic significance. In our study, the codon usage patterns and variation sources of the chloroplast genome of nine important forage legumes were systematically analyzed. Meanwhile, we also constructed a phylogenetic tree based on the whole chloroplast genomes and protein coding sequences of these nine forage legumes. Our results showed that the chloroplast genomes of nine forage legumes end with A/T bases, and seven identical high-frequency (HF) codons were detected among the nine forage legumes. ENC-GC3s mapping, PR2 analysis, and neutral analysis showed that the codon bias of nine forage legumes was influenced by many factors, among which natural selection was the main influencing factor. The codon usage frequency showed that the Nicotiana tabacum and Saccharomyces cerevisiae can be considered as receptors for the exogenous expression of chloroplast genes of these nine forage legumes. The phylogenetic relationships of the chloroplast genomes and protein coding genes were highly similar, and the nine forage legumes were divided into three major clades. Among the clades Melilotus officinalis was more closely related to Medicago sativa, and Galega officinalis was more closely related to Galega orientalis. This study provides a scientific basis for the molecular markers research, species identification and phylogenetic studies of forage legumes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01421-0.

Four new complete mitochondrial genomes of Gobioninae fishes (Teleostei: Cyprinidae) and their phylogenetic implications

The subfamily Gobioninae is one of the most diverse fish groups within Cyprinidae. Their taxonomy and phylogenetic relationships are not completely resolved. In this study, the complete mitochondrial genomes (mitogenome) of four Gobioninae species (Microphysogobio elongatus, Microphysogobio chinssuensis, Gobio rivuloides and Rhinogobio nasutus) were sequenced and compared. The mitogenomes of four species ranges from 16603 bp to 16609 bp in length, consisting of 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a control region. Most PCGs had significant codon usage bias. Except for the tRNASer (GCT), all the nucleotide substitutions of tRNA loops higher than the stems could fold into a stable secondary structure. The nucleotide compositions of Gobioninae mitogenome were biased toward A/T, and NAD4 was subjected to low purification selection and had a faster evolution rate among 13 PCGs. Bayesian inference and maximum likelihood phylogenetic analyses showed the consistent results. The four sequenced species clustered together with their congener species. However, more samples and mitogenome data are needed to untangle the phylogenetic relationships among genera Microphysogobio, Romanogobio, Hugobio, Biwia and Platysmacheilus.

Thirteen complete chloroplast genomes of the costaceae family: insights into genome structure, selective pressure and phylogenetic relationships

BACKGROUND

Costaceae, commonly known as the spiral ginger family, consists of approximately 120 species distributed in the tropical regions of South America, Africa, and Southeast Asia, of which some species have important ornamental, medicinal and ecological values. Previous studies on the phylogenetic and taxonomic of Costaceae by using nuclear internal transcribed spacer (ITS) and chloroplast genome fragments data had low resolutions. Additionally, the structures, variations and molecular evolution of complete chloroplast genomes in Costaceae still remain unclear. Herein, a total of 13 complete chloroplast genomes of Costaceae including 8 newly sequenced and 5 from the NCBI GenBank database, representing all three distribution regions of this family, were comprehensively analyzed for comparative genomics and phylogenetic relationships.

RESULT

The 13 complete chloroplast genomes of Costaceae possessed typical quadripartite structures with lengths from 166,360 to 168,966 bp, comprising a large single copy (LSC, 90,802 - 92,189 bp), a small single copy (SSC, 18,363 - 20,124 bp) and a pair of inverted repeats (IRs, 27,982 - 29,203 bp). These genomes coded 111 - 113 different genes, including 79 protein-coding genes, 4 rRNA genes and 28 - 30 tRNAs genes. The gene orders, gene contents, amino acid frequencies and codon usage within Costaceae were highly conservative, but several variations in intron loss, long repeats, simple sequence repeats (SSRs) and gene expansion on the IR/SC boundaries were also found among these 13 genomes. Comparative genomics within Costaceae identified five highly divergent regions including ndhF, ycf1-D2, ccsA-ndhD, rps15-ycf1-D2 and rpl16-exon2-rpl16-exon1. Five combined DNA regions (ycf1-D2 + ndhF, ccsA-ndhD + rps15-ycf1-D2, rps15-ycf1-D2 + rpl16-exon2-rpl16-exon1, ccsA-ndhD + rpl16-exon2-rpl16-exon1, and ccsA-ndhD + rps15-ycf1-D2 + rpl16-exon2-rpl16-exon1) could be used as potential markers for future phylogenetic analyses and species identification in Costaceae. Positive selection was found in eight protein-coding genes, including cemA, clpP, ndhA, ndhF, petB, psbD, rps12 and ycf1. Maximum likelihood and Bayesian phylogenetic trees using chloroplast genome sequences consistently revealed identical tree topologies with high supports between species of Costaceae. Three clades were divided within Costaceae, including the Asian clade, Costus clade and South American clade. Tapeinochilos was a sister of Hellenia, and Parahellenia was a sister to the cluster of Tapeinochilos + Hellenia with strong support in the Asian clade. The results of molecular dating showed that the crown age of Costaceae was about 30.5 Mya (95% HPD: 14.9 - 49.3 Mya), and then started to diverge into the Costus clade and Asian clade around 23.8 Mya (95% HPD: 10.1 - 41.5 Mya). The Asian clade diverged into Hellenia and Parahellenia at approximately 10.7 Mya (95% HPD: 3.5 - 25.1 Mya).

CONCLUSION

The complete chloroplast genomes can resolve the phylogenetic relationships of Costaceae and provide new insights into genome structures, variations and evolution. The identified DNA divergent regions would be useful for species identification and phylogenetic inference in Costaceae.

Whole mitochondrial and chloroplast genome sequencing of Tunisian date palm cultivars: diversity and evolutionary relationships

BACKGROUND

Date palm (Phoenix dactylifera L.) is the most widespread crop in arid and semi-arid regions and has great traditional and socioeconomic importance, with its fruit well-known for its high nutritional and health value. However, the genetic variation of date palm cultivars is often neglected. The advent of high-throughput sequencing has made possible the resequencing of whole organelle (mitochondria and chloroplast) genomes to explore the genetic diversity and phylogenetic relationships of cultivated plants with unprecedented detail.

RESULTS

Whole organelle genomes of 171 Tunisian accessions (135 females and 36 males) were sequenced. Targeted bioinformatics pipelines were used to identify date palm haplotypes and genome variants, aiming to provide variant annotation and investigate patterns of evolutionary relationship. Our results revealed the existence of unique haplotypes, identified by 45 chloroplastic and 156 mitochondrial SNPs. Estimation of the effect of these SNPs on genes functions was predicted in silico.

CONCLUSIONS

The results of this study have important implications, in the light of ongoing environmental changes, for the conservation and sustainable use of the genetic resources of date palm cultivars in Tunisia, where monoculture threatens biodiversity leading to genetic erosion. These data will be useful for breeding and genetic improvement programs of the date palm through selective cross-breeding.

The First Genome-Wide Survey of Shortbelly Eel (Dysomma anguillare Barnard, 1923) to Provide Genomic Characteristics, Microsatellite Markers and Complete Mitogenome Information

Dysomma anguillare is a demersal eel widespread distributing in tropical waters of the Indo-West Pacific and Atlantic. As an important component of the coastal fishery and marine ecosystem, the lack of genomic information for this species severely restricts the progress of relevant researches. In this study, the abecedarian genome-wide characteristics and phylogenetic relationships analyses were carried out based on next-generation sequencing (NGS) platform. The revised genome size was approximately 1 310 Mb, with the largest scaffold length reaching 23 878 bp through K-mer (K = 17) method. The heterozygosity, repetitive rate and average GC content were about 0.94%, 51.93% and 42.23%, respectively. A total of 1 160 104 microsatellite motifs were identified from the de novo assembled genome of D. anguillare, in which dinucleotide repeats accounted for the largest proportion (592 234, 51.05%), the highest occurrence frequency (14.58%) as well as the largest relative abundance (379.27/Mb). The high-polymorphic and moderate-polymorphic loci composed around 73% of the total single sequence repeats (SSRs), showing a latent capacity for subsequent population genetic structure and genetic diversity appraisal researches. Another byproduct of whole-genome sequencing, the double-stranded and circular mitogenome (16 690 bp) was assembled to investigate the evolutionary relationships of D. anguillare. The phylogenic tree constructed with maximum likelihood (ML) method showed that D. anguillare was closely related to Synaphobranchidae species, and the molecular systematic results further supported classical taxonomy status of D. anguillare.

A new species of Notocotylus (Digenea: Notocotylidae) from the black-necked swan Cygnus melancorhyphus (Molina) of Argentina

Notocotylus cygni n. sp. is described here, taken from the intestine of the black-necked swan Cygnus melancoryphus (Molina) of Patagonia, Argentina. This new species differs from other members of the genus Notocotylus by having the genital pore anterior to the caecal bifurcation (located slightly posterior to oral sucker) and the unequal number and arrangement of ventral papillae (2-3 in the lateral rows and 10-12 in the median row). Phylogenetic analyses of the 28S and ITS1-5.8S ribosomal DNA (rRNA) sequences of the new species and other notocotylid trematodes available in GenBank indicate that N. cygni n. sp. is a sister taxon of Notocotylus fosteri Kinsella et Tkach, 2005, a trematode of the intestine of the rice rat Oryzomys palustris of Florida, United States. The new species differs from N. fosteri in the unequal number and arrangement of ventral papillae, number of uterine loops, size of the egg, definitive hosts (birds vs. mammals), and disparate environment and geographical distribution (freshwater environment in Patagonia vs. salt marsh in North America). This is the eighth species of Notocotylus reported from birds in Argentina, and the ninth species from the family Notocotylidae recorded in black-necked swans in South America.

Plastome evolution in the genus Sium (Apiaceae, Oenantheae) inferred from phylogenomic and comparative analyses

BACKGROUND

Sium L. (Apiaceae) is a small genus distributed primarily in Eurasia, with one species also occurring in North America. Recently, its circumscription has been revised to include 10 species, however, the phylogenetic relationships within its two inclusive clades were poorly supported or collapsed in previous studies based on nuclear ribosomal DNA ITS or cpDNA sequences. To identify molecular markers suitable for future intraspecific phylogeographic and population genetic studies, and to evaluate the efficacy of plastome in resolving the phylogenetic relationships of the genus, the complete chloroplast (cp) genomes of six Sium species were sequenced.

RESULTS

The Sium plastomes exhibited typical quadripartite structures of Apiaceae and most other higher plant plastid DNAs, and were relatively conserved in their size (153,029-155,006 bp), gene arrangement and content (with 114 unique genes). A total of 61-67 SSRs, along with 12 highly divergent regions (trnQ, trnG-atpA, trnE-trnT, rps4-trnT, accD-psbI, rpl16, ycf1-ndhF, ndhF-rpl32, rpl32-trnL, ndhE-ndhG, ycf1a and ycf1b) were discovered in the plastomes. No significant IR length variation was detected showing that plastome evolution was conserved within this genus. Phylogenomic analysis based on whole chloroplast genome sequences produced a highly resolved phylogenetic tree, in which the monophyly of Sium, as well as the sister relationship of its two inclusive clades were strongly supported.

CONCLUSIONS

The plastome sequences could greatly improve phylogenetic resolution, and will provide genomic resources and potential markers useful for future studies of the genus.

Insights into phylogenetic relationships in Pinus inferred from a comparative analysis of complete chloroplast genomes

BACKGROUND

Pinus is the largest genus of Pinaceae and the most primitive group of modern genera. Pines have become the focus of many molecular evolution studies because of their wide use and ecological significance. However, due to the lack of complete chloroplast genome data, the evolutionary relationship and classification of pines are still controversial. With the development of new generation sequencing technology, sequence data of pines are becoming abundant. Here, we systematically analyzed and summarized the chloroplast genomes of 33 published pine species.

RESULTS

Generally, pines chloroplast genome structure showed strong conservation and high similarity. The chloroplast genome length ranged from 114,082 to 121,530 bp with similar positions and arrangements of all genes, while the GC content ranged from 38.45 to 39.00%. Reverse repeats showed a shrinking evolutionary trend, with IRa/IRb length ranging from 267 to 495 bp. A total of 3,205 microsatellite sequences and 5,436 repeats were detected in the studied species chloroplasts. Additionally, two hypervariable regions were assessed, providing potential molecular markers for future phylogenetic studies and population genetics. Through the phylogenetic analysis of complete chloroplast genomes, we offered novel opinions on the genus traditional evolutionary theory and classification.

CONCLUSION

We compared and analyzed the chloroplast genomes of 33 pine species, verified the traditional evolutionary theory and classification, and reclassified some controversial species classification. This study is helpful in analyzing the evolution, genetic structure, and the development of chloroplast DNA markers in Pinus.

Comprehensive phylogenomic analyses reveal that order Armophorida is most closely related to class Armophorea (Protista, Ciliophora)

Ciliate species within the class Armophorea are widely distributed in various anaerobic environments, hence they are of great interest to researchers studying evolution and adaptation of eukaryotes to extreme habitats. However, phylogenetic relationships within the class remain largely elusive, most especially assignment of the order Armophorida and classification within the family Metopidae. In this study, we newly sequenced transcriptomes and the SSU rDNA of five armophorean species, Sulfonecta cf. uniserialis (order Armophorida), Nyctotheroides sp. (order Clevelandellida), and Metopus major, M. paraes, and Brachonella contorta (order Metopida). Comprehensive phylogenomic analyses revealed that Armophorea was most closely related to Muranotrichea and Parablepharismea. Our results indicate that the order Armophorida either belongs to Armophorea or represents a new class within APM (Armophorea-Parablepharismea-Muranotrichea). Analyses combining ecological niches and molecular trees showed that APM species might descend from an anaerobic free-living ciliate species. Existing molecular phylogenomic/phylogenetic and morphological evidence indicate that the family Metopidae is non-monophyletic and should be further classified with inclusion of additional lines of evidences. Our results provide new insights into the long-debated relationships within Armophorea.

Phylogenetic study and comparison of different TbpB obtained from Glaesserella parasuis present in Spanish clinical isolates

Glaesserella parasuis (Gp) is the etiological agent of Glässer's disease (GD), which causes important economic losses for the pig intensive production worldwide. This organism uses a smart protein-based receptor to acquire specifically iron from the porcine transferrin. This surface receptor consists of transferrin-binding protein A (TbpA) and transferrin-binding protein B (TbpB). TbpB has been considered the most promising antigen to formulate a based-protein vaccine with broad-spectrum of protection against GD. The purpose of our study was to determine the capsular diversity of Gp clinical isolates collected in different Spanish regions between 2018 and 2021. A total of 68 Gp isolates were recovered from porcine respiratory or systemic samples. A species-specific PCR based on tbpA gene, followed by multiplex PCR for typing Gp isolates were performed. Serovars 5, 10, 2, 4 and 1 were the most prevalent and involved almost 84% of isolates. TbpB amino acid sequences from 59 of these isolates were analyzed, and a total of ten clades could be established. All of them showed a wide diversity with respect to capsular type, anatomical isolation site and geographical origin, with minor exceptions. Regardless of the serovars, the in silico analysis of TbpB sequences revealed that a vaccine based on a TbpB recombinant protein could potentially prevent Glässer's disease outbreaks in Spain.

Insights into the phylogeny and chloroplast genome evolution of Eriocaulon (Eriocaulaceae)

BACKGROUND

Eriocaulon is a wetland plant genus with important ecological value, and one of the famous taxonomically challenging groups among angiosperms, mainly due to the high intraspecific diversity and low interspecific variation in the morphological characters of species within this genus. In this study, 22 samples representing 15 Eriocaulon species from China, were sequenced and combined with published samples of Eriocaulon to test the phylogenetic resolution using the complete chloroplast genome. Furthermore, comparative analyses of the chloroplast genomes were performed to investigate the chloroplast genome evolution of Eriocaulon.

RESULTS

The 22 Eriocaulon chloroplast genomes and the nine published samples were proved highly similar in genome size, gene content, and order. The Eriocaulon chloroplast genomes exhibited typical quadripartite structures with lengths from 150,222 bp to 151,584 bp. Comparative analyses revealed that four mutation hotspot regions (psbK-trnS, trnE-trnT, ndhF-rpl32, and ycf1) could serve as effective molecular markers for further phylogenetic analyses and species identification of Eriocaulon species. Phylogenetic results supported Eriocaulon as a monophyletic group. The identified relationships supported the taxonomic treatment of section Heterochiton and Leucantherae, and the section Heterochiton was the first divergent group. Phylogenetic tree supported Eriocaulon was divided into five clades. The divergence times indicated that all the sections diverged in the later Miocene and most of the extant Eriocaulon species diverged in the Quaternary. The phylogeny and divergence times supported rapid radiation occurred in the evolution history of Eriocaulon.

CONCLUSION

Our study mostly supported the taxonomic treatment at the section level for Eriocaulon species in China and demonstrated the power of phylogenetic resolution using whole chloroplast genome sequences. Comparative analyses of the Eriocaulon chloroplast genome developed molecular markers that can help us better identify and understand the evolutionary history of Eriocaulon species in the future.

Start codon targeted (SCoT) polymorphism marker in plant genome analysis: current status and prospects

The present review illustrates a comprehensive overview of the start codon targeted (SCoT) polymorphism marker and their utilization in various applications related to genetic and genomic studies. Start codon targeted (SCoT) polymorphism marker, a targeted fingerprinting marker technique, has gained considerable importance in plant genetics, genomics, and molecular breeding due to its many desirable features. SCoT marker targets the region flanking the start codon, a highly conserved region in plant genes. Therefore, it can distinguish genetic variations in a specific gene that link to a specific trait. It is a simple, novel, cost-effective, highly polymorphic, and reproducible molecular marker for which there is no need for prior sequence information. In the recent past, SCoT markers have been employed in many commercially important and underutilized plant species for a variety of applications, including genetic diversity analysis, interspecific/generic genetic relationships, cultivar/hybrid/species identification, sex determination, construction of linkage map, association mapping/analysis, differential gene expression, and genetic fidelity analysis of tissue culture-raised plants. The main aim of this review is to provide up-to-date information on SCoT markers and their application in many commercially important and underutilized plant species, mainly progress made in the last 8-10 years.

The plastid genome of twenty-two species from Ferula, Talassia, and Soranthus: comparative analysis, phylogenetic implications, and adaptive evolution

BACKGROUND

The Ferula genus encompasses 180-185 species and is one of the largest genera in Apiaceae, with many of Ferula species possessing important medical value. The previous studies provided more information for Ferula, but its infrageneric relationships are still confusing. In addition, its genetic basis of its adaptive evolution remains poorly understood. Plastid genomes with more variable sites have the potential to reconstruct robust phylogeny in plants and investigate the adaptive evolution of plants. Although chloroplast genomes have been reported within the Ferula genus, few studies have been conducted using chloroplast genomes, especially for endemic species in China.

RESULTS

Comprehensively comparative analyses of 22 newly sequenced and assembled plastomes indicated that these plastomes had highly conserved genome structure, gene number, codon usage, and repeats type and distribution, but varied in plastomes size, GC content, and the SC/IR boundaries. Thirteen mutation hotspot regions were detected and they would serve as the promising DNA barcodes candidates for species identification in Ferula and related genera. Phylogenomic analyses with high supports and resolutions showed that Talassia transiliensis and Soranthus meyeri were nested in the Ferula genus, and thus they should be transferred into the Ferula genus. Our phylogenies also indicated the monophyly of subgenera Sinoferula and subgenera Narthex in Ferula genus. Twelve genes with significant posterior probabilities for codon sites were identified in the positively selective analysis, and their function may relate to the photosystem II, ATP subunit, and NADH dehydrogenase. Most of them might play an important role to help Ferula species adapt to high-temperatures, strong-light, and drought habitats.

CONCLUSION

Plastome data is powerful and efficient to improve the support and resolution of the complicated Ferula phylogeny. Twelve genes with significant posterior probabilities for codon sites were helpful for Ferula to adapt to the harsh environment. Overall, our study supplies a new perspective for comprehending the phylogeny and evolution of Ferula.

Identification, characterization and prevalence in southern China of a new iflavirus in the leafhopper Recilia dorsalis (Hemiptera: Cicadellidae)

The leafhopper Recilia dorsalis (Hemiptera: Cicadellidae) is not only a significant pest in agriculture but also an important vector involved in transmitting numerous pathogens that are known to cause economic losses by affecting rice crops. Here, a new iflavirus was discovered in the leafhopper R. dorsalis by employing a transcriptomic approach. The complete viral genome was determined to be 10,711 nucleotides (nt) in length and contains a single open reading frame (ORF) encoding a putative polyprotein comprised of 3,161 amino acids (aa), which is flanked by 5' and 3' untranslated regions. The full viral genome nt and the deduced polyprotein aa sequence showed the highest similarity (71.6% and 77.8%, respectively) with Langfang leafhopper iflavirus. Phylogenetic analysis based on the RdRp domain indicated that the isolated virus, which we have tentatively named Recilia dorsalis iflavirus 2 (RdIV2), is clustered with the members of the family Iflaviridae. Moreover, the results of our surveys indicate that RdIV2 predominates in southwestern Guangdong and southeastern Guangxi, China, and was absent in the other three species of leafhoppers; Nephotettix cincticeps, N. virescens and N. nigropictus. Notably, R. dorsalis was found to be co-infected with RdIV2 and rice stripe mosaic virus (RSMV; a well-known rice-infecting virus vectored by R. dorsalis) in rice fields, although the co-infection rate is low.

Two new syntopic species of glassfrogs (Amphibia, Centrolenidae, Centrolene) from the southwestern Andes of Ecuador

We describe two new species of glassfrogs of the genus Centrolene living in syntopy at La Enramada, province of Azuay, southwestern Ecuador. They were found in a small creek in montane evergreen forests at 2,900 m elevation. The first new species is distinguished from all other members of the genus Centrolene by having the following combination of characters: dentigerous process of vomer absent; sloping snout in lateral view; thick, white labial stripe and a faint white line between the lip and anterior ¼ of body; humeral spine in adult males; parietal peritoneum covered by iridophores, visceral peritonea translucent (except pericardium); ulnar and tarsal ornamentation; dorsal skin shagreen with dispersed warts; uniform green dorsum with light yellowish green warts; and green bones. The new species is remarkable by being sister to a species from the opposite Andean versant, C. condor. The second new species is distinguished from all other Centrolene by having the following combination of characters: dentigerous process of vomer absent; round snout in lateral view; thin, yellowish labial stripe with a row of white tubercles between the lip and arm insertion, and a yellowish line between arm insertion and groin; uniform green dorsum; humeral spine in adult males; parietal peritoneum covered by iridophores, visceral peritonea translucent (except pericardium); dorsal skin shagreen with dispersed spicules; ulnar and tarsal ornamentation; and green bones. The second new species is sister to C. sabini and an undescribed species of Centrolene from southeastern Ecuador. Based on nuclear and mitochondrial DNA sequences, we present a new phylogeny for Centrolene and comment on the phylogenetic relationships inside the genus.

Evaluation of genome size and phylogenetic relationships of the Saccharum complex species

"Saccharum complex" is a hypothetical group of species, which is supposed to be involved in the origin of modern sugarcane, and displays large genomes and complex chromosomal alterations. The utilization of restricted parents in breeding programs of modern cultivated sugarcane has resulted in a genetic blockage, which controlled its improvement because of the limited genetic diversity. The use of wild relatives is an effective way to broaden the genetic composition of cultivated sugarcane. Due to the infrequent characterization of genomes, the potential of wild relatives is diffused in improving the cultivated sugarcane. To characterize the genomes of the wild relatives, the genome size and phylogenetic relationships among eight species, including Saccharum spontaneum, Erianthus arundinaceus, E. fulvus, E. rockii, Narenga porphyrocoma, Miscanthus floridulus, Eulalia quadrinervis, and M. sinensis were evaluated based on flow cytometry, genome surveys, K-mer analysis, chloroplast genome sequencing, and whole-genome SNPs analysis. We observed highly heterozygous genomes of S. spontaneum, E. rockii, and E. arundinaceus and the highly repetitive genome of E. fulvus. The genomes of Eulalia quadrinervis, N. porphyrocoma, M. sinensis, and M. floridulus were highly complex. Phylogenetic results of the two approaches were dissimilar, however, both indicate E. fulvus displayed closer relationships to Miscanthus and Saccharum than other species of Saccharum complex. Eulalia quadrinervis was more closely related to M. floridulus than M. sinensis; E. arundinaceus differ significantly from Miscanthus, Narenga, and Saccharum, but was relatively close to Erianthus. We proved the point of E. rockii and E. fulvus should not be classified as one genus, and E. fulvus should be classified as the Saccharum genus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03338-5.

Analysis of complete chloroplast genome sequences and insight into the phylogenetic relationships of Ferula L

Background

Ferula L. is one of the largest and most taxonomically complicated genera as well as being an important medicinal plant resource in the family Apiaceae. To investigate the plastome features and phylogenetic relationships of Ferula and its neighboring genera Soranthus Ledeb., Schumannia Kuntze., and Talassia Korovin, we sequenced 14 complete plastomes of 12 species. 

Results

The size of the 14 complete chloroplast genomes ranged from 165,607 to 167,013 base pairs (bp) encoding 132 distinct genes (87 protein-coding, 37 tRNA, and 8 rRNA genes), and showed a typical quadripartite structure with a pair of inverted repeats (IR) regions. Based on comparative analysis, we found that the 14 plastomes were similar in codon usage, repeat sequence, simple sequence repeats (SSRs), and IR borders, and had significant collinearity. Based on our phylogenetic analyses, Soranthus, Schumannia, and Talassia should be considered synonymous with Ferula. Six highly divergent regions (rps16/trnQ-UUG, trnS-UGA/psbZ, psbH/petB, ycf1/ndhF, rpl32, and ycf1) were also detected, which may represent potential molecular markers, and combined with selective pressure analysis, the weak positive selection gene ccsA may be a discriminating DNA barcode for Ferula species.

Conclusion

Plastids contain abundant informative sites for resolving phylogenetic relationships. Combined with previous studies, we suggest that there is still much room for improvement in the classification of Ferula. Overall, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of this genus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08868-z.

The complete chloroplast genomes of Tetrastigma hemsleyanum (Vitaceae) from different regions of China: molecular structure, comparative analysis and development of DNA barcodes for its geographical origin discrimination

BACKGROUND

Tetrastigma hemsleyanum is a valuable traditional Chinese medicinal plant widely distributed in the subtropical areas of China. It belongs to the Cayratieae tribe, family Vitaceae, and exhibited significant anti-tumor and anti-inflammatory activities. However, obvious differences were observed on the quality of T. hemsleyanum root from different regions, requiring the discrimination strategy for the geographical origins.

RESULT

This study characterized five complete chloroplast (cp) genomes of T. hemsleynum samples from different regions, and conducted a comparative analysis with other representing species from family Vitaceae to reveal the structural variations, informative markers and phylogenetic relationships. The sequenced cp genomes of T. hemsleyanum exhibited a conserved quadripartite structure with full length ranging from 160,124 bp of Jiangxi Province to 160,618 bp of Zhejiang Province. We identified 112 unique genes (80 protein-coding, 28 tRNA and 4 rRNA genes) in the cp genomes of T. hemsleyanum with highly similar gene order, content and structure. The IR contraction/expansion events occurred on the junctions of ycf1, rps19 and rpl2 genes with different degrees, causing the differences of genome sizes in T. hemsleyanum and Vitaceae plants. The number of SSR markers discovered in T. hemsleyanum was 56-57, exhibiting multiple differences among the five geographic groups. Phylogenetic analysis based on conserved cp genome proteins strongly grouped the five T. hemsleyanum species into one clade, showing a sister relationship with T. planicaule. Comparative analysis of the cp genomes from T. hemsleyanum and Vitaceae revealed five highly variable spacers, including 4 intergenic regions and one protein-coding gene (ycf1). Furthermore, five mutational hotspots were observed among T. hemsleyanum cp genomes from different regions, providing data for designing DNA barcodes trnL and trnN. The combination of molecular markers of trnL and trnN clustered the T. hemsleyanum samples from different regions into four groups, thus successfully separating specimens of Sichuan and Zhejiang from other areas.

CONCLUSION

Our study obtained the chloroplast genomes of T. hemsleyanum from different regions, and provided a potential molecular tracing tool for determining the geographical origins of T. hemsleyanum, as well as important insights into the molecular identification approach and and phylogeny in Tetrastigma genus and Vitaceae family.

Plastome structure, phylogenomics and evolution of plastid genes in Swertia (Gentianaceae) in the Qing-Tibetan Plateau

BACKGROUND

The genus Swertia is of great medicinal importance and one of the most taxonomically challenging taxa within Gentianaceae, largely due to the morphological similarities of species within this genus and with its closely related genera. Previous molecular studies confirmed its polyphyly but suffered from low phylogenetic resolutions because only limited sequence loci were used. Thus, we conducted the structural, gene evolutionary, and phylogenetic analyses of 11 newly obtained plastomes of Swertia. Our result greatly improved the phylogenetic resolutions in Swertia, shed new light on the plastome evolution and phylogenetic relationships of this genus.

RESULTS

The 11 Swertia plastomes together with the published seven species proved highly similar in overall size, structure, gene order, and content, but revealed some structural variations caused by the expansion and contraction of the IRb region into the LSC region, due to the heterogeneous length of the ψycf1. The gene rps16 was found to be in a state flux with pseudogenes or completely lost. Similar situation was also documented in other genera of Gentianaceae. This might imply loss of the gene in the common ancestor of Gentianaceae. The distribution plot of ENC vs. GC3 showed all these plastomes arranging very close in the Wright line with an expected ENC value (49-52%), suggesting the codon usage of Swertia was mainly constrained by a GC mutation bias. Most of the genes remained under the purifying selection, however, the cemA was identified under positive selection, possibly reflecting an adaptive response to low CO₂ atmospheric conditions during the Late Miocene. Our phylogenomic analyses, based on 74 protein-coding genes (CDS), supported the polyphyly of Swertia with its close allies in the subtribe Swertiinae, presumably due to recent rapid radiation. The topology inferred from our phylogenetic analyses partly supported the current taxonomic treatment. Finally, several highly variable loci were identified, which can be used in future phylogenetic studies and accurate identification of medicinal genuineness of Swertia.

CONCLUSIONS

Our study confirmed the polyphyly of Swertia and demonstrated the power of plastome phylogenomics in improvement of phylogenetic resolution, thus contributing to a better understanding of the evolutionary history of Swertia.

The species of Oxytropis DC. of section Gloeocephala Bunge (Fabaceae) from Northeast Asia: genetic diversity and relationships based on sequencing of the intergenic spacers of cpDNA and ITS nrDNA

Phylogenetic relationships within Oxytropis DC. sect. Gloeocephala Bunge from Northeast Asia were studied using plastid intergenic spacers (psbA-trnH + trnL-trnF + trnS-trnG) and ITS nrDNA. Populations of O. anadyrensis Vass., O. borealis DC., O. middendorffii Trautv., O. trautvetteri Meinsh., and O. vasskovskyi Jurtz. were monomorphic or characterised by a low level of chloroplast genetic diversity (h varied from 0.143 to 0.692, and π from 0.0001 to 0.0005). Presumably, the low genetic diversity was a result of the severe bottlenecks during Pleistocene glaciation-interglacial cycles. Twenty chlorotypes were identified; species studied had no shared chlorotypes. Chlorotypes of O. anadyrensis, O. borealis, and O. middendorffii formed two lineages each, while the chlorotypes of O. trautvetteri and O. vasskovskyi formed one separate lineage each in the phylogenetic network. There were specific diagnostic markers of cpDNA in each lineage, excluding O. vasskovskyi. The presence of a species-specific diagnostic marker in O. trautvetteri and specific markers in two lineages of O. anadyrensis support circumscribing these taxa as independent species. Regarding ITS nrDNA polymorphism, five ribotypes were detected. The differences revealed in plastid and nuclear genomes of Oxytropis sect. Gloeocephala confirmed that the Asian sector of Megaberingia was the main centre of diversification of arctic legumes.

Genome-wide bioinformatics analysis of Cellulose Synthase gene family in common bean (Phaseolus vulgaris L.) and the expression in the pod development

Background

CesA and Csl gene families, which belong to the cellulose synthase gene superfamily, plays an important role in the biosynthesis of the plant cell wall. Although researchers have investigated this gene superfamily in several model plants, to date, no comprehensive analysis has been conducted in the common bean.

Results

In this study, we identified 39 putative cellulose synthase genes from the common bean genome sequence. Then, we performed a bioinformatics analysis of this gene family involving sequence alignment, phylogenetic analysis, gene structure, collinearity analysis and chromosome location. We found all members possess a cellulose_synt domain. Phylogenetic analysis revealed that these cellulose synthase genes may be classified into five subfamilies, and that members in the same subfamily share conserved exon-intron distribution and motif compositions. Abundant and distinct cis-acting elements in the 2 k basepairs upstream regulatory regions indicate that the cellulose synthase gene family may plays a vital role in the growth and development of common bean. Moreover, the 39 cellulose synthase genes are distributed on 10 of the 11 chromosomes. Additionally expression analysis shows that all CesA/Csl genes selected are constitutively expressed in the pod development.

Conclusions

This research reveals both the putative biochemical and physiological functions of cellulose synthase genes in common bean and implies the importance of studying non-model plants to understand the breadth and diversity of cellulose synthase genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01026-0.

Reinvestigating the phylogeny of Myriapoda with more extensive taxon sampling and novel genetic perspective

Background

There have been extensive debates on the interrelationships among the four major classes of Myriapoda-Chilopoda, Symphyla, Diplopoda, and Pauropoda. The core controversy is the position of Pauropoda; that is, whether it should be grouped with Symphyla or Diplopoda as a sister group. Two recent phylogenomic studies separately investigated transcriptomic data from 14 and 29 Myriapoda species covering all four groups along with outgroups, and proposed two different topologies of phylogenetic relationships.

Methods

Building on these studies, we extended the taxon sampling by investigating 39 myriapods and integrating the previously available data with three new transcriptomic datasets generated in this study. Our analyses present the phylogenetic relationships among the four major classes of Myriapoda with a more abundant taxon sampling and provide a new perspective to investigate the above-mentioned question, where visual genes' identification were conducted. We compared the appearance pattern of genes, grouping them according to their classes and the visual pathways involved. Positive selection was detected for all identified visual genes between every pair of 39 myriapods, and 14 genes showed positive selection among 27 pairs.

Results

From the results of phylogenomic analyses, we propose that Symphyla is a sister group of Pauropoda. This stance has also received strong support from tree inference and topology tests.

Genome-wide and expression analysis of B-box gene family in pepper

BACKGROUND

BBX transcription factors are a kind of zinc finger transcription factors with one or two B-box domains, which partilant in plant growth, development and response to abiotic or biotic stress. The BBX family has been identified in Arabidopsis, rice, tomato and some other model plant genomes.

RESULTS

Here, 24 CaBBX genes were identified in pepper (Capsicum annuum L.), and the phylogenic analysis, structures, chromosomal location, gene expression patterns and subcellular localizations were also carried out to understand the evolution and function of CaBBX genes. All these CaBBXs were divided into five classes, and 20 of them distributed in 11 of 12 pepper chromosomes unevenly. Most duplication events occurred in subgroup I. Quantitative RT-PCR indicated that several CaBBX genes were induced by abiotic stress and hormones, some had tissue-specific expression profiles or differentially expressed at developmental stages. Most of CaBBX members were predicated to be nucleus-localized in consistent with the transient expression assay by onion inner epidermis of the three tested CaBBX members (CaBBX5, 6 and 20).

CONCLUSION

Several CaBBX genes were induced by abiotic stress and exogenous phytohormones, some expressed tissue-specific and variously at different developmental stage. The detected CaBBXs act as nucleus-localized transcription factors. Our data might be a foundation in the identification of CaBBX genes, and a further understanding of their biological function in future studies.

Diversity and phylogenetic relationships of haemosporidian and hemogregarine parasites in Australian lizards

Apicomplexa is a large monophyletic phylum of unicellular, parasitic organisms. Reptiles are hosts to both haemosporidian (Haemosporida) and hemogregarine (Eucoccidiorida) apicomplexan blood parasites. Within reptiles our understanding of their diversity remains limited, with a paucity of information from Australia, despite a high diversity of squamates (snakes and lizards). We provide a preliminary assessment of haemosporidian and hemogregarine diversity occurring in lizards across northern tropical Australia, building on existing data with results from a microscopy and genetic assessment. We screened total of 233 blood slides using microscopy and detected hemogregarines in 25 geckos, 2 skinks and 1 agamid, while haemosporidians were detected in 13 geckos. DNA sequencing of 28 samples of the hemogregarine 18S rRNA (∼900 bp) nuclear gene revealed five lineages of Australian lizard hemogregarines within heteroxenous adeleids. We sequenced 10 samples of Haemosporida mtDNA (cytb & coI: ∼1313 bp) and phylogenetic analysis with 30 previously published sequences revealed that the Australian Haemosporida grouped within the Haemoproteidae but were not supported as a monophyletic clade. Our results demonstrate that there is significant undocumented evolutionary diversity in Australian lizard haemosporidian and hemogregarine parasites, with preliminary evidence of significantly higher infection rates in geckos.

Evolutionary characteristics and phylogeny of cotton chloroplast tRNAs

The novel structural variations were identified in cotton chloroplast tRNAs and gene loss events were more obvious than duplications in chloroplast tRNAs. Transfer RNAs (tRNA) have long been believed an evolutionary-conserved molecular family, which play the key roles in the process of protein biosynthesis in plant life activities. In this study, we detected the evolutionary characteristics and phylogeny of chloroplast tRNAs in cotton plants, an economic and fibered important taxon in the world. We firstly annotated the chloroplast tRNAs of 27 Gossypium species to analyze their genetic composition, structural characteristics and evolution. Compared with the traditional view of evolutionary conservation of tRNA, some novel tRNA structural variations were identified in cotton plants. I.g., tRNA^Val-UAC and tRNA^Ile-GAU only contained one intron in the anti-condon loop region of tRNA secondary structure, respectively. In the variable region, some tRNAs contained a circle structure with a few nucleotides. Interestingly, the calculation result of free energy indicated that the variation of novel tRNAs contributed to the stability of tRNA structure. Phylogenetic analysis suggested that chloroplast tRNAs have evolved from multiple common ancestors, and the tRNA^Met seemed to be an ancestral tRNA, which can be duplicated and diversified to produce other tRNAs. The chloroplast tRNAs contained a group I intron in cotton plants, and the evolutionary analysis of introns indicated that group I intron of chloroplast tRNA originated from cyanobacteria. Analysis of gene duplication and loss events showed that gene loss events were more obvious than duplications in Gossypium chloroplast tRNAs. Additionally, we found that the rate of transition was higher than the ones of transversion in cotton chloroplast tRNAs. This study provided new insights into the structural characteristics and evolution of chloroplast tRNAs in cotton plants.

Complete chloroplast genome sequences of five Bruguiera species (Rhizophoraceae): comparative analysis and phylogenetic relationships

Bruguiera is a genus of true mangroves that are mostly distributed in the Indo-West Pacific region. However, the number of published whole chloroplast genome sequences of Bruguiera species are limited. Here, the complete chloroplast sequences of five Bruguiera species were sequenced and assembled using Illumina data. The chloroplast genomes of B. gymnorhiza, B. hainesii, B. cylindrica, B. parviflora and B. sexangula were assembled into 161,195, 164,295, 164,297, 163,228 and 164,170 bp, respectively. All chloroplast genomes contain 37 tRNA and eight rRNA genes, with either 84 or 85 protein-coding genes. A comparative analysis of these genomes revealed high similarity in gene structure, gene order and boundary position of the LSC, SSC and two IR regions. Interestingly, B. gymnorhiza lost a rpl32 gene in the SSC region. In addition, a ndhF gene in B. parviflora straddles both the SSC and IRB boundary regions. These genes reveal differences in chloroplast evolution among Bruguiera species. Repeats and SSRs in the chloroplast genome sequences were found to be highly conserved between B. cylindrica and B. hainesii as well as B. gymnorhiza and B. sexangula indicating close genetic relationships based on maternal inheritance. Notably, B. hainesii, which is considered a hybrid between B. gymnorhiza and B. cylindrica, appears to have inherited the chloroplast from B. cylindrica. Investigating the effects of selection events on shared protein-coding genes showed a positive selection in rps7 and rpl36 genes in all species compared to land-plant species. A phylogenetic analysis, based on 59 conserved chloroplast protein-coding genes, showed strong support that all Bruguiera species are in the clade Rhizophoraceae. This study provides valuable genetic information for the study of evolutionary relationships and population genetics in Bruguiera and other mangrove species.

In silico Structural, Functional and Phylogenetic Analyses of cellulase from Ruminococcus albus

Background

Cellulose is the primary component of the plant cell wall and an important source of energy for the ruminant and microbial protein synthesis in the rumen. Cell wall content is digested by anaerobic fermentation activity mainly of bacteria belonging to species Fibrobacter succinogenes, Ruminicoccus albus, Ruminococcus flavefaciens, and Butyrivibrio fibrisolvens. Bacteria belonging to the species Ruminococcus albus contain cellulosomes that enable it to adhere to and digest cellulose, and its genome encodes cellulases and hemicellulases.

This study aimed to perform an in silico comparative characterization and functional analysis of cellulase from Ruminococcus albus to explore physicochemical properties and to estimate primary, secondary, and tertiary structure using various bio-computational tools.

The protein sequences of cellulases belonging to 6 different Ruminococcus albus strains were retrieved using UniProt. In in silico composition of amino acids, basic physicochemical characteristics were analyzed using ProtParam and Protscale. Multiple sequence alignment of retrieved sequences was performed using Clustal Omega and the phylogenetic tree was constructed using Mega X software. Bioinformatics tools are used to better understand and determine the 3D structure of cellulase. The predicted model was refined by ModRefiner. Structure alignment between the best-predicted model and the template is applied to evaluate the similarity between structures.

Results

In this study are demonstrated several physicochemical characteristics of the cellulase enzyme. The instability index values indicate that the proteins are highly stable. Proteins are dominated by random coils and alpha helixes. The aliphatic index was higher than 71 providing information that the proteins are highly thermostable. No transmembrane domain was found in the protein, and the enzyme is extracellular and moderately acidic. The best tertiary structure model of the enzyme was obtained by the use of Raptor X, which was refined by ModRefiner. Raptor X suggested the 6Q1I_A as one of the best homologous templates for the predicted 3D protein structure. Ramachandran plot analysis showed that 90.1% of amino acid residues are within the most favored regions.

Conclusions

This study provides for the first time insights about the physicochemical properties, structure, and function of cellulase, from Ruminococcus albus, that will help for detection and identification of such enzyme in vivo or in silico.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-021-00162-x.

Inter-Simple Sequence Repeats (ISSR), Microsatellite-Primed Genomic Profiling Using Universal Primers

Inter-simple sequence repeat (ISSR) markers are highly polymorphic, relatively easy to develop, and inexpensive compared to other methods and have numerous applications. Importantly, the same ISSR primers can potentially be used universally across plant phylogenetic diversity. The basic technique of ISSRs is flexible and can be modified with options for implementation for a broad range of projects and budgets. Ranked in increasing order of technical demand and costs, these are manual agarose and manual polyacrylamide with silver staining and automated using fluorescently labeled primers and capillary electrophoresis. Overall manual agarose-based ISSRs are a sound, safe, easy, and low-cost method for reliably inferring plant genetic diversity. Here, we provide detailed protocols to undertake this fingerprinting method and provide guidance to the literature for the many options available for this technique.

Ophiotaenia echidis n. sp. (Cestoda: Proteocephalidae) from the saw-scaled viper, Echis carinatus sochureki Stemmler (Ophidia: Viperidae), one of the world's deadliest snakes, from the United Arab Emirates

Ophiotaenia echidis n. sp. (Cestoda: Proteocephalidae) is described from the intestine of one of the world's deadliest snakes, the saw-scaled viper Echis carinatus sochureki Stemmler (Ophidia: Viperidae) in the United Arab Emirates. The new species differs from other species of the non-monophyletic Ophiotaenia by the position of testes in two longitudinal lines on both sides of the uterus, and by the large size of an embryophore (diameter of 44–55 μm versus less than 40 μm in other species). Phylogenetic reconstructions based on lsrDNA and concatenated lsrDNA + COI datasets place the new species among proteocephalids from unrelated zoogeographical realms but mostly infecting venomous snakes. In all analyses, O. echidis n. sp. exhibited a strongly supported sister relationship with O. lapata Rambeloson, Ranaivoson et de Chambrier, 2012, a parasite of a pseudoxyrhophiid snake endemic to Madagascar. Despite a shared close evolutionary history between these taxa, morphological synapomorphies remain unclear, which impedes the erection of a new genus to accommodate them. A list of the 71 tapeworms of the former, non-monophyletic subfamily Proteocephalinae, parasitising snakes and lizards, including species inquirendae, and the phylogenetically closely related Thaumasioscolex didelphidis from opossum, with selected characteristics, is also provided, together with a checklist of helminth parasites reported from E. carinatus.

•

We describe the first proteocephalid cestode from the saw-scaled viper Echis carinatus.

•

The new ‘Ophiotaenia’ is typified by the position of testes and embryophore size.

•

The new species phylogenetically clusters with snake cestodes from unrelated regions.

•

The transmission mechanism remains unknown but likely involves some arthropod prey.

Mitochondrial DNA revealed the validation of Quasipaa robertingeri (Amphibia: Anura: Dicroglossidae) and its population genetic diversity

The spiny frog Quasipaa robertingeri is endemic to a narrow region of southwest China and its taxonomic validation is still controversial. Based on COI gene sequences of 110 individuals from seven populations of Q. robertingeri and its related species, we investigated the phylogenetic position and population genetic structure of the species. Phylogenetic analyses indicated that Q. robertingeri was deeply genetically separated from its closely related species Q. boulengeri, indicating the validation of the species. All samples of Q. robertingeri were clustered into two divergent lineages. Haplotype network, AMOVA, and genetic distance estimations also supported the separations of the two groups. Neutrality tests indicated that one lineage has been likely independently experienced a recent population expansion, leading to a secondary contact area between the two lineages.

Paramyxovirus circulation in bat species from French Guiana

Bats are recognized as reservoirs of numerous viruses. Among them, paramyxoviruses, for example, Hendra and Nipah viruses, are highly pathogenic to humans. Nothing is known regarding the circulation of this viral family in bats from French Guiana. To search for the presence of paramyxoviruses in this territory, 103 bats of seven different species were sampled and screened using a molecular approach. Four distinct paramyxovirus sequences were detected from three bat species (Desmodus rotundus, Carollia perspicillata, and Pteronotus alitonus) at high prevalence rates. In D. rotundus, two types of paramyxovirus co-circulate, with most of the bats co-infected. The phylogenetic analysis of these sequences revealed that three of them were closely related to previously characterized sequences from D. rotundus, C. perspicillata, and P. parnellii from Brazil and Costa Rica. The fourth sequence, identified in D. rotundus, was closely related to the one detected in P. alitonus in French Guiana and to previously described sequences detected in P. parnellii in Costa Rica. All paramyxovirus sequences detected in this study are close to the Jeilongvirus genus. Altogether, our results and those of previous studies indicate a wide geographical distribution of these paramyxoviruses (from Central to South America) and suggest potential cross-species transmissions of paramyxoviruses between two different bat families: Mormoopidae (P. alitonus) and Phyllostomidae (D. rotundus). In addition, their closeness to paramyxoviruses identified in rodents emphasizes the need to investigate the role of these animals as potential reservoirs or incidental hosts. Finally, the high prevalence rates of some paramyxoviruses in certain bat species, associated with the presence of large bat colonies and, in some cases, their potential proximity with humans are all parameters that can contribute to the risk of cross-species transmission between bat species and to the emergence of new paramyxoviruses in humans, a risk that deserves further investigation.

In silico characterization of bacterial chitinase: illuminating its relationship with archaeal and eukaryotic cousins

Background

Chitin is one of the most abundant biopolymers on Earth, only trailing second after cellulose. The enzyme chitinase is responsible for the degradation of chitin. Chitinases are found to be produced by wide range of organisms ranging from archaea to higher plants. Though chitin is a major component of fungal cell walls and invertebrate exoskeletons, bacterial chitinase can be industrially generated at low cost, in facile downstream processes at high production rate. Microbial chitinases are more stable, active, and economically practicable compared to the plant- and animal-derived enzymes.

Results

In the present study, computationally obtained results showed functional characteristics of chitinase with particular emphasis on bacterial chitinase which is fulfilling all the required qualities needed for commercial production. Sixty-two chitinase sequences from four different groups of organisms were collected from the RCSB Protein Data Bank. Considering one suitable exemplary sequence from each group is being compared with others. Primary, secondary, and tertiary structures are determined by in silico models. Different physical parameters, viz., pI, molecular weight, instability index, aliphatic index, GRAVY, and presence of functional motifs, are determined, and a phylogenetic tree has been constructed to elucidate relationships with other groups of organisms.

Conclusions

This study provides novel insights into distribution of chitinase among four groups and their characterization. The results represent valuable information toward bacterial chitinase in terms of the catalytic properties and structural features, can be exploited to produce a range of chitin-derived products.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-021-00121-6.

Determination of phylogenetic relationships among methicillin-resistant Staphylococcus aureus recovered from infected humans and Companion Animals

Companion animals carry different microorganism of severely public health hazard for human; the kindness relation and contact between humans and companion animals may the route in the transmission of most zoonotic bacteria, including Methicillin-Resistant Staphylococcus aureus (MRSA). Therefore, the current study investigate the companion animals mainly dogs and cat as a reservoir for MRSA and the genetic similarity between the recovered strains of MRSA from such companion animals and their owners. One hundred swabs were collected under aseptic condition from companion animals and seventy swabs were collected from nasal and soft tissue of the infected owners in contact. All samples were examined with standard microbiological techniques, antimicrobial sensitivity, molecular typing and genetic finger printing using RAPD-PCR to determine the genetic finger printing of the recovered strains from humans and companion animals. The prevalence of the MRSA was higher in dog's swabs than human swabs. Dog swabs showed a rate of (44.4%), cat's revealed (27.3%), while the owner swabs could detect (42.8%). The antibiotics profiles were 69.2% and all MRSA strains were positive for mecA gene (100%), while only 25 strains (38.5%) were positive for Panton Valentine Leukocidin (PVL gene). Phylogenetic tree revealed 4 clusters with complete genetic relatedness and higher identity between the strains recovered from humans and companion animals. Our results revealed that there is great similarity between the recovered strains, indicating that pets play an important role in colonization and transmitting MRSA to humans, and vice versa.

Characterization of the complete mitochondrial genome of Macrotocinclus affinis (Siluriformes; Loricariidae) and phylogenetic studies of Siluriformes

The vertebrate mitochondrial genome is typically circular molecules made up of 14,000 to 16,000 bp, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (12 s rRNA and 16 s rRNA) and a control region. Compared with nuclear DNA, mitochondrial DNA has a higher mutation rate, so it is one of the most effective and reliable molecular markers in fish phylogeny. Macrotocinclus affinis was the only species in Macrotocinclus (it was classified as Otocinclus in the past) and currently lacks genetic information. Most of the current researches are based on the mitochondrial Cytb gene and RAG1 and RAG2 nuclear genes to study the phylogenetic analysis of Siluriformes. So, the study provides the characteristic features of the Macrotocinclus affinis mitochondrial genome and this is the first time that the phylogenetic relationship of Siluriformes has been reconstructed based on COI. We aimed to sequence the entire mitochondrial genome of Macrotocinclus affinis using conventional PCR techniques and to clarify its phylogenetic status in Siluriformes by using the COI sequence of mitochondria. In this study, we sequenced the whole mitochondrial genome of this species yielding a 16,632 bp circular assembly composed of the typical vertebrate mitochondrial features. It contains 13 protein-coding genes, two rRNA genes, 22 tRNA genes, a putative control region, and one origin of replication on the light-strand. The overall base composition includes A (30.07%), T (24.43%), C (29.43%) and G (16.01%). The genome composition is A + T biased (54.5%), and exhibits AT-skew (0.1036) and GC-skew (-0.2962). Moreover, the 13 PCGs encode 3850 amino acids in total. The result of the phylogenetic tree supports Macrotocinclus affinis has a closest relationship with Otocinclus cf. hoppei far. These results will help to understand the characteristics of the mitochondrial genome of Macrotocinclus affinis and provide molecular basis for the evolutionary relationship of Loricariidae.

Multilocus Variable-Number Tandem-Repeat Analysis for Geno-typing of Escherichia coli Strains Isolated from Hospital Wastewater, Tehran, Iran

Background:

Escherichia coli is one of the most frequent causes of many common bacterial infections. As a potential reservoir, hospital wastewater is considered for the dissemination of bacterial pathogens such as E. coli. Therefore, research on hospital waste’s bacteria by low-cost, rapid and easy molecular typing methods such as multilocus variable-number tandem-repeat analysis (MLVA) can be helpful for the study of epidemics.

Methods:

E. coli strains were isolated from hospital wastewater sources in Tehran, Iran, over a 24-month sampling period (Jun 2014- Jun 2016) and identified by standard bacteriological methods. The diversity of repeated sequences of seven variable-number tandem-repeat (VNTR) loci was studied by MLVA method base on polymerase chain reaction (PCR).

Results:

Overall, 80 E. coli isolates were discriminated into 51 different genotypes. Analysis of the MLVA profiles using a minimum spanning tree (MST) algorithm showed two clonal complexes with 71 isolates and only nine isolates were stayed out of clonal complexes in the form of a singleton. High genotypic diversity was seen among E. coli strains isolated from hospital wastewaters; however, a large number of isolates showed a close genetic relationship.

Conclusion:

MLVA showed to be a rapid, inexpensive and useful tool for the analysis of the phylogenetic relationships between E. coli strains under the study.

Plastomes of eight Ligusticum species: characterization, genome evolution, and phylogenetic relationships

BACKGROUND

The genus Ligusticum consists of approximately 60 species distributed in the Northern Hemisphere. It is one of the most taxonomically difficult taxa within Apiaceae, largely due to the varied morphological characteristics. To investigate the plastome evolution and phylogenetic relationships of Ligusticum, we determined the complete plastome sequences of eight Ligusticum species using a de novo assembly approach.

RESULTS

Through a comprehensive comparative analysis, we found that the eight plastomes were similar in terms of repeat sequence, SSR, codon usage, and RNA editing site. However, compared with the other seven species, L. delavayi exhibited striking differences in genome size, gene number, IR/SC borders, and sequence identity. Most of the genes remained under the purifying selection, whereas four genes showed relaxed selection, namely ccsA, rpoA, ycf1, and ycf2. Non-monophyly of Ligusticum species was inferred from the plastomes and internal transcribed spacer (ITS) sequences phylogenetic analyses.

CONCLUSION

The plastome tree and ITS tree produced incongruent tree topologies, which may be attributed to the hybridization and incomplete lineage sorting. Our study highlighted the advantage of plastome with mass informative sites in resolving phylogenetic relationships. Moreover, combined with the previous studies, we considered that the current taxonomy system of Ligusticum needs to be improved and revised. In summary, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of Ligusticum species.

Sequence analyses of RT-PCR products obtained from seven infected leaf samples revealed existence of three potyvirus species in Indian narcissus (Narcissus tazetta L.)

Potyvirus species associated with yellow leaf stripe disease of Indian narcissus (Narcissus tazetta L.) var. Paperwhite has been studied by sequence analyses of ~ 1.5 kb genomic fragments obtained from seven RT-PCR amplifications of infected samples. Sequence analysis revealed the occurrence of three potyvirus species: cyrtanthus elatus virus-A (CEVA: KF430815, KF430816, KM066973, KM066974); narcissus yellow stripe virus (NYSV: KM066972, JQ686724) and narcissus degeneration virus (NDV: MK572806). The existence of three potyvirus species: CEVA, NYSV and NDV are being reported in Indian narcissus.

Repeated intercontinental migrations and recurring hybridizations characterise the evolutionary history of yew (Taxus L.)

The genus Taxus (Taxaceae) consists of 16 genetically well-defined lineages that are predominantly distributed across the Northern hemisphere. We investigated its biogeographic origin and evolutionary history by sampling 13 chloroplast gene sequences, the nuclear internal transcribed spacers (ITS) and NEEDLY sequences for all 16 lineages. We applied Maximum Parsimony and Bayesian Inference analyses to infer their phylogenetic relationships, time-calibrated phylogenies using BEAST and inferred the ancestral area of occupancy with BioGeoBEARS. We found strong evidence for the hybrid origin of three lineages and dated these events to a rather narrow time window of 6.8-4.9 million years ago (Mya). The dated phylogenies inferred an Upper Cretaceous origin of the genus, with the extant lineages diversifying in North America much later during the Oligocene/early Miocene. Repeated migrations via the Bering land bridge to Eurasia and back were further inferred, with the return to North America as a possible result of vicariance. The diversification in Eurasia (from ~8 Mya onwards) coincided with the orogeny of the Hengduan Mountains, the intensification of the East Asian summer monsoon and the occupancy of ecological niches by lineages that experienced secondary contacts and hybridizations in the Hengduan Mountains and Qinling Mountain, especially around the Sichuan basin. We provide a hypothesis for the evolution of extant lineages of Taxus, a genus with an old and complex evolutionary history. The study highlights that the history of complex species can be unravelled with a careful dissection of phylogenetic signals.

The mitochondrial genome contribution to the phylogeny and identification of Metarhizium species and strains

The genus Metarhizium is composed of entomopathogenic fungal biological control agents (BCAs) used for invertebrate pest control. The phylogenetic relationships of species within this genus are still under scrutiny as several cryptic species can be found. In this work, the mitochondrial (mt) genome of Metarhizium brunneum ARSEF 4556 was fully sequenced and a comparative genome analysis was conducted with 7 other available mt genomes, belonging to 5 Metarhizium species: M. anisopliae, M. brunneum, M. robertsii, M. guizhouense and M. majus. Results showed that Metarhizium demonstrates greater conserved stability than other fungal mt genomes. Furthermore, this analysis located 7 diverse regions in both intergenic domains and gene fragments which were ideal for species/strain discrimination. The sequencing of these regions revealed several SNPs among 38 strains tested, 11 of which were uncharacterized. Single gene phylogenies presented variable results which may be used further for intra-species discrimination. Phylogenetic trees based on the concatenation of mt domains and the nuclear ITS1-5.8S-ITS2 region showed discrimination of the species studied and allowed the identification of uncharacterized strains. These were mostly placed within species M. anisopliae and M. brunneum. Five strains clustered together in a clade related to M. brunneum, suggesting that they comprise a cryptic species.

Comparative microscopic examination of scales in 21 clupeid species from the Caspian Sea and the Indo-Pacific regions

The clupeid fishes are large filter-feeding organisms and have diverse groups of trophic guilds and habitats. They are known for their special body scale, which is unique among the other teleost fishes. Here, we examined microscopic characteristics of scales in 21 clupeid species from four marine resources, including the Caspian Sea, the Persian Gulf and the Gulf of Oman (Indian Ocean) and South China Sea (Haian Islands and Shanwei, Pacific Ocean) to evaluate if the scale features could properly contribute to understanding the taxonomic relationships in these fishes. Typically, the studied scales were morphologically categorized into two types; circular (i.e., true circular, cordate, discoidal) and pentagonal. Also, it was found that the number and the orientation of radii and continuous striae, the presence of pores at the posterior field, overall shape, striations on the posterior field, and scale sizes could be used as suitable features to classify the clupeids in species and genus levels. The result of morphological analysis based on the scale characteristics was largely consistent with the phylogenetic relationships of the studied species. This highlighted that morphological characteristics of the clupeid scale could adequately contribute to the understanding of taxonomic relationships in this large group.

Characterization of the mitochondrial genomes of two toads, Anaxyrus americanus (Anura: Bufonidae) and Bufotes pewzowi (Anura: Bufonidae), with phylogenetic and selection pressure analyses

Mitogenomes are useful in analyzing phylogenetic relationships and also appear to influence energy metabolism, thermoregulation and osmoregulation. Much evidence has accumulated for positive selection acting on mitochondrial genes associated with environmental adaptation. Hence, the mitogenome is a likely target for environmental selection. The family Bufonidae (true toads) has only nine complete and four partial mitogenomes published compared to the 610 known species of this family. More mitogenomes are needed in order to obtain a clearer understanding of the phylogenetic relationships within Bufonidae that are currently controversial. To date, no mitogenomes have been reported from the genera Anaxyrus and Bufotes. Anaxyrus americanus can live in low temperature environments and Bufotes pewzowi can live in high salinity environments. We sequenced the mitogenomes of these two species to discuss the phylogenetic relationships within Bufonidae and the selection pressures experienced by specimens living in low temperature or saline environments. Like other toads, the circular mitogenomes of both species contained the typical 37 genes. Anaxyrus americanus had the highest A+T content of the complete mitogenome among the Bufonidae. In addition, A. americanus showed a negative AT-skew in the control region, whereas Bufotes pewzowi showed a positive AT-skew. Additionally, both toad species had unique molecular features in common: an ND1 gene that uses TTG as the start codon, an extra unpaired adenine (A) in the anticodon arm of trnS (AGY), and the loss of the DHU loop in trnC. The monophyly of Bufonidae was corroborated by both BI and ML trees. An analysis of selective pressure based on the 13 protein coding genes was conducted using the EasyCodeML program. In the branch model analysis, we found two branches of A. americanus and Bufotes pewzowi that were under negative selection. Additionally, we found two positively selected sites (at positions 115 and 119, BEB value > 0.90) in the ND6 protein in the site model analysis. The residue D (119) was located only in A. americanus and may be related to adaptive evolution in low temperature environments. However, there was no evidence of a positively selected site in Bufotes pewzowi in this study.

Narcinecotyle longifilamentus n. gen., n. sp. (Monogenea: Hexabothriidae), gill parasite of the numbfish Narcine entemedor (Torpediniformes: Narcinidae) from the Mexican Pacific coast

A new genus and species of Hexabothriidae Price, 1942 is described from specimens infecting the gill of the numbfish Narcine entemedor Jordan and Starks from the Pacific coast of Mexico. In addition to the new taxon described here, species of 8 genera also display symmetrical haptors and have eggs with two polar filaments. However, they differ in the combination of the following features: distal portion of the male copulatory organ tubular, unarmed and proximally dilated confining an internal coiled duct, as well as dorsal origin of haptoral appendix. Based on molecular data derived from 3 loci, the mitochondrial Cytochrome C oxidase subunit 1 as well as the nuclear ribosomal 18S and 28S, the new species was found nested within Hexabothriidae together with the other 4 genera with representatives in GenBank. This is the first species of Hexabothriidae reported from a species of Narcinidae.

The first mitochondrial genome for Phaudidae (Lepidoptera) with phylogenetic analyses of Zygaenoidea

Phauda flammans Walker belongs to Phaudidae (Lepidoptera), which is a holometabolous and leaf-eating pest that harms trees. So far, there is no mitochondrial (mt) genome reported of Phaudidae. Herein, we sequenced and annotated the complete mt genome of P. flammans representing the first mt genome of Phaudidae and predicted the secondary structures of its RNAs in this study. This mt genome is 15470 bp long consisting of 13 protein coding genes (PCGs), 22 tRNAs, 2 rRNAs and the control region, which are usually conserved in insects. Most PCGs used the standard ATN start codons and complete TAA/TAG termination codons. Almost all of tRNA genes exhibited cloverleaf secondary structures except that the dihydorouridine (DHU) arm of tRNA^Ser(AGN) was absent. The phylogenetic analyses using both Bayesian inference (BI) and maximum likelihood (ML) methods all supported that Phaudidae was a single family being the sister group to Zygaenidae. More mt genomes are needed to better understand the phylogenetic relationships within Zygaenoidea in the future.

Characterization of wheat homeodomain-leucine zipper family genes and functional analysis of TaHDZ5-6A in drought tolerance in transgenic Arabidopsis

BACKGROUND

Many studies in Arabidopsis and rice have demonstrated that HD-Zip transcription factors play important roles in plant development and responses to abiotic stresses. Although common wheat (Triticum aestivum L.) is one of the most widely cultivated and consumed food crops in the world, the function of the HD-Zip proteins in wheat is still largely unknown.

RESULTS

To explore the potential biological functions of HD-Zip genes in wheat, we performed a bioinformatics and gene expression analysis of the HD-Zip family. We identified 113 HD-Zip members from wheat and classified them into four subfamilies (I-IV) based on phylogenic analysis against proteins from Arabidopsis, rice, and maize. Most HD-Zip genes are represented by two to three homeoalleles in wheat, which are named as TaHDZX_ZA, TaHDZX_ZB, or TaHDZX_ZD, where X denotes the gene number and Z the wheat chromosome on which it is located. TaHDZs in the same subfamily have similar protein motifs and intron/exon structures. The expression profiles of TaHDZ genes were analysed in different tissues, at different stages of vegetative growth, during seed development, and under drought stress. We found that most TaHDZ genes, especially those in subfamilies I and II, were induced by drought stress, suggesting the potential importance of subfamily I and II TaHDZ members in the responses to abiotic stress. Compared with wild-type (WT) plants, transgenic Arabidopsis plants overexpressing TaHDZ5-6A displayed enhanced drought tolerance, lower water loss rates, higher survival rates, and higher proline content under drought conditions. Additionally, the transcriptome analysis identified a number of differentially expressed genes between 35S::TaHDZ5-6A transgenic and wild-type plants, many of which are involved in stress response.

CONCLUSIONS

Our results will facilitate further functional analysis of wheat HD-Zip genes, and also indicate that TaHDZ5-6A may participate in regulating the plant response to drought stress. Our experiments show that TaHDZ5-6A holds great potential for genetic improvement of abiotic stress tolerance in crops.

Genetic structure and geographical variation of Bithynia siamensis goniomphalos sensu lato (Gastropoda: Bithyniidae), the snail intermediate host of Opisthorchis viverrini sensu lato (Digenea: Opisthorchiidae) in the Lower Mekong Basin revealed by mitochondrial DNA sequences

The freshwater snail Bithynia siamensis goniomphalos sensu lato is widely distributed in the Lower Mekong Basin where it acts as the first intermediate host of the liver fluke Opisthorchis viverrini, a group 1 carcinogen causing cholangiocarcinoma. This study explores the genetic structure and geographical variation of B. s. goniomphalos from eight previously studied catchments and eight new catchments. These catchments belong to five previously studied catchment systems and one new catchment system (Tonlesap) in the Lower Mekong Basin. Two new catchment systems, Prachin Buri and Bang Pakong from eastern and central Thailand, respectively, were also examined. We collected 289 specimens of B. s. goniomphalos from 15 previously studied localities and 18 new localities in Thailand, Lao PDR (People's Democratic Republic), and Cambodia. The mitochondrial cytochrome c oxidase subunit 1 and 16S ribosomal DNA sequences were used to determine genetic variation. Classification of haplotypes specified 100 at the cox1 locus and 15 at the rrnL locus. Comparison between 16 catchment populations found significant genetic differences (Ф_ST) between all populations. The phylogenetic tree and haplotype network analyses classified B. s. goniomphalos into three evolutionary lineages (lineage I-III). Lineage I contained B. s. goniomphalos from the Mekong, Chi, Mun, Prachin Buri and Bang Pakong catchments in Thailand, including the Nam Ngum catchment in Lao PDR. Lineage II contained all specimens from the Tonlesap catchment, whereas lineage III contained specimens from the Mekong and Sea Bang Heang catchments in Thailand and Lao PDR, respectively. Interestingly, Bithynia siamensis siamensis was placed between lineages I and II of B. s. goniomphalos. This study supports the hypothesis that B. s. goniomphalos is a species complex containing at least three distinct evolutionary lineages in the Lower Mekong Basin, and that comprehensive molecular genetic analyses need to be conducted to further our understanding of the evolutionary and systematic relationships of these Bithynia snail taxa.

Comparison of chloroplast genomes of Gynura species: sequence variation, genome rearrangement and divergence studies

Background

Some Gynura species have been reported to be natural anti-diabetic plants. Improvement of their traits towards application relies on hybridization. Clearly, phylogenetic relationships could optimize compatible hybridizations. For flowerings plants, chloroplast genomes have been used to solve many phylogenetic relationships. To date, the chloroplast genome sequences of 4 genera of the tribe Senecioneae have been uploaded to GenBank. The internal relationships within the genus Gynura and the relationship of the genus Gynura with other genera in the tribe Senecioneae need further research.

Results

The chloroplast genomes of 4 Gynura species were sequenced, assembled and annotated. In comparison with those of 12 other Senecioneae species, the Gynura chloroplast genome features were analysed in detail. Subsequently, differences in the microsatellite and repeat types in the tribe were found. From the comparison, it was found that IR expansion and contraction are conserved in the genera Gynura, Dendrosenecio and Ligularia. Compared to other regions on the chloroplast genome, the region from 25,000 to 50,000 bp was not conserved. Seven ndh genes in this region are under purifying selection, with small changes in amino acids. The whole chloroplast genome sequences of 16 Senecioneae species were used to build a phylogenetic tree. Based on the oldest Artemisia pollen fossil, the divergence time was estimated.

Conclusions

Sequencing the chloroplast genomes of 4 Gynura species helps us to solve many problems. The phylogenetic relationships and divergence time among 4 Gynura and 16 Senecioneae species were evaluated by comparing their chloroplast genomes. The phylogenetic relationship of the genera Gynura and Ligularia was different from that observed previous work. In a previous phylogenetic tree, the genus Ligularia belonged to the Tussilagininae subtribe, which was in a lineage that diverged earlier than other genera. Further morphology and genome-wide analyses are needed to clarify the genus relationships.