Coding Domain Sequence Research Articles

Expression, purification and functional validation of phage depolymerase from hypervirulent Klebsiella pneumoniae serotype K1

Objective: To express and purify the phage depolymerase from hypervirulent Klebsiella pneumoniae (hvKp) serotype K1 and validate its function. Methods: Phage that infected serotype K1-type hvKp was isolated from hospital sewage. The biology and morphology of the phage were determined by plaque assay and transmission electron microscopy. The whole genome of the phage was sequenced by the Illumina HiSeq 2500 platform. The presence of depolymerase was determined by observing the plaque halo. Bioinformatic analysis and prokaryotic protein expression system were further used to predict and identify phage depolymerase. The depolymerase gene fragment was obtained by PCR and cloned into the pET28a expression vector, and the expression and purification of the depolymerase were completed in strain BL21. The depolymerase activities on the capsular polysaccharide of serotype K1-type hvKp clinical isolates were detected by plaque assay and low-speed centrifugation assay. Results: A lytic phage (phiA2) that infected serotype K1-type hvKp clinical isolate was isolated from hospital sewage. It was typical of the Caudovirales order and Autographiviridae family, and its whole genome was 43 526 bp in length and contained 51 coding domain sequences. The phage phiA2-derived depolymerase phiA2-dep was predicted, expressed and purified. The plaque assay and low-speed centrifugation assay indicated that the depolymerase phiA2-dep had good lytic activity on the capsular polysaccharide of serotype K1-type hvKp clinical isolates. Conclusion: Depolymerase phiA2-dep can specifically degrade the capsular polysaccharide of serotype K1-type hvKp, which has potential application value in treating bacterial infection.

Tidings from the Tides–De novo transcriptome assembly of the endemic estuarine bivalve Villorita cyprinoides

The Indian black clam Villorita cyprinoides Gray, 1825, is an economically valuable estuarine bivalve that faces challenges from multiple stressors and anthropogenic pressures. However, limited genomic resources have hindered molecular investigations into the impact of these stressors on clam populations. Here, we have generated the first transcriptomic reference datasets for V. cyprinoides to address this knowledge gap. A total of 25,040,592 and 22,486,217 million Illumina paired-end reads generated from two individuals were assembled using Trinity and rnaSPAdes. From the 47,607 transcripts identified as Coding Domain Sequences, 37,487 returned positive BLAST hits against six different databases. Additionally, a total of 14,063 Single Sequence Repeats were identified using GMATA. This study significantly enhances the genetic understanding of V. cyprinoides, a potential candidate for aquaculture that supports the livelihoods of many people dependent on small-scale fisheries. The data generated provides insights into broader genealogical connections within the family Cyrenidae through comparative transcriptomics. Furthermore, this transcriptional profile serves as baseline data for future studies in toxicological and conservation genetics.

Brassica napus BnaA09.MYB52 enhances seed coat mucilage accumulation and tolerance to osmotic stress during seed germination in Arabidopsis thaliana.

Seed coat mucilage plays an important role in promoting seed germination under adversity. Previous studies have shown that Arabidopsis thaliana MYB52 (AtMYB52) can positively regulate seed coat mucilage accumulation. However, the role of Brassica napus MYB52 (BnaMYB52) in accumulation of seed coat mucilage and tolerance to osmotic stress during seed germination remains largely unknown. We cloned the BnaA09.MYB52 coding domain sequence from B. napus cv ZS11, identified its conserved protein domains and elucidated its relationship with homologues from a range of plant species. Transgenic plants overexpressing BnaA09.MYB52 in the A. thaliana myb52-1 mutant were generated through Agrobacterium-mediated transformation and used to assess the possible roles of BnaA09.MYB52 in accumulation of seed coat mucilage and tolerance to osmotic stress during seed germination. Subcellular localization and transcriptional activity assays demonstrated that BnaA09.MYB52 functions as a transcription factor. RT-qPCR results indicate that BnaA09.MYB52 is predominantly expressed in roots and developing seeds of B. napus cv ZS11. Introduction of BnaA09.MYB52 into myb52-1 restored thinner seed coat mucilage in this mutant to levels in the wild type. Consistently, expression levels of three key genes participating in mucilage formation in developing seeds of myb52-1 were also restored to wild type levels by overexpressing BnaA09.MYB52. Furthermore, BnaA09.MYB52 was induced by osmotic stress during seed germination in B. napus, and ectopic expression of BnaA09.MYB52 successfully corrected sensitivity of the myb52-1 mutant to osmotic stress during seed germination. These findings enhance our understanding of the functions of BnaA09.MYB52 and provide a novel strategy for future B. napus breeding.

