Metagenomic profiling of gut microbiota in Fall Armyworm (Spodoptera frugiperda) larvae fed on different host plants.

BACKGROUND: The fall armyworm (FAW, Spodoptera frugiperda) is a polyphagous pest known for causing significant crop damage. The gut microbiota plays a pivotal role in influencing the biology, physiology and adaptation of the host. However, understanding of the taxonomic composition and functional characteristics of the gut microbiota in FAW larvae fed on different host plants remains limited. METHODS: This study utilized metagenomic sequencing to explore the structure, function and antibiotic resistance genes (ARGs) of the gut microbiota in FAW larvae transferred from an artificial diet to four distinct host plants: maize, sorghum, tomato and pepper. RESULTS: The results demonstrated significant variations in gut microbiota structure among FAW larvae fed on different host plants. Firmicutes emerged as the dominant phylum, with Enterococcaceae as the dominant family and Enterococcus as the prominent genus. Notably, Enterococcus casseliflavus was frequently observed in the gut microbiota of FAW larvae across host plants. Metabolism pathways, particularly those related to carbohydrate and amino acid metabolism, played a crucial role in the adaptation of the FAW gut microbiota to different host plants. KEGG orthologs associated with the regulation of the peptide/nickel transport system permease protein in sorghum-fed larvae and the 6-phospho-ß-glucosidase gene linked to glycolysis/gluconeogenesis as well as starch and sucrose metabolism in pepper-fed larvae were identified. Moreover, the study identified the top 20 ARGs in the gut microbiota of FAW larvae fed on different host plants, with the maize-fed group exhibiting the highest abundance of vanRC. CONCLUSIONS: Our metagenomic sequencing study reveals significant variations in the gut microbiota composition and function of FAW larvae across diverse host plants. These findings underscore the intricate co-evolutionary relationship between hosts and their gut microbiota, suggesting that host transfer profoundly influences the gut microbiota and, consequently, the adaptability and pest management strategies for FAW.

Inferring Historical Introgression with Deep Learning.

Resolving phylogenetic relationships among taxa remains a challenge in the era of big data due to the presence of genetic admixture in a wide range of organisms. Rapidly developing sequencing technologies and statistical tests enable evolutionary relationships to be disentangled at a genome-wide level, yet many of these tests are computationally intensive and rely on phased genotypes, large sample sizes, restricted phylogenetic topologies, or hypothesis testing. To overcome these difficulties, we developed a deep learning-based approach, named ERICA, for inferring genome-wide evolutionary relationships and local introgressed regions from sequence data. ERICA accepts sequence alignments of both population genomic data and multiple genome assemblies, and efficiently identifies discordant genealogy patterns and exchanged regions across genomes when compared with other methods. We further tested ERICA using real population genomic data from Heliconius butterflies that have undergone adaptive radiation and frequent hybridization. Finally, we applied ERICA to characterize hybridization and introgression in wild and cultivated rice, revealing the important role of introgression in rice domestication and adaptation. Taken together, our findings demonstrate that ERICA provides an effective method for teasing apart evolutionary relationships using whole genome data, which can ultimately facilitate evolutionary studies on hybridization and introgression.

Biochemical Property Based Positional Matrix: A New Approach Towards Genome Sequence Comparison.

The growth of the genome sequence has become one of the emerging areas in the study of bioinformatics. It has led to an excessive demand for researchers to develop advanced methodologies for evolutionary relationships among species. The alignment-free methods have been proved to be more efficient and appropriate related to time and space than existing alignment-based methods for sequence analysis. In this study, a new alignment-free genome sequence comparison technique is proposed based on the biochemical properties of nucleotides. Each genome sequence can be distributed in four parameters to represent a 21-dimensional numerical descriptor using the Positional Matrix. To substantiate the proposed method, phylogenetic trees are constructed on the viral and mammalian datasets by applying the UPGMA/NJ clustering method. Further, the results of this method are compared with the results of the Feature Frequency Profiles method, the Positional Correlation Natural Vector method, the Graph-theoretic method, the Multiple Encoding Vector method, and the Fuzzy Integral Similarity method. In most cases, it is found that the present method produces more accurate results than the prior methods. Also, in the present method, the execution time for computation is comparatively small.

Multichromosomal mitochondrial genome of Punica granatum: comparative evolutionary analysis and gene transformation from chloroplast genomes.

