Transposon Elements Research Articles

Pangenome and pantranscriptome as the new reference for gene family characterisation – a case study of basic helix-loop-helix (bHLH) genes in barley

Genome-wide identification and comparative gene family analyses have been commonly performed to investigate species-specific evolution linked to various traits and molecular pathways. However, most previous studies were limited to gene screening in a single reference genome, failing to account for the gene presence/absence variations (gPAVs) in a species. Here, we propose an innovative pangenome-based approach of gene family analyses based on orthologous gene groups (OGGs). Using the basic helix-loop-helix (bHLH) transcription factor family in barley as an example, we identified 161 ∼ 176 bHLHs in 20 barley genomes, which could be classified into 201 OGGs. These 201 OGGs were further classified into 140 core, 12 soft-core, 29 shell, and 20 line-specific/cloud bHLHs, revealing a complete profile of bHLH in barley. Using a genome-scan approach, we overcome the genome annotation bias and identified on average 1.5 un-annotated core bHLHs per barley genome. We found that all core bHLHs belong to whole genome/segmental duplicates whilst dispensable bHLHs were more likely to result from small scale duplication events. Interestingly, we noticed that the dispensable bHLHs tended to enrich in specific subfamilies SF13, SF27, and SF28, implying the potential biased expansion of specific bHLHs in barley. We found that 50% of the bHLHs contain at least one intact transposon element within the 2kb upstream-to-downstream region. bHLHs with CNV have 1.48 TEs on average, significantly higher than 1.36 for core bHLH without CNV, supporting TEs’ potential role in bHLH expansion. Selection pressure analyses showed that dispensable bHLHs had experienced clear relaxed selection compared to core bHLHs, consistent with their conservation patterns. We further integrate pangenome with recently available barley pantranscriptome data in 5 tissues and discovered apparent transcriptional divergence within and across bHLH subfamilies. We conclude that pangenome-based gene family analyses can better describe the genuine evolution status of bHLHs untapped before and provided novel insights into bHLH evolution in barley. We expect this study will inspire similar analyses in many other gene families and species.

C G composition in transposon-derived genes is increased in FXD with perturbed immune system.

Increasing incidence of Fragile X disorders (FXD) and of immune-mediated disorders in FXD suggests that additional factors besides FMR1 mutations contribute to the pathogenesis. Here, we discovered that the expression levels or splicing of specific transposon element (TE)-derived genes, regulating purine metabolism and immune responses against viral infections are altered in FXD. These genes include HLA genes clustered in chr6p21.3 and viral responsive genes in chr5q15. Remarkably, these TE-derived genes contain a low A T/C G suggesting base substitutions of A T to C G. The TE-derived genes with changed expression levels contained a higher content of 5'-CG-3' dinucleotides in FXD compared to healthy donors. This resembles the genomes of some RNA viruses, which maintain high contents of CG dinucleotides to sustain their latent infection exploiting antiviral responses. Thus, past viral infections may have persisted as TEs, provoking immune-mediated disorders in FXD.

Identification and molecular detection of the pathogen of Phalaenopsis leaf yellowing through genome analysis.

Moth orchids (Phalaenopsis spp.) are globally popular ornamental flowers. However, effective management strategies for Phalaenopsis leaf yellowing remain elusive, making the disease a challenging obstacle affecting moth orchids at various growth stages. This disease manifests as collar rot, leaf yellowing, leaf abscission, and eventually, plant death. The lack of effective management strategies is likely attributed to a limited understanding of the disease pathogenesis and pathogen dissemination pathways. Fusarium phalaenopsidis sp. nov. was established in this study to stabilize the classification status of Phalaenopsis leaf yellowing pathogens using molecular and morphological features. The genome of the holotype strain was sequenced and assembled, revealing its genome structures. Analyses of virulence-related elements, including transposon elements, secondary metabolite biosynthetic gene clusters, effectors, and secreted carbohydrate-active enzymes, shed light on the potential roles of three fast core chromosomes in virulence. Two species-specific primers were designed based on unique gene sequences of two virulence-related proteins through comparative genomics and BLAST screening. The specificity of these primers was validated using isolates of F. phalaenopsidis, non-target species in the Fusarium solani species complex, other Fusarium species complexes, and saprophytic fungi. These results are intended to accelerate the identification of the pathogens, facilitate the study of disease pathogenesis, and pave the way for elucidating pathogen dissemination pathways. Ultimately, they aim to contribute to the formulation of effective control strategies against Phalaenopsis leaf yellowing.