A novel phage carrying capsule depolymerase effectively relieves pneumonia caused by multidrug-resistant Klebsiella aerogenes

BackgroundKlebsiella aerogenes can cause ventilator-associated pneumonia by forming biofilms, and it is frequently associated with multidrug resistance. Phages are good antibiotic alternatives with unique advantages. There has been a lack of phage therapeutic explorations, kinetic studies, and interaction mechanism research targeting K. aerogenes.MethodsPlaque assay, transmission electron microscopy and whole-genome sequencing were used to determine the biology, morphology, and genomic characteristics of the phage. A mouse pneumonia model was constructed by intratracheal/endobronchial delivery of K. aerogenes to assess the therapeutic effect of phage in vivo. Bioinformatics analysis and a prokaryotic protein expression system were used to predict and identify a novel capsule depolymerase. Confocal laser scanning microscopy, Galleria mellonella larvae infection models and other experiments were performed to clarify the function of the capsule depolymerase.ResultsA novel lytic phage (pK4-26) was isolated from hospital sewage. It was typical of the Podoviridae family and exhibited serotype specificity, high lytic activity, and high environmental adaptability. The whole genome is 40,234 bp in length and contains 49 coding domain sequences. Genomic data show that the phage does not carry antibiotic resistance, virulence, or lysogenic genes. The phage effectively lysed K. aerogenes in vivo, reducing mortality and alleviating pneumonia without promoting obvious side effects. A novel phage-derived depolymerase was predicted and proven to be able to digest the capsule, remove biofilms, reduce bacterial virulence, and sensitize the bacteria to serum killing.ConclusionsThe phage pK4-26 is a good antibiotic alternative and can effectively relieve pneumonia caused by multidrug-resistant K. aerogenes. It carries a depolymerase that removes biofilms, reduces virulence, and improves intrinsic immune sensitivity.

Biological characteristics and genomic analysis of a novel Escherichia phage Kayfunavirus CY1.

As the problem of bacterial resistance becomes serious day by day, bacteriophage as a potential antibiotic substitute attracts more and more researchers' interest. In this study, Escherichia phage Kayfunavirus CY1 was isolated from sewage samples of swine farms and identified by biological characteristics and genomic analysis. One-step growth curve showed that the latent period of phage CY1 was about 10min, the outbreak period was about 40min and the burst size was 35 PFU/cell. Analysis of the electron microscopy and whole-genome sequence showed that the phage should be classified as a member of the Autographiviridae family, Studiervirinae subfamily. Genomic analysis of phage CY1 (GenBank accession no. OM937123) revealed a genome size of 39,173bp with an average GC content of 50.51% and 46 coding domain sequences (CDSs). Eight CDSs encoding proteins involved in the replication and regulation of phage DNA, 2 CDSs encoded lysis proteins, 14 CDSs encoded packing and morphogenesis proteins. Genomic and proteomic analysis identified no sequence that encoded for virulence factor, integration-related proteins or antibiotic resistance genes. In summary, morphological and genomics suggest that phage CY1 is more likely a novel Escherichia phage.

High prevalence of fluconazole resistant Candida tropicalis among candiduria samples in China: An ignored matter of concern.