BACKGROUND: Punica granatum is a fundamentally important fruit tree that has important economic, medicinal and ornamental properties. At present, there are few reports on the mitochondrial genome of pomegranate. Hence, in this study the P. granatum mitogenome was sequenced and assembled to further understanding of organization, variation, and evolution of mitogenomes of this tree species. RESULTS: The genome structure was multi-chromosomes with seven circular contigs, measuring 382,774 bp in length with a 45.91% GC content. It contained 74 genes, including 46 protein-coding genes, 25 tRNA genes, and three rRNA genes. There were 188 pairs of dispersed repeats with lengths of 30 or greater, primarily consisting of reverse complementary repeats. The mitogenome analysis identified 114SSRs and 466 RNA editing sites. Analyses of codon usage, nucleotide diversity and gene migration from chloroplast to mitochondrial were also conducted. The collinear and comparative analysis of mitochondrial structures between P. granatum and its proximal species indicated that P. granatum 'Taishanhong' was closely related to P. granatum 'Qingpitian' and Lagerstroemia indica. Phylogenetic examination based on the mitogenome also confirmed the evolutionary relationship. CONCLUSION: The results offered crucial information on the evolutionary biology of pomegranate and highlighted ways to promote the utilization of the species' germplasm.

Range-wide evolutionary relationships and historical demography of brown bears (Ursus arctos) revealed by whole-genome sequencing of isolated central Asian populations.

Phylogeographic studies uncover hidden pathways of divergence and inform conservation. Brown bears (Ursus arctos) have one of the broadest distributions of all land mammals, ranging from Eurasia to North America, and are an important model for evolutionary studies. Although several whole genomes were available for individuals from North America, Europe and Asia, limited whole-genome data were available from Central Asia, including the highly imperilled brown bears in the Gobi Desert. To fill this knowledge gap, we sequenced whole genomes from nine Asian brown bears from the Gobi Desert of Mongolia, Northern Mongolia and the Himalayas of Pakistan. We combined these data with published brown bear sequences from Europe, Asia and North America, as well as other bear species. Our goals were to determine the evolutionary relationships among brown bear populations worldwide, their genetic diversity and their historical demography. Our analyses revealed five major lineages of brown bears based on a filtered set of 684,081 single nucleotide polymorphisms. We found distinct evolutionary lineages of brown bears in the Gobi, Himalayas, northern Mongolia, Europe and North America. The lowest level of genetic diversity and the highest level of inbreeding were found in Pakistan, the Gobi Desert and Central Italy. Furthermore, the effective population size (Ne ) for all brown bears decreased over the last 70,000 years. Our results confirm the genetic distinctiveness and ancient lineage of brown bear subspecies in the Gobi Desert of Mongolia and the Himalayas of Pakistan and highlight their importance for conservation.

Peach DELLA Protein PpeDGYLA Is Not Degraded in the Presence of Active GA and Causes Dwarfism When Overexpressed in Poplar and Arabidopsis.

Controlling the tree size of fruit species such as peach can reduce the amount of labor and input needed for orchard management. The phytohormone gibberellin (GA) positively regulates tree size by inducing degradation of the GA signaling repressor DELLA. The N-terminal DELLA domain in this protein is critical for its GA-dependent interaction with the GA receptor GID1 and the resulting degradation of the DELLA protein, which allows for growth-promoting GA signaling. In this study, a DELLA family member, PpeDGYLA, contains a DELLA domain but has amino acid changes in three conserved motifs (DELLA into DGYLA, LEQLE into LERLE, and TVHYNP into AVLYNP). In the absence or presence of GA3, the PpeDGYLA protein did not interact with PpeGID1c and was stable in 35S-PpeDGYLA peach transgenic callus. The overexpression of PpeDGYLA in both polar and Arabidopsis showed an extremely dwarfed phenotype, and these transgenic plants were insensitive to GA3 treatment. PpeDGYLA could interact with PpeARF6-1 and -2, supposed growth-promoting factors. It is suggested that the changes in the DELLA domain of PpeDGYLA may, to some extent, account for the severe dwarf phenotype of poplar and Arabidopsis transgenic plants. In addition, our study showed that the DELLA family contained three clades (DELLA-like, DELLA, and DGLLA). PpeDGYLA clustered into the DGLLA clade and was expressed in all of the analyzed tissues. These results lay the foundation for the further study of the repression of tree size by PpeDGYLA.

Genome-Wide Analysis of the Polygalacturonase Gene Family Sheds Light on the Characteristics, Evolutionary History, and Putative Function of Akebia trifoliata.