Genome-wide analysis of DNA methylation and transcriptional changes associated with overwintering memory in Brassica rapa L. grown in the field

BackgroundWinter rapeseed, the sole overwintering oilseed crop in northern China, emphasizes winter resilience, yet epigenetic regulatory mechanisms governing overwintering memory remain poorly understood.ResultsIn this study, the root collar tissues from the robust cold-resistant variety Longyou-7 were sampled during the pre-winter period (S1), overwintering periods (S2–S5), and re-greening period (S6), to analyze overall genomic DNA methylation levels using high-performance liquid chromatography (HPLC). The result showed that DNA methylation level exceeded 80% in the S1 stage. Throughout the overwintering periods, methylation levels displayed a decreasing trend in S3, followed by an increase in S5, and a pronounced decrease in S6. Consequently, S1, S3, S5, and S6 periods were chosen for whole-genome bisulfite sequencing analyses to elucidate the overwintering memory mechanisms of Longyou-7. The result revealed that DNA methylation primarily occurs in the CG context in Longyou-7. However, methylation of mC sites is most prevalent in the CHH type, gradually decreasing during overwintering periods. Analysis of methylation patterns in specific genomic regions of Longyou-7 showed that the highest methylation levels in the intergenic region. Moreover, mC sites in repeats and transposon elements are distributed differently across the three contexts. Subsequently, differentially methylated regions and promoters of Longyou-7 were identified during various periods compared to the S1 stage, followed by joint analysis with transcriptome sequencing. Functional enrichment analysis highlighted the involvement of most overlapping genes in the MAPK signaling pathway, plant hormone signal transduction, and starch and sucrose metabolism pathways. Changes in candidate gene expression within these three pathways correlated closely with DNA methylation levels.ConclusionsOur findings underscored the critical role of DNA methylation in regulating the expression of overwintering memory genes in winter rapeseed. These results offer a comprehensive insights into the epigenetic regulatory mechanisms governing winter rapeseed's overwintering memory, while identified overwintering memory genes served as crucial genetic resources for multifaceted breeding of winter-resistant varieties.Graphical

Upregulation of PECTATE LYASE5 by a NAC transcription factor promotes fruit softening in apple.

Flesh firmness is a critical breeding trait that determines consumer selection, shelf life, and transportation. The genetic basis controlling firmness in apple (Malus×domestica Borkh.) remains to be fully elucidated. We aimed to decipher genetic variance for firmness at harvest and develop potential molecular markers for marker-assisted breeding. Maturity firmness for 439 F1 hybrids from a cross of 'Cripps Pink' and 'Fuji' was determined in 2016 and 2017. The phenotype segregated extensively, with a Gaussian distribution. In a combined bulked segregant analysis (BSA) and RNA-sequencing analysis, eighty-four differentially expressed genes were screened from the 10 QTL regions. Interestingly, next-generation re-sequencing analysis revealed a Harbinger-like transposon element insertion upstream of the candidate gene PECTATE LYASE5 (MdPL5); the genotype was associated with flesh firmness at harvest. The presence of this transposon repressed MdPL5 expression and was closely linked to the extra-hard phenotype. MdPL5 was demonstrated to promote softening in apples and tomatoes. Subsequently, using the MdPL5 promoter as bait, MdNAC1-L was identified as a transcription activator that positively regulates ripening and softening in the developing fruit. We also demonstrated that MdNAC1-L could induce the up-regulation of MdPL5, MdPG1, and the ethylene-related genes MdACS1 and MdACO1. Our findings provide insight into TE-related genetic variation and the PL-mediated regulatory network for the firmness of apple fruit.

A high-quality chromosome-scale genome assembly of blood orange, an important pigmented sweet orange variety

Blood orange (BO) is a rare red-fleshed sweet orange (SWO) with a high anthocyanin content and is associated with numerous health-related benefits. Here, we reported a high-quality chromosome-scale genome assembly for Neixiu (NX) BO, reaching 336.63 Mb in length with contig and scaffold N50 values of 30.6 Mb. Furthermore, 96% of the assembled sequences were successfully anchored to 9 pseudo-chromosomes. The genome assembly also revealed the presence of 37.87% transposon elements and 7.64% tandem repeats, and the annotation of 30,395 protein-coding genes. A high level of genome synteny was observed between BO and SWO, further supporting their genetic similarity. The speciation event that gave rise to the Citrus species predated the duplication event found within them. The genome-wide variation between NX and SWO was also compared. This first high-quality BO genome will serve as a fundamental basis for future studies on functional genomics and genome evolution.