The rapid rise of azole resistance in Candida tropicalis causing invasive infections has become a public health concern; however, the prevalence of resistant isolates in urine samples was not well studied, because the clinical significance of candiduria was not unambiguous due to possible host colonization. We performed a 12-year laboratory-based surveillance study of C. tropicalis causing either invasive infection or candiduria and studied their susceptibility profiles to common antifungal drugs. The complete coding domain sequence of the ERG11 gene was amplified in all fluconazole resistant isolates, and aligned with the wild-type sequence to detect nucleotide mutations. A total of 519 unique C. tropicalis strains isolates, 69.9% of which were isolated from urine samples and remaining 30.1% were invasive strains. Overall, 16.5% isolates were confirmed to be resistant to fluconazole, of which 91.9% were cross-resistant voriconazole. Of note, at the beginning of surveillance (2010-2011), the fluconazole resistance rates were low in both candiduria and invasive groups (6.8% and 5.9%, respectively). However, the resistant rate in the candiduria group significantly increased to 29.5% since 2012-2013 (p = 0.001) and stayed high since then, whilst the resistance rate in the invasive group only showed a gradually increasing trends till 2021 (p > 0.05). Sequence analysis of ERG11 from fluconazole-resistant strains revealed the prevalence of A395T/W mutations were relatively low (16.7%) in the beginning but reached 87.5-100% after 2014. Moreover, the A395W heterozygous mutation isolates became predominant (>60% of resistant strains) after 2016, and indeed isolates carrying corresponding amino acid substitution (Y132F) was highly resistant to fluconazole with MIC50 exceeded 256 μg/ml. Our study revealed high azole resistant rate in candiduria with its increasing trends observed much earlier than stains causing invasive infections. Given antimicrobial resistance as a critical "One Health" issue, the emergence of antifungal resistance in Candida species that are common commensal colonizers in the human body should be concerned.

PiRNAs may regulate expression of candidate genes of esophageal adenocarcinoma.

Elucidation of ways to regulate the expression of candidate cancer genes will contribute to the development of methods for cancer diagnosis and therapy. The aim of the present study was to show the role of piRNAs as efficient regulators of mRNA translation of esophageal adenocarcinoma (EAC) candidate genes. We used bioinformatic methods to study the interaction characteristics of up to 200 thousand piRNAs with mRNAs of 38 candidate EAC genes. The piRNAs capable of binding to mRNA of AR, BTG3, CD55, ERBB3, FKBP5, FOXP1, LEP, SEPP1, SMAD4, and TP53 genes with high free energy by the formation of hydrogen bonds between canonical (G-C, A-U) and noncanonical (G-U, A-C) piRNA and mRNA nucleotide pairs were revealed. The organization of piRNA binding sites (BSs) in the mRNA of candidate genes was found to overlap nucleotide sequences to form clusters. Clusters of piRNA BSs were detected in the 5'-untranslated region, coding domain sequence, and 3'-untranslated region of mRNA. Due to the formation of piRNA binding site clusters, compaction of BSs occurs and competition between piRNAs for binding to mRNA of candidate EAC genes occurs. Associations of piRNA and candidate genes were selected for use as markers for the diagnosis of EAC.

Isolation and characterization of a novel Escherichia coli phage Kayfunavirus ZH4.

Escherichia coli, a gram-negative bacterium, was generally considered conditional pathogenic bacteria and the proportion of bacteria resistant to commonly used specified antibacterial drugs exceeded 50%. Phage therapeutic application has been revitalized since antibiotic resistance in bacteria was increasing. Compared with antibiotics, phage is the virus specific to bacterial hosts. However, further understanding of phage-host interactions is required. In this study, a novel phage specific to a E. coli strain, named as phage Kayfunavirus ZH4, was isolated and characterized. Transmission electron microscopy showed that phage ZH4 belongs to the family Autographiviridae. The whole-genome analysis showed that the length of phage ZH4 genome was 39,496 bp with 49 coding domain sequence (CDS) and no tRNA was detected. Comparative genome and phylogenetic analysis demonstrated that phage ZH4 was highly similar to phages belonging to the genus Kayfunavirus. Moreover, the highest average nucleotide identity (ANI) values of phage ZH4 with all the known phages was 0.86, suggesting that ZH4 was a relatively novel phage. Temperature and pH stability tests showed that phage ZH4 was stable from 4° to 50°C and pH range from 3 to 11. Host range of phage ZH4 showed that there were only 2 out of 17 strains lysed by phage ZH4. Taken together, phage ZH4 was considered as a novel phage with the potential for applications in the food and pharmaceutical industries.

MicroRNA profiling reveals miR-145-5p inhibits goat myoblast differentiation by targeting the coding domain sequence of USP13.