Polygalacturonase (PG) is one of the largest families of hydrolytic enzymes in plants. It is involved in the breakdown of pectin in the plant cell wall and even contributes to peel cracks. Here, we characterize PGs and outline their expression profiles using the available reference genome and transcriptome of Akebia trifoliata. The average length and exon number of the 47 identified AktPGs, unevenly assigned on 14 chromosomes and two unassembled contigs, were 5399 bp and 7, respectively. The phylogenetic tree of 191 PGs, including 47, 57, 51, and 36 from A. trifoliata, Durio zibethinus, Actinidia chinensis, and Vitis vinifera, respectively, showed that AktPGs were distributed in all groups except group G and that 10 AktPGs in group E were older, while the remaining 37 AktPGs were younger. Evolutionarily, all AktPGs generally experienced whole-genome duplication (WGD)/segmental repeats and purifying selection. Additionally, the origin of conserved domain III was possibly associated with a histidine residue (H) substitute in motif 8. The results of both the phylogenetic tree and expression profiling indicated that five AktPGs, especially AktPG25, could be associated with the cracking process. Detailed information and data on the PG family are beneficial for further study of the postharvest biology of A. trifoliata.

Genome-Wide Analysis of the LATERAL ORGAN BOUNDARIES Domain (LBD) Members in Alfalfa and the Involvement of MsLBD48 in Nitrogen Assimilation.

The LATERAL ORGAN BOUNDARIES DOMAIN (LBD) proteins, a transcription factor family specific to the land plants, have been implicated in multiple biological processes including organ development, pathogen response and the uptake of inorganic nitrogen. The study focused on LBDs in legume forage Alfalfa. The genome-wide analysis revealed that in Alfalfa 178 loci across 31 allelic chromosomes encoded 48 unique LBDs (MsLBDs), and the genome of its diploid progenitor M. sativa spp. Caerulea encoded 46 LBDs. Synteny analysis indicated that the expansion of AlfalfaLBDs was attributed to the whole genome duplication event. The MsLBDs were divided into two major phylogenetic classes, and the LOB domain of the Class I members was highly conserved relative to that of the Class II. The transcriptomic data demonstrated that 87.5% of MsLBDs were expressed in at least one of the six test tissues, and Class II members were preferentially expressed in nodules. Moreover, the expression of Class II LBDs in roots was upregulated by the treatment of inorganic nitrogen such as KNO3 and NH4Cl (0.3 mM). The overexpression of MsLBD48, a Class II member, in Arabidopsis resulted in growth retardance with significantly declined biomass compared with the non-transgenic plants, and the transcription level of the genes involved in nitrogen uptake or assimilation, including NRT1.1, NRT2.1, NIA1 and NIA2 was repressed. Therefore, the LBDs in Alfalfa are highly conserved with their orthologs in embryophytes. Our observations that ectopic expression of MsLBD48 inhibited Arabidopsis growth by repressing nitrogen adaption suggest the negative role of the transcription factor in plant uptake of inorganic nitrogen. The findings imply the potential application of MsLBD48 in Alfalfa yield improvement via gene editing.

Comprehensive Analysis of bHLH Transcription Factors in Ipomoea aquatica and Its Response to Anthocyanin Biosynthesis.

The basic helix-loop-helix (bHLH) proteins compose one of the largest transcription factor (TF) families in plants, which play a vital role in regulating plant biological processes including growth and development, stress response, and secondary metabolite biosynthesis. Ipomoea aquatica is one of the most important nutrient-rich vegetables. Compared to the common green-stemmed I. aquatica, purple-stemmed I. aquatica has extremely high contents of anthocyanins. However, the information on bHLH genes in I. aquatica and their role in regulating anthocyanin accumulation is still unclear. In this study, we confirmed a total of 157 bHLH genes in the I. aquatica genome, which were classified into 23 subgroups according to their phylogenetic relationship with the bHLH of Arabidopsis thaliana (AtbHLH). Of these, 129 IabHLH genes were unevenly distributed across 15 chromosomes, while 28 IabHLH genes were spread on the scaffolds. Subcellular localization prediction revealed that most IabHLH proteins were localized in the nucleus, while some were in the chloroplast, extracellular space, and endomembrane system. Sequence analysis revealed conserved motif distribution and similar patterns of gene structure within IabHLH genes of the same subfamily. Analysis of gene duplication events indicated that DSD and WGD played a vital role in the IabHLH gene family expansion. Transcriptome analysis showed that the expression levels of 13 IabHLH genes were significantly different between the two varieties. Of these, the IabHLH027 had the highest expression fold change, and its expression level was dramatically higher in purple-stemmed I. aquatica than that in green-stemmed I. aquatica. All upregulated DEGs in purple-stemmed I. aquatica exhibited the same expression trends in both qRT-PCR and RNA-seq. Three downregulated genes including IabHLH142, IabHLH057, and IabHLH043 determined by RNA-seq had opposite expression trends of those detected by qRT-PCR. Analysis of the cis-acting elements in the promoter region of 13 differentially expressed genes indicated that light-responsive elements were the most, followed by phytohormone-responsive elements and stress-responsive elements, while plant growth and development-responsive elements were the least. Taken together, this work provides valuable clues for further exploring IabHLH function and facilitating the breeding of anthocyanin-rich functional varieties of I. aquatica.