First Detection of International High-Risk blaKPC-2-Harbouring Escherichia coli Pandemic Lineage ST648 in Pet Food Packages

The continued worldwide increase in pet ownership has significantly boosted the growth of the pet food industry accompanied by new food safety risks and challenges. This study was designed to determine the occurrence and molecularly characterize multidrug-resistant (MDR) Enterobacterales in pet food. Eighty-six (86) packages of dry and wet pet food purchased in different retail stores were screened for carbapenem-resistant Enterobacterales (CRE). Antimicrobial susceptibility testing was performed by agar dilution technique using EUCAST/BrCAST recommendations. Blue-Carba test was further used to screen for carbapenemase-producing isolates. Isolated CRE strains were identified at the species level using matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry (MALDI-TOF MS). Detection of carbapenemase-encoding genes was carried out by PCR, Sanger sequencing, and whole genome sequencing (WGS). A total of 15 (17.4%) MDR-CRE (Escherichia coli (n = 2), Enterobacter cloacae (n = 10), Leclercia adecarboxylata (n = 2), and Cronobacter spp. (n = 1)) were isolated from 86 pet food samples. In addition to being resistant to beta-lactams, the Gram-negative bacterial isolates were also resistant to aminoglycosides, fluoroquinolones, and tigecycline. Interestingly, two carbapenem-resistant E. coli isolates harboured blaKPC-2 gene. WGS analysis of the two blaKPC-2-producing E. coli isolates revealed that they both belong to ST648 and serotype O153:H2 group. The genetic context of the blaKPC-2 showed that they were carried by an IncN plasmid on a Tn4401b transposon element. To the best of our knowledge, this is the first description of blaKPC-2-harbouring E. coli ST648 pathogens in pet food. The detection of blaKPC-2-harbouring E. coli ST648 pandemic high-risk lineage in pet food is worrisome and a serious “One Health” issue. Therefore, pet food should be considered as a potential vehicle for the transmission of MDR pathogens to companion animals, and a risk factor for the dissemination of these bacterial pathogens to pet animals and their human guardians.

Chromosome-level genome assembly of chub mackerel (Scomber japonicus) from the Indo-Pacific Ocean

Chub mackerels (Scomber japonicus) are a migratory marine fish widely distributed in the Indo-Pacific Ocean. They are globally consumed for their high Omega-3 content, but their population is declining due to global warming. Here, we generated the first chromosome-level genome assembly of chub mackerel (fScoJap1) using the Vertebrate Genomes Project assembly pipeline with PacBio HiFi genomic sequencing and Arima Hi-C chromosome contact data. The final assembly is 828.68 Mb with 24 chromosomes, nearly all containing telomeric repeats at their ends. We annotated 31,656 genes and discovered that approximately 2.19% of the genome contained DNA transposon elements repressed within duplicated genes. Analyzing 5-methylcytosine (5mC) modifications using HiFi reads, we observed open/close chromatin patterns at gene promoters, including the FADS2 gene involved in Omega-3 production. This chromosome-level reference genome provides unprecedented opportunities for advancing our knowledge of chub mackerels in biology, industry, and conservation.

Trinary Fusion with Ligand Binding Domain Castration Is Crucial in Rarg-Driven APL and Render Unresponsive to ATRA Via Allosteric Disability Mechanism