MicroRNAs (miRNAs) are evolutionarily conserved endogenous small non-coding RNAs that play critical roles in skeletal muscle development. In this study, we identified putative miRNAs that were differentially expressed in the longissimus dorsi muscle between fetus (75days of pregnancy) and lamb (1day of age). We detected 1208miRNAs, 313 of which were differentially expressed. In particular, we found that miR-145-5p was differentially and highly expressed in lamb skeletal muscle. In addition, our results demonstrated that overexpression of miR-145-5p inhibited the differentiation and apoptosis of goat primary myoblasts (GPMs), whereas knockdown of miR-145-5p had the opposite effect. The coding domain sequence (CDS) of ubiquitin-specific peptidase 13 (USP13) was predicted and validated as a target of miR-145-5p. We also demonstrated that the influence as a key regulator of GPMs differentiation is primarily mediated by targeting and inhibiting USP13. Taken together, these results revealed a novel pathway in skeletal muscle development in which miR-145-5p targets CDS region of USP13 to regulate differentiation and apoptosis of GPMs.

Isolation and expression of a new gene <italic>HvMEL1 AGO</italic> in Tibetan hulless barley under leaf stripe stress

<p id="C3">To explore <italic>AGO</italic> genes related to leaf stripe in Tibetan hulless barley (BLS), we obtained a differentially expressed <italic>HvMEL1 AGO </italic>gene from the transcriptional sequencing in normal and diseased leaves from resistant Tibetan hulless barley variety Kunlun 14 and the susceptible variety 1141. The <italic>HvMEL1 AGO</italic> gene was 3462 bp in length, of which the CDS (coding domain sequence) was 100% consistent in Kunlun 14 and 1141 varieties. The length of <italic>HvMEL1 AGO</italic> gene was 3161 bp without intron, contained a 3129 bp open reading frame, encoded 1043 amino acids, had a theoretical isoelectric point of 9.33, and has a predicted molecular weight of 115,865.58 Da. Protein sequence analysis showed that HvMEL1 AGO was a hydrophilic unstable acidic protein with highly conserved structural domains of DUF1785, PAZ and PIWI, belonging to <italic>AGO</italic> gene family. Phylogenetic tree analysis showed that HvMEL1 AGO belonged to the AGO1 class in <italic>Arabidopsis</italic> AGO family phylogenetic tree with HvAGO12, HvAGO18, HvAGO1D and HvAGO1B in barley AGO family, and were closely related to HvAGO12. The predicted protein interactions showed that the known proteins that acted closely with MEL1 in rice were the DCL classes, DCL1, DCL2A, DCL3A, DCL3B and DCL4, respectively. Semi-quantitative and quantitative PCR results indicated that the relative expression levels of <italic>HvMEL1 AGO</italic> gene were significantly down-regulated in 1141 and Kunlun 14 under BLS. We hypothesized that <italic>HvMEL1 AGO</italic> gene played an important role in the streak resistance mechanism of barley. This study lays the foundation for exploring the role and regulatory mechanism of <italic>HvMEL1 AGO</italic> in the process of BLS resistance.

Cloning of FITM2 gene and investigating its expression levels in Banna miniature inbred pig (Sus scrofa) tissues

Fat storage-inducing transmembrane protein 2 (FITM2) plays an important role in regulating lipid storage and could be regarded as a candidate gene for intramuscular fat deposition in pigs. The aim of this study was to clone the coding domain sequence (CDS) of FITM2 gene, to compare the nucleotide acid and deduced amino acid sequences between breeds and species, to analyze the structure and characteristics of protein and to detect the expression profile of gene. The results exhibited that the CDS of FITM2 gene was 789 bp in length. The mutation of nucleotide acids led to the mutation of deduced amino acids between Banna miniature inbred pigs and other two breeds (Yorkshire × Landrace pigs and Duroc × (Landrace × Yorkshire) pigs). It was indicated that high identities of nucleotide acid and deduced amino acid sequences between Banna miniature inbred pigs and other species. The deduced amino acids were composed of loops and alpha helices in the structure. FITM2 protein may be a 30 kDa hydrophobic protein with 26 phosphorylation sites and one potential N-glycosylated site. FITM2 gene was widely expressed in various tissues, and the highest expression level was in adipose tissue.

Sequence polymorphisms of PR39 cathelicidins and extensive copy variations in commercial pig breeds