Farmland nutrient pollution and its evolutionary relationship with plantation economic development in China.

Contradiction between growing plantation economic demand and agro-ecological degradation has always restricted sustainable development of agricultural countries. This study applied the unit inventory analysis to evaluate the productions and discharges of farmland non-point source (FNPS) nitrogen (TN) and phosphorus (TP) among China's nine national-level agricultural districts over 1999-2019. On this basis, we quantified the evolutionary relationship between plantation economic output and FNPS pollution based on optimal regression fitting. The results showed that over 1999-2019, farmland cumulative TN and TP discharges for the whole China were approximately 15807 × 104 t and 1312 × 104 t, with prominent district heterogeneity. According to FNPS discharge magnitudes, China's agricultural districts can be classified into three categories: high, moderate and slight discharge zones. Huang-Huai-Hai Plain and Middle-lower Yangtze Plain were identified as the main severely-polluted districts. Mineral fertilizer is the primary contributor to FNPS pollution. Annual FNPS load showed a trend of increasing followed by decreasing, and the peak interval was recorded in 2014-2016. Spatiotemporal dynamics in FNPS discharge intensities were disparate from that in discharge magnitudes. SC has the highest TN discharge intensity, with an annual average intensity of 0.068 t/ha, followed by MLYP (0.044 t/ha) and HHHP (0.041 t/ha). HHHP has the highest TP discharge intensity, with an annual average intensity of 0.0051 t/ha, followed by SC (0.0038 t/ha) and MLYP (0.0031 t/ha). District-based agro-ecological restoration strategies were accordingly proposed considering FNPS discharge magnitude and intensity concurrently. In most agricultural districts, with the growing economic output in plantation, the FNPS load showed an increase followed by a decrease or to leveling off. Furthermore, with the increasing TN/TP economic partial productivity, the FNPS TN/TP discharge intensities reached the climax, then declined or tended to be flattening out.

Codon usage bias and genetic diversity in chloroplast genomes of Elaeagnus species (Myrtiflorae: Elaeagnaceae).

Codon usage bias (CUB) reveals the characteristics of species and can be utilized to understand their evolutionary relationship, increase the target genes' expression in the heterologous receptor plants, and further provide theoretic assistance for correlative study on molecular biology and genetic breeding. The chief aim of this work was to analyze the CUB in chloroplast (cp.) genes in nine Elaeagnus species to provide references for subsequent studies. The codons of Elaeagnus cp. genes preferred to end with A/T bases rather than with G/C bases. Most of the cp. genes were prone to mutation, while the rps7 genes were identical in sequences. Natural selection was inferred to have a powerful impact on the CUB in Elaeagnus cp. genomes, and their CUB was extremely strong. In addition, the optimal codons were identified in the nine cp. genomes based on the relative synonymous codon usage (RSCU) values, and the optimal codon numbers were between 15 and 19. The clustering analyses based on RSCU were contrasted with the maximum likelihood (ML)-based phylogenetic tree derived from coding sequences, suggesting that the t-distributed Stochastic Neighbor Embedding clustering method was more appropriate for evolutionary relationship analysis than the complete linkage method. Moreover, the ML-based phylogenetic tree based on the conservative matK genes and the whole cp. genomes had visible differences, indicating that the sequences of specific cp. genes were profoundly affected by their surroundings. Following the clustering analysis, Arabidopsis thaliana was considered the optimal heterologous expression receptor plant for the Elaeagnus cp. genes. Supplementary Information: The online version contains supplementary material available at 10.1007/s12298-023-01289-6.