Introduction All-trans retinoic acid (ATRA) has been used with great success in acute promyelocytic leukemia (APL) cases with PML::RARA fusion gene (FG). There are still PML::RARA negative cases that manifest as APL, and they often carry FGs involved in RARA, RARB, or RARG. X::RARG positive APL ( RARG-APL) cases have been increasingly reported in recent years, thanks to the utilization of transcriptome sequencing (RNA-seq) and whole genome sequencing (WGS). Almost all evaluable RARG-APL cases were resistant to ATRA, and no mechanism has been elucidated. Intriguingly, cells transformed by artificial X::RARG were extremely sensitive to ATRA (PMID: 25510432), indicating intricate mechanisms unrevealed. This study first and systematically unveils the distinctive features and pivotal molecular mechanism of RARG-driven APL. Methods A total of 22 RARG-APL cases with RNA-seq and WGS data from 11 centers were enrolled (Fig. 1A). We used Arriba for routine FG calling and manually investigated the non-coding fusion sequence. Reverse transcription PCR (RT-PCR) and Sanger sequencing were used to validate the FGs and to determine whether the RARG 5′ and 3′ fusion events were located on the same cistron. The responsiveness of the fusion proteins to ATRA was evaluated using an optimized UAS/GAL4 reporter system. Results Among the cases enrolled, there were 15 males and 7 females, aged 0.9-69 (median 38) years. Morphological and immunophenotype analyses all showed features of APL. Routine FGs calling reported RARG fusion to a 5′ partner in each case, with variable RARG splicing sites (Fig. 1A). There were 7 RARG 5' partner genes, with CPSF6, HNRNPC, and NUP98 observed in multiple cases. In the 3 index cases, we identified that they all had RARG 3‘ fusion events besides the 5‘ one (Fig. 1A). RT-PCR confirmed that the 5′ and 3′ RARG fusion events were in the same cistron in each case. Such as, the fusion transcript is the insertion of RARG-e4_9 (e refers to exon) between NPM1-e4 and e11 in case #G1, rather than two separate transcripts of NPM1-e4:: RARG-e4 and RARG-9:: NPM1-e11. We named this novel form of tandem splicingtrinary fusion. Analysis of the 19 validation cases further confirmed that all RARG-APL cases had RARG 3′ fusion events (Fig. 1A). We also validated that the 5′ and 3′ RARG fusion events were in the same cistron in 3 more cases, of which archived cDNA were available. Remarkably, the 3' splicing sits of all cases were consistently at the end of RARG-e9, leading to RARG-e10 truncation. The RARG 3' partner gene was the same as its 5' partner in 11 cases. But in the other 11 cases, the RARG 3' partner was a transposon element (TE) sequence, the most common being a LINE-L2a (8 cases). Sequences analysis indicated that the committed locus of the involved LINE-L2a was at 4.2 Kbp upstream of the RARG gene. Gene expression analysis confirmed the aberrant activation of the involved TEs. All TEs involved in the RARG 3' fusions confer a poly_A signal sequence to the fusion transcript, which was essential for a mRNA. In-depth analysis indicated that TEs participate in the formation of RARG-FGs through a transposition mechanism rather than a translocation mechanism. RARG-e10 encodes helix 11_12 (H11_12) of its ligand binding domain (LBD), which plays a pivotal role in allosteric response to ATRA (PMID: 31178221). The trinary fusion protein with RARG LBD-H11_12 castration will lose responsiveness to ATRA via an allosteric disability mechanism (Fig. 1B, 1C). Experiments confirmed that cells transfected with X::RARG fusion with intact LBD responded well to ATRA, while fusion with LBD-H11_12 castration showed no response to physiological and pharmacological concentrations of ATRA (Fig. 1D, 1E). Conclusions This study clearly revealed that trinary fusion and LBD castration are crucial molecular etiologies in RARG-APL. The trinary fusion orchestrates the aberrant activation of RARG (5‘ splicing) and the complete allosteric disability in responding to ATRA (3‘ splicing). Thus, this study elegantly explains the mechanism of leukemogenesis and the extensive ATRA resistance of RARG-APL. The formation of trinary fusion requires two steps of genome splicing or transposition events, which explains the rarity of RARG-APL. This also clarifies that RARG-APL has a significant disparity in molecular mechanism from PML::RARA-positive APL, and we suggest that it should be considered as a separate entity of acute myeloid leukemia.

Endogenous Retroviral Elements Constitute a Promising Vaccine Antigen Source for Haematological Malignancies