Copy number polymorphisms (CNPs) of antimicrobial peptides (AMPs) in livestock can influence the innate immune response of individuals. We conducted a high-resolution analysis of the genomic variations of porcine cathelicidin PR39 using cloned PR39 amplicons corresponding to the 5′ untranslated region (UTR) to 3′ UTR from four individuals of three different pig breeds. We identified 15 different sequences corresponding to 9 different coding domain sequences (CDSs), encoding 7 different protein sequences consisting of 3 functional and 4 non-functional forms. Subsequently, we developed a PR39 CNP typing method using real-time polymerase chain reaction (PCR) and analyzed the PR39 copy numbers from 44 pigs of six breeds. Significant variations in PR39 copies ranging from 2 to 10 copies, with a mean copy number of 5, were observed among all commercial breeds, except the wild boar. Among the different breeds, the PR39 copy number was highest (10) in Korean native pigs. Gene expression analysis showed that PR39 expression was correlated with the copy number. Moreover, the comparative analysis of the cathelicidin cluster-containing region among eight mammalian species showed the complete evolutionary conservation of the region, except for differences in the degree of cathelicidin expansion in each species. Therefore, characterization of CNPs in AMP genes could aid in improving the genetic potential of innate immune responses in livestock animals.

Genome-wide discovery of G-quadruplexes in barley

G-quadruplexes (G4s) are four-stranded nucleic acid structures with closely spaced guanine bases forming square planar G-quartets. Aberrant formation of G4 structures has been associated with genomic instability. However, most plant species are lacking comprehensive studies of G4 motifs. In this study, genome-wide identification of G4 motifs in barley was performed, followed by a comparison of genomic distribution and molecular functions to other monocot species, such as wheat, maize, and rice. Similar to the reports on human and some plants like wheat, G4 motifs peaked around the 5′ untranslated region (5′ UTR), the first coding domain sequence, and the first intron start sites on antisense strands. Our comparative analyses in human, Arabidopsis, maize, rice, and sorghum demonstrated that the peak points could be erroneously merged into a single peak when large window sizes are used. We also showed that the G4 distributions around genic regions are relatively similar in the species studied, except in the case of Arabidopsis. G4 containing genes in monocots showed conserved molecular functions for transcription initiation and hydrolase activity. Additionally, we provided examples of imperfect G4 motifs.

Global alternative splicing landscape of skeletal muscle atrophy induced by hindlimb unloading.

BackgroundLong-term exposure to microgravity will cause skeletal muscle atrophy, which can cause serious harm to astronauts in space travel. Therefore, it is important to explore skeletal muscle atrophy’s molecular mechanism for its prevention and treatment. However, as an important regulatory approach of skeletal muscle physiology, the role of alternative splicing in skeletal muscle atrophy, especially skeletal muscle atrophy caused by disuse, is unclear.MethodsWe established a rat hindlimb unloading model and performed RNA sequencing on soleus muscle, which was seriously atrophied during unloading. Several bioinformatics methods were used to identify alternative splicing events and determine their gene functions.ResultsMany alternative splicing events were found to occur at different time points (12 h, 24 h, 36 h, 3 days, and 7 days) after hindlimb unloading. These differential alternative splicing events mainly occurred in the gene's coding domain sequence region, and 59% of the alternative splicing events caused open reading frameshift. Bioinformatics analysis results showed that genes with different alternative splicing events were enriched in multiple pathways related to muscle atrophy, including the insulin signaling pathway, endocytosis, mitophagy, and ubiquitin-proteasome pathway. Moreover, alternative splicing of several deubiquitinase genes persisted during skeletal muscle atrophy induced by unloading. Additionally, we identified 10 differentially expressed RNA binding proteins during skeletal muscle atrophy induced by unloading, mainly containing Xpo4, Eif4e2, P4ha1, Lrrfip1, Zc3h14, Emg1, Hnrnp h1, Mbnl2, RBfox1, and Mbnl1. Hnrnp h1 and Mbnl2 were significantly downregulated, and RBfox1 and Mbnl1 were significantly upregulated during skeletal muscle atrophy caused by unloading.ConclusionsTo the best of our knowledge, the present study is the first to propose alternative splicing alterations related to disuse-induced muscle atrophy, emphasizing that alternative splicing is a new focus of attention in the occurrence of muscle atrophy.

Characterization and expression analysis of ascorbate peroxidase from the moss Dicranum scoparium during abiotic stresses