Expansion and Diversification of the 14-3-3 Gene Family in Camellia sinensis.

14-3-3 proteins are signal moderators in sensing various stresses and play essential functions in plant growth and development. Although, 14-3-3 gene families have been identified and characterized in many plant species, its evolution has not been studied systematically. In this study, the plant 14-3-3 family was comprehensively analyzed from green algae to angiosperm. Our result indicated that plant 14-3-3 originated during the early evolutionary history of green algae and expanded in terricolous plants. Twenty-six 14-3-3 genes were identified in the tea genome. RNA-seq analysis showed that tea 14-3-3 genes display different expression patterns in different organs. Moreover, the expression of most tea 14-3-3 genes displayed variable expression patterns under different abiotic and biotic stresses. In conclusion, our results elucidate the evolutionary origin of plant 14-3-3 genes, and beneficial for understanding their biological functions and improving tea agricultural traits in the future.

Population genomics of Zea species identifies selection signatures during maize domestication and adaptation.

BACKGROUND: Maize (Zea mays L. ssp. mays) was domesticated from teosinte (Zea mays ssp. parviglumis) about 9000 years ago in southwestern Mexico and adapted to a range of environments worldwide. Researchers have depicted the maize domestication and adaptation processes over the past two decades, but efforts have been limited either in sample size or genetic diversity. To better understand these processes, we conducted a genome-wide survey of 982 maize inbred lines and 190 teosinte accessions using over 40,000 single-nucleotide polymorphism markers. RESULTS: Population structure, principal component analysis, and phylogenetic trees all confirmed the evolutionary relationship between maize and teosinte, and determined the evolutionary lineage of all species within teosinte. Shared haplotype analysis showed similar levels of ancestral alleles from Zea mays ssp. parviglumis and Zea mays ssp. mexicana in maize. Scans for selection signatures identified 394 domestication sweeps by comparing wild and cultivated maize and 360 adaptation sweeps by comparing tropical and temperate maize. Permutation tests revealed that the public association signals for flowering time were highly enriched in the domestication and adaptation sweeps. Genome-wide association study identified 125 loci significantly associated with flowering-time traits, ten of which identified candidate genes that have undergone selection during maize adaptation. CONCLUSIONS: In this study, we characterized the history of maize domestication and adaptation at the population genomic level and identified hundreds of domestication and adaptation sweeps. This study extends the molecular mechanism of maize domestication and adaptation, and provides resources for basic research and genetic improvement in maize.

A rapid sustainability dynamic assessment of the USA and China 1995-2018.

The USA and China have the world's largest GDP and greenhouse gas emissions, as well as having a significant proportion of the world's total population. These form part of the triumvirate of influencing factors which potentially detrimentally impact upon sustainability. How the triumvirate factors have impacted upon the sustainability of the USA and China has not been adequately determined quantitatively or holistically. Understanding this can strongly indicate how significant levels of population, economic growth, and greenhouse gas emissions impacts and influences the indicated levels, nature, and dynamics of sustainability occurring. Therefore, to contribute towards developing understanding of the issues outlined, the Sustainability Dynamics Framework is applied to conduct a sustainability dynamic assessment for the USA and China 1995-2018. The results indicate both countries are consistently unsustainable due to significant levels of triumvirate-based factors, with differing dynamic trends, characteristics, and potential outcomes. The current static nature of the USA's environment-human system provides the opportunity for a managed co-evolutionary strategy. China's dynamic trends indicate a significant increasingly stressed environment-human system, which is at risk of a potential systemic crisis or collapse unless stringent corrective measures are implemented. In order for both countries to return to sustainability, there needs to be population mitigation to return them to optimum levels, and the use of contractionary policies and rapid decarbonation of their economies and societies. Unless this occurs, then their continued reticence to undertaking the urgent and necessary measures required not only risking their own future, but the world as well.

Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir).