Introduction Immunotherapy has revolutionized modern cancer treatment, both within the fields of solid and haematology malignancies. The recent success of checkpoint inhibitor (CPI) therapies (e.g., aPD-1, aCTLA-4) in so called “high” tumour mutational burden (TMB) cancers, such as melanoma and lung cancer (Alexandrov et al. Nature 2013), have found only limited transferability to the haematological setting outside classical Hodgkins lymphoma (Pianko et al. Cancer J. 2018). Recently, the effect of CPIs in high TMB cancers has been further augmented by the combination of CPIs with personalized cancer vaccines (PCVs) in patients with resectable melanoma (NCT03897881). However, recent progress in the identification of novel alternative antigen targets, besides classical missense mutations targeted by PCVs, provides hope for the targeting of haematological cancers with low TMB by immunotherapy-based treatments such as PCVs or even precision-based cancer vaccines. Such sources include endogenous retroviral elements (ERVs) (Saini et al. Nat Commun. 2020), structural rearrangements (Robbe et al. Nat Gen 2022), or activation of transposon elements (Wolf et al. Nat Rev Nephrol. 2023). To investigate the applicability of such approaches in a haematological setting, prospective datasets curated from the database of Genotypes and Phenotypes (dbGaP) and in-house developed antigen prediction platforms, were employed. Methods Data was downloaded from dbGAP and missense mutations (SNVs, INDELs) were determined either from datasets (firehose, TCGA), when available, or through GATK-Mutect2. ERVs were identified using an in-house developed platform (ObsERV), as previously described (Garde et al. BioRxiv 2023). Simulated PCVs were designed for seven selected patients covering high/low TMB and ERV burdens from datasets pertaining to Acute Myeloid Leukaemia (AML), Acute Lymphoblastic Leukaemia (ALL), Multiple Myeloma (MM), Diffuse Large B Cell Lymphoma (DLBCL) and Chronic Lymphocytic Leukaemia (CLL) using in-house platforms PIONEER (neo-epitope identification and ranking for PCV design) and/or ObsERV. PIONEER ranking was based on MHC class I and II ligand propensity, epitope expression level, clonality and tumour specificity Cut-off for inclusion of “efficacious” neoepitopes were based on PIONEER ranking from a recent Phase 1 clinical trial (Mørk et al. ASCO 2023-9551). Due to the high degree of ERV conservation amongst patients within each haematological indications, a precision-based approach was attempted relying on MHC I/II ligand propensity prediction and epitope combination for optimal patient coverage (HLA haplotype) - a similar machine learning approach was proven useful for viral vaccine design (Persson et al Front. Immunol. 2023). Shared vaccines were designed based on HLA frequencies in the human population, which was then applied to the dataset of each patient. The resulting patient specific coverage was then used to compute the total coverage for a given indication. Results and Discussion We found a high presence of ERVs in several haematological malignancies suggesting the potential application of PCVs in these indications, relying on this antigen source (Fig 1A). The combination of missense mutations and ERV sources allowed PCVs to be designed for six of seven patients on average for AML, CLL, and MM compared to only three of seven on average with missense mutations alone. The high level of ERV expression in AML, CLL, and ALL allowed for the generation of vaccines based on shared antigens within each indication. The resulting vaccines had an average coverage, with at least three epitopes, for 84% of AML, 45% of ALL and 35% of CLL patients. In comparison, a missense mutation approach, for the highest TMB cancer (melanoma) would result in only 2.6 % patient coverage with three epitopes. We believe that these findings suggest a novel approach to the treatment of haematological malignancies, either through a PCVs approach, combining multiple sources of antigens, or a precision-based approach with “warehouse” HLA haplotype-specific vaccines for a given indication. This may allow for treatment regimens with greater treatment effect, more durable response, and improved quality of life in indications with a high unmet medical need such as AML and CLL.

ZmELF6-ZmPRR37 module regulates maize flowering and salt response.

The control of flowering time in maize is crucial for reproductive success and yield, and it can be influenced by environmental stresses. Using the approaches of Ac/Ds transposon and transposable element amplicon sequencing techniques, we identified a Ds insertion mutant in the ZmPRR37 gene. The Ds insertion showed a significant correlation with days to anthesis. Further research indicated that ZmPRR37-CR knockout mutants exhibited early flowering, whereas ZmPRR37-overexpression lines displayed delayed flowering compared to WT under long-day (LD) conditions. We demonstrated that ZmPRR37 repressed the expression of ZmNF-YC2 and ZmNF-YA3 to delay flowering. Association analysis revealed a significant correlation between flowering time and a SNP2071-C/T located upstream of ZmPRR37. The SNP2071-C/T impacted the binding capacity of ZmELF6 to the promoter of ZmPRR37. ZmELF6 also acted as a flowering suppressor in maize under LD conditions. Notably, our study unveiled that ZmPRR37 can enhance salt stress tolerance in maize by directly regulating the expression of ABA-responsive gene ZmDhn1. ZmDhn1 negatively regulated maize salt stress resistance. In summary, our findings proposed a novel pathway for regulating photoperiodic flowering and responding to salt stress based on ZmPRR37 in maize, providing novel insights into the integration of abiotic stress signals into floral pathways.

Identification of IncA Plasmid, Harboring blaVIM-1 Gene, in S. enterica Goldcoast ST358 and C. freundii ST62 Isolated in a Hospitalized Patient.