Ascorbate peroxidase (APX) is an important antioxidant enzyme responsible for the conversion of H2O2 to H2O and O2. In this study, APX was studied in the widespread boreal cushion moss Dicranum scoparium. Native PAGE of crude extracts of moss thalli revealed the presence of APX isoforms in D. scoparium with a range of molecular masses. An APX complementary DNA (cDNA) gene of 771-bp length was cloned and designated as DsAPX. The cloned coding domain sequence (CDS) encoded a 256 amino acid polypeptide, and the predicted protein product was calculated to have a molecular mass of 28.4 kDa with an isoelectric point of 5.7. Several highly conserved sites important for enzyme activity were predicted using bioinformatic tools. The DsAPX protein has similarities of 91%, 78%, 67% and 66% with APX homologs from Grimmia pilifera, Physcomitrella patens, Zea mays and Nicotiana tomentosiformis, respectively. The high homology with cytosolic APX from G. pilifera strongly suggests that the cloned DsAPX gene encodes a cytosolic APX. We studied the role of APX in the tolerance of this moss to abiotic stresses. Changes in both APX activity and DsAPX gene expression, estimated using reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR), showed that the gene was up-regulated in response to desiccation/rehydration, and expression was maintained during a heat stress of +50°C. By contrast, DsAPX was down-regulated by a freezing (–20°C) treatment. Results obtained extend our knowledge of the diversity of APX isoforms that occur in bryophytes and suggest that APX may be involved in the tolerance of D. scoparium to abiotic stresses.

Transcriptome profiling of venom gland from wasp species: de novo assembly, functional annotation, and discovery of molecular markers

BackgroundVespa velutina, one of the most aggressive and fearful wasps in China, can cause grievous allergies and toxic reactions, leading to organ failure and even death. However, there is little evidence on molecular data regarding wasps. Therefore, we aimed to provide an insight into the transcripts expressed in the venom gland of wasps.ResultsIn our study, high-throughput RNA sequencing was performed using the venom glands of four wasp species. First, the mitochondrial cytochrome C oxidase submit I (COI) barcoding and the neighbor joining (NJ) tree were used to validate the unique identity and lineage of each individual species. After sequencing, a total of 127,630 contigs were generated and 98,716 coding domain sequences (CDS) were predicted from the four species. The Gene ontology (GO) enrichment analysis of unigenes revealed their functional role in important biological processes (BP), molecular functions (MF) and cellular components (CC). In addition, c-type, p1 type, p2 type and p3 type were the most commonly found simple sequence repeat (SSR) types in the four species of wasp transcriptome. There were differences in the distribution of SSRs and single nucleotide polymorphisms (SNPs) among the four wasp species.ConclusionsThe transcriptome data generated in this study will improve our understanding on bioactive proteins and venom-related genes in wasp venom gland and provide a basis for pests control and other applications. To our knowledge, this is the first study on the identification of large-scale genomic data and the discovery of microsatellite markers from V. tropica ducalis and V. analis fabricius.

Characterization of Novel Broad-Host-Range Bacteriophage DLP3 Specific to Stenotrophomonas maltophilia as a Potential Therapeutic Agent.

A novel Siphoviridae phage specific to the bacterial species Stenotrophomonas maltophilia was isolated from a pristine soil sample and characterized as a second member of the newly established Delepquintavirus genus. Phage DLP3 possesses one of the broadest host ranges of any S. maltophilia phage yet characterized, infecting 22 of 29 S. maltophilia strains. DLP3 has a genome size of 96,852 bp and a G+C content of 58.4%, which is significantly lower than S. maltophilia host strain D1571 (G+C content of 66.9%). The DLP3 genome encodes 153 coding domain sequences covering 95% of the genome, including five tRNA genes with different specificities. The DLP3 lysogen exhibits a growth rate increase during the exponential phase of growth as compared to the wild type strain. DLP3 also encodes a functional erythromycin resistance protein, causing lysogenic conversion of the host D1571 strain. Although a temperate phage, DLP3 demonstrates excellent therapeutic potential because it exhibits a broad host range, infects host cells through the S. maltophilia type IV pilus, and exhibits lytic activity in vivo. Undesirable traits, such as its temperate lifecycle, can be eliminated using genetic techniques to produce a modified phage useful in the treatment of S. maltophilia bacterial infections.