BACKGROUND: R2R3-MYB is a class of transcription factor crucial in regulating secondary cell wall development during wood formation. The regulation of wood formation in gymnosperm has been understudied due to its large genome size. Using Single-Molecule Real-Time sequencing, we obtained full-length transcriptomic libraries from the developmental stem of Cunninghamia lanceolata, a perennial conifer known as Chinese fir. The R2R3-MYB of C. lanceolata (hereafter named as ClMYB) associated with secondary wall development were identified based on phylogenetic analysis, expression studies and functional study on transgenic line. RESULTS: The evolutionary relationship of 52 ClMYBs with those from Arabidopsis thaliana, Eucalyptus grandis, Populus trichocarpa, Oryza sativa, two gymnosperm species, Pinus taeda, and Picea glauca were established by neighbour-joining phylogenetic analysis. A large number of ClMYBs resided in the woody-expanded subgroups that predominated with the members from woody dicots. In contrast, the woody-preferential subgroup strictly carrying the members of woody dicots contained only one candidate. The results suggest that the woody-expanded subgroup emerges before the gymnosperm/angiosperm split, while most of the woody-preferential subgroups are likely lineage-specific to woody dicots. Nine candidates shared the same subgroups with the A. thaliana orthologs, with known function in regulating secondary wall development. Gene expression analysis inferred that ClMYB1/2/3/4/5/26/27/49/51 might participate in secondary wall development, among which ClMYB1/2/5/26/27/49 were significantly upregulated in the highly lignified compression wood region, reinforcing their regulatory role associated with secondary wall development. ClMYB1 was experimentally proven a transcriptional activator that localised in the nucleus. The overexpression of ClMYB1 in Nicotiana benthamiana resulted in an increased lignin deposition in the stems. The members of subgroup S4, ClMYB3/4/5 shared the ERF-associated amphiphilic repression motif with AtMYB4, which is known to repress the metabolism of phenylpropanoid derived compounds. They also carried a core motif specific to gymnosperm lineage, suggesting divergence of the regulatory process compared to the angiosperms. CONCLUSIONS: This work will enrich the collection of full-length gymnosperm-specific R2R3-MYBs related to stem development and contribute to understanding their evolutionary relationship with angiosperm species.

Genetic polymorphism analysis of 23 STR loci in the Tujia population from Chongqing, Southwest China.

To evaluate the applicability of 23 autosomal STR loci (D10S1248, D11S4463, D12ATA63, D14S1434, D17S1301, D18S853, D1GATA113, D1S1627, D6S1017, D20S1082, D20S482, D17S974, D22S1045, D1S1677, D2S1776, D2S441, D3S4529, D4S2408, D9S1122D5S2500, D6S474, D18S51, D9S2157) included in DNA Typer™ 25 Kit for individual identification and parentage testing, allele frequencies and forensic efficiency parameters were first obtained from healthy, unrelated 506 Chongqing Tujia individuals. A total of 1012 alleles were identified in 23 STR loci, and allele frequencies ranged from 0.001 to 0.5761. The combined power of discrimination (CPD) and the combined power of exclusion (CPE) of the 23 STR loci were 0.999999999999999999999753 and 0.99999967, respectively. These results suggested that 23 autosomal STR loci could be used as an effective tool for forensic application in Chongqing Tujia population. Comprehensive comparisons were conducted based on the analysis of genetic distance, principal component analysis (PCA), multidimensional scaling plot (MDS), and phylogenetic tree to explore the interpopulation genetic relationship. Our results revealed that Chongqing Tujia keeps the more relatively genetic similarity with Hunan Han, Hubei Tujia, and Sichuan Han, which could be interpreted by that those populations were originated from the same ethnic ancestor or genetic communication were happened in adjacent areas.

Mitochondrial genome of Theileria uilenbergi endemic in sheep and goats in China.

Mitochondrial genomes provide new insights that help elucidating biological features, genetic evolution, and classification of protozoans. Theileria uilenbergi (T. uilenbergi), transmitted by Haemaphysalis qinghaiensis and H. longicornis, is considered as highly pathogenic to sheep and goats in China. This study reports and outlines features of its mitochondrial genome. The T. uilenbergi mitochondrial genome is a linear monomeric molecule of 6.0 kb length, which encodes three protein-coding genes named cytochrome c oxidase I (cox1), cytochrome b (cob), and cytochrome c oxidase III (cox3), as well as six large subunit (LSU) rRNA gene fragments, and ends in terminal inverted repeats (TIRs). The array structure and organization of the mitochondrial genome of T. uilenbergi is identical to that of T. parva. Phylogenetic analysis based on the amino acid sequences of cox1, cob, and cox3 genes suggests that T. uilenbergi is distantly related to the group of transforming Theileria species such as T. parva. This study contributes to a comprehensive understanding of the phylogeny and evolution of the mitochondrial genome of piroplasms and provides useful information of diagnostic marker for T. uilenbergi.