In the present study, we analyzed the genome of two S. enterica strains TS1 and TS2 from stool and blood cultures, respectively, and one strain of C. freundii TS3, isolated from a single hospitalized patient with acute myeloid leukemia. The S. enterica Goldcoast ST358 (O:8 (C2-C3) serogroup), sequenced by the MiSeq Illumina system, showed the presence of β-lactamase genes (blaVIM-1, blaSHV-12 and blaOXA-10), aadA1, ant(2″)-Ia, aac(6')-Iaa, aac(6')-Ib3, aac(6')-Ib-cr, qnrVC6, parC(T57S), and several incompatibility plasmids. A wide variety of insertion sequences (ISs) and transposon elements were identified. In C. freundii TS3, these were the blaVIM-1, blaCMY-150, and blaSHV-12, aadA1, aac(6')-Ib3, aac(6')-Ib-cr, mph(A), sul1, dfrA14, ARR-2, qnrVC6, and qnrB38. IncA plasmid isolated from E.coli/K12 transconjugant and C. freundii exhibited a sequence identity >99.9%. The transfer of IncA plasmid was evaluated by conjugation experiments.

Systematic dissection of genomic features determining the vast diversity of conotoxins

BackgroundConus, a highly diverse species of venomous predators, has attracted significant attention in neuroscience and new drug development due to their rich collection of neuroactive peptides called conotoxins. Recent advancements in transcriptome, proteome, and genome analyses have facilitated the identification of conotoxins within Conus’ venom glands, providing insights into the genetic features and evolutionary patterns of conotoxin genes. However, the underlying mechanism behind the extraordinary hypervariability of conotoxins remains largely unknown.ResultsWe analyzed the transcriptomes of 34 Conus species, examining various tissues such as the venom duct, venom bulb, and salivary gland, leading to the identification of conotoxin genes. Genetic variation analysis revealed that a subset of these genes (15.78% of the total) in Conus species underwent positive selection (Ka/Ks > 1, p < 0.01). Additionally, we reassembled and annotated the genome of C. betulinus, uncovering 221 conotoxin-encoding genes. These genes primarily consisted of three exons, with a significant portion showing high transcriptional activity in the venom ducts. Importantly, the flanking regions and adjacent introns of conotoxin genes exhibited a higher prevalence of transposon elements, suggesting their potential contribution to the extensive variability observed in conotoxins. Furthermore, we detected genome duplication in C. betulinus, which likely contributed to the expansion of conotoxin gene numbers. Interestingly, our study also provided evidence of introgression among Conus species, indicating that interspecies hybridization may have played a role in shaping the evolution of diverse conotoxin genes.ConclusionsThis study highlights the impact of adaptive evolution and introgressive hybridization on the genetic diversity of conotoxin genes and the evolution of Conus. We also propose a hypothesis suggesting that transposable elements might significantly contribute to the remarkable diversity observed in conotoxins. These findings not only enhance our understanding of peptide genetic diversity but also present a novel approach for peptide bioengineering.

MegaLTR: a web server and standalone pipeline for detecting and annotating LTR-retrotransposons in plant genomes.

LTR-retrotransposons (LTR-RTs) are a class of RNA-replicating transposon elements (TEs) that can alter genome structure and function by moving positions, repositioning genes, shifting exons, and causing chromosomal rearrangements. LTR-RTs are widespread in many plant genomes and constitute a significant portion of the genome. Their movement and activity in eukaryotic genomes can provide insight into genome evolution and gene function, especially when LTR-RTs are located near or within genes. Building the redundant and non-redundant LTR-RTs libraries and their annotations for species lacking this resource requires extensive bioinformatics pipelines and expensive computing power to analyze large amounts of genomic data. This increases the need for online services that provide computational resources with minimal overhead and maximum efficiency. Here, we present MegaLTR as a web server and standalone pipeline that detects intact LTR-RTs at the whole-genome level and integrates multiple tools for structure-based, homologybased, and de novo identification, classification, annotation, insertion time determination, and LTR-RT gene chimera analysis. MegaLTR also provides statistical analysis and visualization with multiple tools and can be used to accelerate plant species discovery and assist breeding programs in their efforts to improve genomic resources. We hope that the development of online services such as MegaLTR, which can analyze large amounts of genomic data, will become increasingly important for the automated detection and annotation of LTR-RT elements.

Globally elevated levels of histone H3 lysine 9 trimethylation in early infancy are associated with poor growth trajectory in Bangladeshi children

BackgroundStunting is a global health problem affecting hundreds of millions of children worldwide and contributing to 45% of deaths in children under the age of five. Current therapeutic interventions have limited efficacy. Understanding the epigenetic changes underlying stunting will elucidate molecular mechanisms and likely lead to new therapies.ResultsWe profiled the repressive mark histone H3 lysine 9 trimethylation (H3K9me3) genome-wide in peripheral blood mononuclear cells (PBMCs) from 18-week-old infants (n = 15) and mothers (n = 14) enrolled in the PROVIDE study established in an urban slum in Bangladesh. We associated H3K9me3 levels within individual loci as well as genome-wide with anthropometric measurements and other biomarkers of stunting and performed functional annotation of differentially affected regions. Despite the relatively small number of samples from this vulnerable population, we observed globally elevated H3K9me3 levels were associated with poor linear growth between birth and one year of age. A large proportion of the differentially methylated genes code for proteins targeting viral mRNA and highly significant regions were enriched in transposon elements with potential regulatory roles in immune system activation and cytokine production. Maternal data show a similar trend with child’s anthropometry; however, these trends lack statistical significance to infer an intergenerational relationship.ConclusionsWe speculate that high H3K9me3 levels may result in poor linear growth by repressing genes involved in immune system activation. Importantly, changes to H3K9me3 were detectable before the overt manifestation of stunting and therefore may be valuable as new biomarkers of stunting.

Online data resource for exploring transposon insertion polymorphisms in public soybean germplasm accessions.

Soybean (Glycine max L. Merrill) is one of the most important economical crops. A large number of whole-genome resequencing datasets have been generated and are increasingly expanded for exploring genetic diversity and mining important quantitative trait loci. Most genome-wide association studies have focused on single-nucleotide polymorphisms, short insertions, and deletions. Nevertheless, structure variants mainly caused by transposon element mobilization are not fully considered. To fill this gap, we uniformly processed the publicly available whole-genome resequencing data from 5,521 soybean germplasm accessions and built an online soybean transposon insertion polymorphisms database named Soybean Transposon Insertion Polymorphisms Database (SoyTIPdb) (https://biotec.njau.edu.cn/soytipdb). The collected germplasm accessions derived from more than 45 countries and 160 regions representing the most comprehensive genetic diversity of soybean. SoyTIPdb implements easy-to-use query, analysis, and browse functions to help understand and find meaningful structural variations from TE insertions. In conclusion, SoyTIPdb is a valuable data resource and will help soybean breeders/researchers take advantage of the whole-genome sequencing datasets available in the public depositories.

Decoupling subgenomes within hybrid lavandin provide new insights into speciation and monoterpenoid diversification of Lavandula.

Polyploidization and transposon elements contribute to shape plant genome diversity and secondary metabolic variation in some edible crops. However, the specific contribution of these variations to the chemo-diversity of Lamiaceae, particularly in economic shrubs, is still poorly documented. The rich essential oils (EOs) of Lavandula plants are distinguished by monoterpenoids among the main EO-producing species, L. angustifolia (LA), L. × intermedia (LX) and L. latifolia (LL). Herein, the first allele-aware chromosome-level genome was assembled using a lavandin cultivar 'Super' and its hybrid origin was verified by two complete subgenomes (LX-LA and LX-LL). Genome-wide phylogenetics confirmed that LL, like LA, underwent two lineage-specific WGDs after the γ triplication event, and their speciation occurred after the last WGD. Chloroplast phylogenetic analysis indicated LA was the maternal source of 'Super', which produced premium EO (higher linalyl/lavandulyl acetate and lower 1,8-cineole and camphor) close to LA. Gene expression, especially the monoterpenoid biosynthetic genes, showed bias to LX-LA alleles. Asymmetric transposon insertions in two decoupling 'Super' subgenomes were responsible for speciation and monoterpenoid divergence of the progenitors. Both hybrid and parental evolutionary analysis revealed that LTR (long terminal repeat) retrotransposon associated with AAT gene loss cause no linalyl/lavandulyl acetate production in LL, and multi-BDH copies retained by tandem duplication and DNA transposon resulted in higher camphor accumulation of LL. Advances in allelic variations of monoterpenoids have the potential to revolutionize future lavandin breeding and EO production.

Study on the role of \(\textit{TDRD1}\) variants in male infertility among 310 Vietnamese individuals

Infertility is a global reproductive health burden affected by various genetic factors, including spermatogenic defects. The instability of the germ cells’ genome caused by the unregulated replication of transposable elements is one of the causes of spermatogenic impairment. Tudor domain-containing 1 (TDRD1) expressed only in germ cells plays a significant role in the piRNA (PIWI-interacting RNA) pathway to maintain genome integrity via suppressing transposon elements during spermatogenesis. Despite the protein’s role in male germline development, TDRD1 has not been studied intensively. In this study, we established the relationship between male infertility and single nucleotide polymorphisms (SNPs) of TDRD1 (rs541192490, rs77559927) among 310 Vietnamese men (160 infertile patients and 150 healthy controls). Genotypes of single nucleotide polymorphisms of the TDRD1 gene (SNPs) were identified using the PCR-RFLP method. The results showed that TDRD1 SNPs were not associated with male infertility in all three test models (additive, dominant, and recessive) (p-value &gt; 0.05). Haplotype analysis of the two SNPs also showed similar findings. This study would contribute to the knowledge of TDRD1’s association with male infertility in the Vietnamese population.

Phenotypic and Genetic Analysis of KPC-49, a KPC-2 Variant Conferring Resistance to Ceftazidime-Avibactam and Maintaining Resistance to Imipenem and Meropenem.

Klebsiella pneumoniae, a gram-negative bacterium, poses a severe hazard to public health, with many bacterial hosts having developed resistance to most antibiotics in clinical use. The goal of this study was to look into the development of resistance to both ceftazidime-avibactam and carbapenems, including imipenem and meropenem, in a K. pneumonia strain expressing a novel K. pneumoniae carbapenemase-2 (KPC-2) variant, referred to as KPC-49. After 1 day of incubation of K1 on agar containing ceftazidime-avibactam (MIC = 16/4 mg/L), a second KPC-producing K. pneumoniae strain (K2) was recovered. Antimicrobial susceptibility assays, cloning assays, and whole genome sequencing were performed to analyse and evaluate antibiotic resistance phenotypes and genotypes. K. pneumoniae strain (K1), that produced KPC-2, was susceptible to ceftazidime-avibactam but resistant to carbapenems. The K2 isolate harboured a novel bla KPC-49 variant, which differs from bla KPC-2 by a single nucleotide (C487A), and results in an arginine-serine substitution at amino acid position 163 (R163S). The mutant K2 strain was resistant to both ceftazidime-avibactam and carbapenems. We demonstrated the ability of KPC-49 to hydrolyse carbapenems, which may be attributed to high KPC-49 expression or presence of an efflux pump and/or absence of membrane pore proteins in K2. Furthermore, blaKPC-like was carried on an IncFII (pHN7A8)/IncR-type plasmid within a TnAs1-orf-orf-orf-orf-orf-orf-ISKpn6-bla KPC-ISKpn27 structure. The bla KPC-like gene was flanked by various insertion sequences and transposon elements, including the Tn3 family transposon, such as TnAs1, TnAs3, IS26, and IS481-ISKpn27. New KPC variants are emerging owing to sustained exposure to antimicrobials and modifications in their amino acid sequences. We demonstrated the drug resistance mechanisms of the new mutant strains through experimental whole genome sequencing combined with bioinformatics analysis. Enhanced understanding of laboratory and clinical features of infections due to K. pneumoniae of the new KPC subtype is key to early and accurate anti-infective therapy.

Mining metagenomes reveals diverse antibiotic biosynthetic genes in uncultured microbial communities

Pathogens resistant to antimicrobials form a significant threat to public health worldwide. Tackling multidrug-resistant pathogens via screening metagenomic libraries has become a common approach for the discovery of new antibiotics from uncultured microorganisms. This study focuses on capturing nonribosomal peptide synthase (NRPS) gene clusters implicated in the synthesis of many natural compounds of industrial relevance. A NRPS PCR assay was used to screen 2976 Escherichia coli clones in a soil metagenomic library to target NRPS genes. DNA extracts from 4 clones were sequenced and subjected to bioinformatic analysis to identify NRPS domains, their phylogeny, and substrate specificity.Successfully, 17 NRPS-positive hits with a biosynthetic potential were identified. DNA sequencing and BLAST analysis confirmed that NRPS protein sequences shared similarities with members of the genus Delftia in the Proteobacteria taxonomic position. Multiple alignment and phylogenetic analysis demonstrated that clones no. 15cd35 and 15cd37 shared low bootstrap values (54%) and were distantly far from close phylogenetic neighbors. Additionally, NRPS domain substrate specificity has no hits with the known ones; hence, they are more likely to use different substrates to produce new diverse antimicrobials. Further analysis confirmed that the NRPS hits resemble several transposon elements from other bacterial taxa, confirming its diversity. We confirmed that the analyses of the soil metagenomic library revealed a diverse set of NRPS related to the genus Delftia. An in-depth understanding of those positive NRPS hits is a crucial step for genetic manipulation of NRPS, shedding light on alternative novel antimicrobial compounds that can be used in drug discovery and hence supports the pharmaceutical sector.