Genome-Wide Discovery of G-Quadruplexes in Wheat: Distribution and Putative Functional Roles.

G-quadruplexes are nucleic acid secondary structures formed by a stack of square planar G-quartets. G-quadruplexes were implicated in many biological functions including telomere maintenance, replication, transcription, and translation, in many species including humans and plants. For wheat, however, though it is one of the world’s most important staple food, no G-quadruplex studies have been reported to date. Here, we computationally identify putative G4 structures (G4s) in wheat genome for the first time and compare its distribution across the genome against five other genomes (human, maize, Arabidopsis, rice, and sorghum). We identified close to 1 million G4 motifs with a density of 76 G4s/Mb across the whole genome and 93 G4s/Mb over genic regions. Remarkably, G4s were enriched around three regions, two located on the antisense and one on the sense strand at the following positions: 1) the transcription start site (TSS) (antisense), 2) the first coding domain sequence (CDS) (antisense), and 3) the start codon (sense). Functional enrichment analysis revealed that the gene models containing G4 motifs within these peaks were associated with specific gene ontology (GO) terms, such as developmental process, localization, and cellular component organization or biogenesis. We investigated genes encoding MADS-box transcription factors and showed examples of G4 motifs within critical regulatory regions in the VRN-1 genes in wheat. Furthermore, comparison with other plants showed that monocots share a similar distribution of G4s, but Arabidopsis shows a unique G4 distribution. Our study shows for the first time the prevalence and possible functional roles of G4s in wheat.

Genome-wide DNA polymorphisms in four Actinidia arguta genotypes based on whole-genome re-sequencing.

Among the genus Actinidia, Actinidia arguta possesses the strongest cold resistance and produces fresh fruit with an intense flavor. To investigate genomic variation that may contribute to variation in phenotypic traits, we performed whole-genome re-sequencing of four A. arguta genotypes originating from different regions in China and identified the polymorphisms using InDel markers. In total, 4,710,650, 4,787,750, 4,646,026, and 4,590,616 SNPs and 1,481,002, 1,534,198, 1,471,304, and 1,425,393 InDels were detected in the ‘Ruby-3’, ‘Yongfeng male’, ‘Kuilv male’, and ‘Hongbei male’ genomes, respectively, compared with the reference genome sequence of cv ‘Hongyang’. A subset of 120 InDels were selected for re-sequencing validation. Additionally, genes related to non-synonymous SNPs and InDels in coding domain sequences were screened for functional analysis. The analysis of GO and KEGG showed that genes involved in cellular responses to water deprivation, sucrose transport, decreased oxygen levels and plant hormone signal transduction were significantly enriched in A. arguta. The results of this study provide insight into the genomic variation of kiwifruit and can inform future research on molecular breeding to improve cold resistance in kiwifruit.

Identification of QTL TGW12 responsible for grain weight in rice based on recombinant inbred line population crossed by wild rice (Oryza minuta) introgression line K1561 and indica rice G1025

BackgroundLimited genetic resource in the cultivated rice may hinder further yield improvement. Some valuable genes that contribute to rice yield may be lost or lacked in the cultivated rice. Identification of the quantitative trait locus (QTL) for yield-related traits such as thousand-grain weight (TGW) from wild rice speices is desired for rice yield improvement.ResultsIn this study, sixteen TGW QTL were identified from a recombinant inbred line (RIL) population derived from the cross between the introgression line K1561 of Oryza minuta and the rice cultivar G1025. TGW12, One of most effective QTL was mapped to the region of 204.12 kb between the marker 2,768,345 and marker 2,853,491 of the specific locus amplified fragment (SLAF). The origin of TGW12 was tested using three markers nearby or within the TGW12 region, but not clarified yet. Our data indicated thirty-two open reading fragments (ORFs) were present in the region. RT-PCR analysis and sequence alignment showed that the coding domain sequences of ORF12, one MADS-box gene, in G1025 and K1561 were different due to alternative slicing, which caused premature transcription termination. The MADS-box gene was considered as a candidate of TGW12.ConclusionThe effective QTL, TGW12, was mapped to a segment of 204.12 kb using RILs population and a MADS-box gene was identified among several candidate genes in the segment. The region of TGW12 should be further narrowed and creation of transgenic lines will reveal the gene function. TGW12 could be applied for improvement of TGW in breeding program.