A new graph-theoretic approach to determine the similarity of genome sequences based on nucleotide triplets.

Methods of finding sequence similarity play a significant role in computational biology. Owing to the rapid increase of genome sequences in public databases, the evolutionary relationship of species becomes more challenging. But traditional alignment-based methods are found inappropriate due to their time-consuming nature. Therefore, it is necessary to find a faster method, which applies to species phylogeny. In this paper, a new graph-theory based alignment-free sequence comparison method is proposed. A complete-bipartite graph is used to represent each genome sequence based on its nucleotide triplets. Subsequently, with the help of the weights of edges of the graph, a vector descriptor is formed. Finally, the phylogenetic tree is drawn using the UPGMA algorithm. In the present case, the datasets for comparison are related to mammals, viruses, and bacteria. In most of the cases, the phylogeny in the present case is found to be more satisfactory as compared to earlier methods.

Genome-Wide Identification and Analysis of the Polycomb Group Family in Medicago truncatula.

Polycomb group (PcG) proteins, which are important epigenetic regulators, play essential roles in the regulatory networks involved in plant growth, development, and environmental stress responses. Currently, as far as we know, no comprehensive and systematic study has been carried out on the PcG family in Medicago truncatula. In the present study, we identified 64 PcG genes with distinct gene structures from the M. truncatula genome. All of the PcG genes were distributed unevenly over eight chromosomes, of which 26 genes underwent gene duplication. The prediction of protein interaction network indicated that 34 M. truncatula PcG proteins exhibited protein-protein interactions, and MtMSI1;4 and MtVRN2 had the largest number of protein-protein interactions. Based on phylogenetic analysis, we divided 375 PcG proteins from 27 species into three groups and nine subgroups. Group I and Group III were composed of five components from the PRC1 complex, and Group II was composed of four components from the PRC2 complex. Additionally, we found that seven PcG proteins in M. truncatula were closely related to the corresponding proteins of Cicer arietinum. Syntenic analysis revealed that PcG proteins had evolved more conservatively in dicots than in monocots. M. truncatula had the most collinearity relationships with Glycine max (36 genes), while collinearity with three monocots was rare (eight genes). The analysis of various types of expression data suggested that PcG genes were involved in the regulation and response process of M. truncatula in multiple developmental stages, in different tissues, and for various environmental stimuli. Meanwhile, many differentially expressed genes (DEGs) were identified in the RNA-seq data, which had potential research value in further studies on gene function verification. These findings provide novel and detailed information on the M. truncatula PcG family, and in the future it would be helpful to carry out related research on the PcG family in other legumes.

Roles of the 2-Oxoglutarate-Dependent Dioxygenase Superfamily in the Flavonoid Pathway: A Review of the Functional Diversity of F3H, FNS I, FLS, and LDOX/ANS.

The 2-oxoglutarate-dependent dioxygenase (2-OGD) superfamily is one of the largest protein families in plants. The main oxidation reactions they catalyze in plants are hydroxylation, desaturation, demethylation, epimerization, and halogenation. Four members of the 2-OGD superfamily, i.e., flavonone 3ß-hydroxylase (F3H), flavones synthase I (FNS I), flavonol synthase (FLS), and anthocyanidin synthase (ANS)/leucoanthocyanidin dioxygenase (LDOX), are present in the flavonoid pathway, catalyzing hydroxylation and desaturation reactions. In this review, we summarize the recent research progress on these proteins, from the discovery of their enzymatic activity, to their functional verification, to the analysis of the response they mediate in plants towards adversity. Substrate diversity analysis indicated that F3H, FNS Ⅰ, ANS/LDOX, and FLS perform their respective dominant functions in the flavonoid pathway, despite the presence of functional redundancy among them. The phylogenetic tree classified two types of FNS Ⅰ, one mainly performing FNS activity, and the other, a new type of FNS present in angiosperms, mainly involved in C-5 hydroxylation of SA. Additionally, a new class of LDOXs is highlighted, which can catalyze the conversion of (+)-catechin to cyanidin, further influencing the starter and extension unit composition of proanthocyanidins (PAs). The systematical description of the functional diversity and evolutionary relationship among these enzymes can facilitate the understanding of their impacts on plant metabolism. On the other hand, it provides molecular genetic evidence of the chemical evolution of flavonoids from lower to higher plants, promoting plant adaptation to harsh environments.