Sequence Repeat Research Articles

Sequencing-guided design of genetically encoded small RNAs targeting CAG repeats for selective inhibition of mutant huntingtin

Huntington's disease (HD) is an incurable neurodegenerative disorder caused by genetic expansion of a CAG repeat sequence in one allele of the huntingtin (HTT) gene. Reducing expression of the mutant HTT (mutHTT) protein has remained a clear therapeutic goal, but reduction of wild-type HTT (wtHTT) is undesirable, as it compromises gene function and potential therapeutic efficacy. One promising allele-selective approach involves targeting the CAG repeat expansion with steric binding small RNAs bearing central mismatches. However, successful genetic encoding requires consistent placement of mismatches to the target within the small RNA guide sequence, which involves 5' processing precision by cellular enzymes. Here, we used small RNA sequencing (RNA-seq) to monitor the processing precision of a limited set of CAG repeat-targeted small RNAs expressed from multiple scaffold contexts. Small RNA-seq identified expression constructs with high-guide strand 5' processing precision and promising allele-selective inhibition of mutHTT. Transcriptome-wide mRNA-seq also identified an allele-selective small RNA with a favorable off-target profile. These results support continued investigation and optimization of genetically encoded repeat-targeted small RNAs for allele-selective HD gene therapy and underscore the value of sequencing methods to balance specificity with allele selectivity during the design and selection process.

Elucidating the multichromosomal structure within the Brasenia schreberi mitochondrial genome through assembly and analysis

Brasenia schreberi, a plant species traditionally utilized in Chinese medicine and cuisine, represents an early evolutionary stage among flowering plants (angiosperms). While the plastid genome of this species has been published, its mitochondrial genome (mitogenome) has not been extensively explored, with a notable absence of thorough comparative analyses of its organellar genomes. In our study, we had assembled the entire mitogenome of B. schreberi utilizing the sequencing data derived from both Illumina platform and Oxford Nanopore. The B. schreberi mitogenome mostly exists as six circular DNA molecules, with the largest being 628,257 base pairs (bp) and the smallest 110,220 bp, amounting to 1.49 megabases (Mb). Then we annotated the mitogenome of B. schreberi. The mitogenome encompasses a total of 71 genes: 40 of these are coding proteins genes (PCGs), 28 are genes for transfer RNA (tRNA), and the remaining 3 are genes for ribosomal RNA (rRNA). In the analysis of codon usage, we noted a unique codon preference specific to each amino acid. The most commonly used codons exhibited an average RSCU of 1.36, indicating a noticeable bias in codon selection. In the repeat sequence analysis, a total of 553 simple sequence repeats (SSRs) were identified, 1,822 dispersed repeats (comprising 1,015 forward and 807 palindromic repeats), and 608 long terminal repeats (LTRs). Additionally, in the analysis of homologous sequences between organelle genomes, we detected 38 homologous sequences derived from the plastid genome, each exceeding 500 bp, within the B. schreberi mitochondrial genome. Notably, ten tRNA genes (trnC-GCA, trnM-CAU, trnI-CAU, trnQ-UUG, trnN-GUU, trnT-GGU, trnW-CCA, trnA-UGC, trnI-GAU, and trnV-GAC) appear to have been completely transferred from the chloroplast to the mitogenome. Utilizing the Deepred-mt to predict the RNA editing sites in the mitogenome, we have identified 675 high-quality RNA editing sites in the 40 mitochondrial PCGs. In the final stage of our study, we performed an analysis of colinearity and inferred the phylogenetic relationship of B. schreberi with other angiosperms, utilizing the mitochondrial PCGs as a basis. The results showed that the non-coding regions of the B. schreberi mitogenome are characterized by an abundance of repetitive sequences and exogenous sequences, and B. schreberi is more closely related with Euryale ferox.

Roles of Nuclear Orphan Receptors TR2 and TR4 during Hematopoiesis.

TR2 and TR4 (NR2C1 and NR2C2, respectively) are evolutionarily conserved nuclear orphan receptors capable of binding direct repeat sequences in a stage-specific manner. Like other nuclear receptors, TR2 and TR4 possess important roles in transcriptional activation or repression with developmental stage and tissue specificity. TR2 and TR4 bind DNA and possess the ability to complex with available cofactors mediating developmental stage-specific actions in primitive and definitive erythrocytes. In erythropoiesis, TR2 and TR4 are required for erythroid development, maturation, and key erythroid transcription factor regulation. TR2 and TR4 recruit and interact with transcriptional corepressors or coactivators to elicit developmental stage-specific gene regulation during hematopoiesis.

Solanum aculeatissimum and Solanum torvum chloroplast genome sequences: a comparative analysis with other Solanum chloroplast genomes

BackgroundSolanum aculeatissimum and Solanum torvum belong to the Solanum species, and they are essential plants known for their high resistance to diseases and adverse conditions. They are frequently used as rootstocks for grafting and are often crossbred with other Solanum species to leverage their resistance traits. However, the phylogenetic relationship between S. aculeatissimum and S. torvum within the Solanum genus remains unclear. Therefore, this paper aims to sequence the complete chloroplast genomes of S. aculeatissimum and S. torvum and analyze them in comparison with 29 other previously published chloroplast genomes of Solanum species.ResultsWe observed that the chloroplast genomes of S. aculeatissimum and S. torvum possess typical tetrameric structures, consisting of one Large Single Copy (LSC) region, two reverse-symmetric Inverted Repeats (IRs), and one Small Single Copy (SSC) region. The total length of these chloroplast genomes ranged from 154,942 to 156,004 bp, with minimal variation. The highest GC content was found in the IR region, while the lowest was in the SSC region. Regarding gene content, the total number of chloroplast genes and CDS genes remained relatively consistent, ranging from 128 to 134 and 83 to 91, respectively. Nevertheless, there was notable variability in the number of tRNA genes and rRNAs. Relative synonymous codon usage (RSCU) analysis revealed that both S. aculeatissimum and S. torvum preferred codons that utilized A and U bases. Analysis of the IR boundary regions indicated that contraction and expansion primarily occurred at the junction between SSC and IR regions. Nucleotide polymorphism analysis and structural variation analysis demonstrated that chloroplast variation in Solanum species mainly occurred in the LSC and SSC regions. Repeat sequence analysis revealed that A/T was the most frequent base pair in simple repeat sequences (SSR), while Palindromic and Forward repeats were more common in long sequence repeats (LSR), with Reverse and Complement repeats being less frequent. Phylogenetic analysis indicated that S. aculeatissimum and S. torvum belonged to the same meristem and were more closely related to Cultivated Eggplant.ConclusionThese findings enhance our comprehension of chloroplast genomes within the Solanum genus, offering valuable insights for plant classification, evolutionary studies, and potential molecular markers for species identification.

Complete Plastomes of Ten Rorippa Species (Brassicaceae): Comparative Analysis and Phylogenetic Relationships

The genus Rorippa belongs to the family Brassicaceae, and its members usually have high medicinal value. The genus consists of approximately 75 species and mainly grows in the Northern Hemisphere, occurring in every continent except Antarctica. The taxonomy and phylogenetic relationships of Rorippa are still unsettled, largely due to complex morphological variations in Rorippa, which were caused by frequent hybridization events. Here, we sequenced four complete plastid genomes of Rorippa species by Illumina paired-end sequencing. The four new plastid genomes of Rorippa ranged in total size from 154,671 bp for R. palustris to 154,894 bp for R. sylvestris. There are 130 genes in the four plastomes, embodying 8 rRNA, 37 tRNA, and 85 protein-coding genes. Combining with six published plastid genomes, we carried on comparative and phylogenetic analyses. We found that the ten Rorippa plastid genomes were conservative in gene number and order, total size, genomic structure, codon usage, long repeat sequence, and SSR. Fourteen mutational hotspot regions could be selected as candidate DNA barcoding to distinguish Rorippa plants. The phylogenetic trees clearly identified that ten Rorippa species displayed monophyletic relationships within the tribe Cardamineae based on plastomes and nrDNA ITS sequences. However, there are significant cytonuclear discordances in the interspecific relationships within Rorippa, as well as the intergeneric relationships between Rorippa and its related genera. We inferred that the cytonuclear discordance is most likely a result of interspecific hybridization within Rorippa, as well as intergeneric hybridization with its related genera. These plastid genomes can offer precious information for studies of species authentication, evolutionary history, and the phylogeny of Rorippa.

Inhibition of sugar-binding activity of Galectins-8 by thiogalactoside (TDG) attenuates secondary brain damage and improves long-term prognosis following intracerebral hemorrhage

Galectins-8 (Gal-8), the tandem repeat sequences of the galectin family, can influence the pathophysiologic processes in neurological disorders. However, its effect on intracerebral hemorrhage and related mechanisms remains nebulous. Using collagenase VII-S-induced ICH in the left striatum of mice, we investigated the effects of Gal-8 on cellular and molecular immune inflammatory responses in hemorrhagic brain and evaluated the severity of short- and long-term brain injury. Our results showed that activated microglia in the periphery of hematoma in mice with intracerebral hemorrhage expressed Gal-8, while Gal-8 could regulate the expression of cytokines, such as HMGB-1 (P = 0.0032), TNF-α (P = 0.0158), and IL-10 (P = 0.0379). Inhibition of the glucose-binding activity of Gal-8 by thiogalactoside (TDG) significantly reduced the volume of cerebral hematoma (P = 0.0241) and hydrocephalus (P = 0.0112) during the acute phase of cerebral hemorrhage and improved the long-term prognosis. TDG can reduce acute-phase brain tissue injury and improve the prognosis by inhibiting the activation of immune-inflammatory cells in the periphery of hematoma and reducing the release of pro-inflammatory factors.

Precise editing of pathogenic nucleotide repeat expansions in iPSCs using paired prime editor.

Nucleotide repeat expansion disorders, a group of genetic diseases characterized by the expansion of specific DNA sequences, pose significant challenges to treatment and therapy development. Here, we present a precise and programmable method called prime editor-mediated correction of nucleotide repeat expansion (PE-CORE) for correcting pathogenic nucleotide repeat expansion. PE-CORE leverages a prime editor and paired pegRNAs to achieve targeted correction of repeat sequences. We demonstrate the effectiveness of PE-CORE in HEK293T cells and patient-derived induced pluripotent stem cells (iPSCs). Specifically, we focus on spinal and bulbar muscular atrophy and spinocerebellar ataxia type, two diseases associated with nucleotide repeat expansion. Our results demonstrate the successful correction of pathogenic expansions in iPSCs and subsequent differentiation into motor neurons. Specifically, we detect distinct downshifts in the size of both the mRNA and protein, confirming the functional correction of the iPSC-derived motor neurons. These findings highlight PE-CORE as a precision tool for addressing the intricate challenges of nucleotide repeat expansion disorders, paving the way for targeted therapies and potential clinical applications.

A point-of-care testing assay for clonorchiasis using a EuNPs-CsTR1 fluorescent probe-based immunoassay.

Clonorchis sinensis is one of the most important fish-borne zoonotic parasitic worms in humans, and is distributed in several countries with more than 15 million people infected globally. However, the lack of a point-of-care testing (POCT) method is still the critical barrier to effectively prevent clonorchiasis. With the application of novel fluorescent nanomaterials, the development of on-site testing methods with high signal enhancement can provide a simple, precise and inexpensive tool for disease detection. In this study, Eu-(III) nanoparticles (EuNPs) were used as indicative probes, combined with C. sinensis tandem repeat sequence 1 (CSTR1) antigen to capture specific antibodies. Afterward, the complex binds to mouse anti-human IgG immobilized on the test line (T-line) producing a fluorescent signal under UV light. The EuNPs-fluorescent immunoassay (EuNPs-FIA) was successfully constructed, allowing sample detection within 10 min. It enabled both qualitative determination with the naked eye under UV light and quantitative detection by scanning the fluorescence intensity on the test line and control line (C-line). A total of 133 clinical human sera (74 negative, 59 clonorchiasis, confirmed by conventional Kato-Katz (KK) methods and PCR via testing fecal samples corresponding to each serum sample) were used in this study. For qualitative analysis, the cut-off value of fluorescence for positive serum was 31.57 by testing 74 known negative human samples. The assay had no cross-reaction with other 9 parasite-infected sera, and could recognize the mixed infection sera of C. sinensis and other parasites. The sensitivity and specificity of EuNPs-FIA were both 100% compared with KK smear method. Taking advantage of its high precision and user-friendly procedure, the established EuNPs-FIA provides a powerful tool for the diagnosis and epidemiological survey of clonorchiasis.

Genome organization of satellite DNA beta molecule associated with Mungbean yellow mosaic India virus Gujarat, India

Recent advancements in genomics have unveiled fascinating insights into the intricate world of plant pathogens. The discovery of an intriguing association between the DNA beta and MYMIV-[IN:Ana: CpMBKA25:04], a strain containing 1367 base pairs of nucleotides, has sparked a wave of scientific curiosity. Through sequence identity matrix analysis, it was revealed that this strain shares a striking 96% similarity with the Potato apical leaf curl virus, closely followed by the Tomato leaf curl beta satellite virus [Accession No. AY230138] at 95.3%. In contrast, its relationship with the beta satellite Bhendi yellow vein beta satellite [Accession No. AJ308425] stood at a comparative 45.7% when compared against seventeen other DNA beta isolates infecting various crops. Noteworthy features of the beta satellite sequences include a distinctive putative stem-loop structure TAATATT/AC in the common region and a recurring GCTACGC repeat sequence. An A-rich region spanning from nucleotides 910 to 1120 further adds to the genetic complexity of this novel strain. Intriguingly, an analysis of the amino acid sequence unveiled additional residues V, C, and N at the carboxyl end of the βC1 protein, setting it apart from its counterparts. The identity matrix analysis established a 95.9% similarity between βC1 and the beta satellite associated with Tomato leaf curl Karnataka virus and Potato apical leaf curl virus, phylogenetic analysis placed the beta satellite linked with MYMIV-[IN:Ana: CpMBKA25:04] in a distinct category within the Papaya leaf curl virus family, shedding light on its unique evolutionary lineage and biological significance in the realm of plant virology.

A telomere-to-telomere gap-free reference genome assembly of avocado provides useful resources for identifying genes related to fatty acid biosynthesis and disease resistance.

Avocado (Persea americana Mill.) is an economically valuable plant because of the high fatty acid content and unique flavor of its fruits. Its fatty acid content, especially the relatively high unsaturated fatty acid content, provides significant health benefits. We herein present a telomere-to-telomere gapless genome assembly (841.6Mb) of West Indian avocado. The genome contains 40 629 predicted protein-coding genes. Repeat sequences account for 57.9% of the genome. Notably, all telomeres, centromeres, and a nucleolar organizing region are included in this genome. Fragments from these three regions were observed via fluorescence in situ hybridization. We identified 376 potential disease resistance-related nucleotide-binding leucine-rich repeat genes. These genes, which are typically clustered on chromosomes, may be derived from gene duplication events. Five NLR genes (Pa11g0262, Pa02g4855, Pa07g3139, Pa07g0383, and Pa02g3196) were highly expressed in leaves, stems, and fruits, indicating they may be involved in avocado disease responses in multiple tissues. We also identified 128 genes associated with fatty acid biosynthesis and analyzed their expression patterns in leaves, stems, and fruits. Pa02g0113, which encodes one of 11 stearoyl-acyl carrier protein desaturases mediating C18 unsaturated fatty acid synthesis, was more highly expressed in the leaves than in the stems and fruits. These findings provide valuable insights that enhance our understanding of fatty acid biosynthesis in avocado.

The first chromosome-level genome of the stag beetle Dorcus hopei Saunders, 1854 (Coleoptera: Lucanidae)

Stag beetles (Coleoptera: Lucanidae) represent a significant saproxylic assemblage in forest ecosystems and are noted for their enlarged mandibles and male polymorphism. Despite their relevance as ideal models for the study of exaggerated mandibles that aid in attracting mates, the regulatory mechanisms associated with these traits remain understudied, and restricted by the lack of high-quality reference genomes for stag beetles. To address this limitation, we successfully assembled the first chromosome-level genome of a representative species Dorcus hopei. The genome was 496.58 Mb in length, with a scaffold N50 size of 54.61 Mb, BUSCO values of 99.8%, and 96.8% of scaffolds anchored to nine pairs of chromosomes. We identified 285.27 Mb (57.45%) of repeat sequences and annotated 11,231 protein-coding genes. This genome will be a valuable resource for further understanding the evolution and ecology of stag beetles, and provides a basis for studying the mechanisms of exaggerated mandibles through comparative analysis.

Proximal telomeric decompaction due to telomere shortening drives FOXC1-dependent myocardial senescence.

Telomeres, TTAGGGn DNA repeat sequences located at the ends of eukaryotic chromosomes, play a pivotal role in aging and are targets of DNA damage response. Although we and others have demonstrated presence of short telomeres in genetic cardiomyopathic and heart failure cardiomyocytes, little is known about the role of telomere lengths in cardiomyocyte. Here, we demonstrate that in heart failure patient cardiomyocytes, telomeres are shortened compared to healthy controls. We generated isogenic human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) with short telomeres (sTL-CMs) and normal telomeres (nTL-CMs) as model. Compared to nTL-CMs, short telomeres result in cardiac dysfunction and expression of senescent markers. Using Hi-C and RNASeq, we observe that short telomeres induced TAD insulation decrease near telomeric ends and this correlated with a transcription upregulation in sTL-CMs. FOXC1, a key transcription factor involved in early cardiogenesis, was upregulated in sTL-CMs and its protein levels were negatively correlated with telomere lengths in heart failure patients. Overexpression of FOXC1 induced hiPSC-CM aging, mitochondrial and contractile dysfunction; knockdown of FOXC1 rescued these phenotypes. Overall, the work presented demonstrate that increased chromatin accessibility due to telomere shortening resulted in the induction of FOXC1-dependent expression network responsible for contractile dysfunction and myocardial senescence.

Whole-Genome Sequencing Analyses Reveal the Evolution Mechanisms of Typical Biological Features of Decapterus maruadsi.

Decapterus maruadsi is a typical representative of small pelagic fish characterized by fast growth rate, small body size, and high fecundity. It is a high-quality marine commercial fish with high nutritional value. However, the underlying genetics and genomics research focused on D. maruadsi is not comprehensive. Herein, a high-quality chromosome-level genome of a male D. maruadsi was assembled. The assembled genome length was 716.13 Mb with contig N50 of 19.70 Mb. Notably, we successfully anchored 95.73% contig sequences into 23 chromosomes with a total length of 685.54 Mb and a scaffold N50 of 30.77 Mb. A total of 22,716 protein-coding genes, 274.90 Mb repeat sequences, and 10,060 ncRNAs were predicted, among which 22,037 (97%) genes were successfully functionally annotated. The comparative genome analysis identified 459 unique, 73 expanded, and 52 contracted gene families. Moreover, 2804 genes were identified as candidates for positive selection, of which some that were related to the growth and development of bone, muscle, cardioid, and ovaries, such as some members of the TGF-β superfamily, were likely involved in the evolution of typical biological features in D. maruadsi. The study provides an accurate and complete chromosome-level reference genome for further genetic conservation, genomic-assisted breeding, and adaptive evolution research for D. maruadsi.

Assembly of continuous high-resolution draft genome sequence of Hemicentrotus pulcherrimus using long-read sequencing.

The update of the draft genome assembly of sea urchin, Hemicentrotus pulcherrimus, which is widely studied in East Asia as a model organism of early development, was performed using Oxford nanopore long-read sequencing. The updated assembly provided ~600-Mb genome sequences divided into 2,163 contigs with N50 = 516 kb. BUSCO completeness score and transcriptome model mapping ratio (TMMR) of the present assembly were obtained as 96.5% and 77.8%, respectively. These results were more continuous with higher resolution than those by the previous version of H. pulcherrimus draft genome, HpulGenome_v1, where the number of scaffolds = 16,251 with a total of ~100 Mb, N50 = 143 kb, BUSCO completeness score = 86.1%, and TMMR = 55.4%. The obtained genome contained 36,055 gene models that were consistent with those in other echinoderms. Additionally, two tandem repeat sequences of early histone gene locus containing 47 copies and 34 copies of all histone genes, and 185 of the homologous sequences of the interspecifically conserved region of the Ars insulator, ArsInsC, were obtained. These results provide further advance for genome-wide research of development, gene regulation, and intranuclear structural dynamics of multicellular organisms using H. pulcherrimus.

Cdc13 exhibits dynamic DNA strand exchange in the presence of telomeric DNA.

Telomerase is the enzyme that lengthens telomeres and is tightly regulated by a variety of means to maintain genome integrity. Several DNA helicases function at telomeres, and we previously found that the Saccharomyces cerevisiae helicases Hrq1 and Pif1 directly regulate telomerase. To extend these findings, we are investigating the interplay between helicases, single-stranded DNA (ssDNA) binding proteins (ssBPs), and telomerase. The yeast ssBPs Cdc13 and RPA differentially affect Hrq1 and Pif1 helicase activity, and experiments to measure helicase disruption of Cdc13/ssDNA complexes instead revealed that Cdc13 can exchange between substrates. Although other ssBPs display dynamic binding, this was unexpected with Cdc13 due to the reported in vitro stability of the Cdc13/telomeric ssDNA complex. We found that the DNA exchange by Cdc13 occurs rapidly at physiological temperatures, requires telomeric repeat sequence DNA, and is affected by ssDNA length. Cdc13 truncations revealed that the low-affinity binding site (OB1), which is distal from the high-affinity binding site (OB3), is required for this intermolecular dynamic DNA exchange (DDE). We hypothesize that DDE by Cdc13 is the basis for how Cdc13 'moves' at telomeres to alternate between modes where it regulates telomerase activity and assists in telomere replication.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

The complete Chloroplast genome of Stachys geobombycis and comparative analysis with related Stachys species

Herb genomics, at the forefront of traditional Chinese medicine research, combines genomics with traditional practices, facilitating the scientific validation of ancient remedies. This integration enhances public understanding of traditional Chinese medicine’s efficacy and broadens its scope in modern healthcare. Stachys species encompass annual or perennial herbs or small shrubs, exhibiting simple petiolate or sessile leaves. Despite their wide-ranging applications across various fields, molecular data have been lacking, hindering the precise identification and taxonomic elucidation of Stachys species. To address this gap, we assembled the complete chloroplast (CP) genome of Stachys geobombycis and conducted reannotation and comparative analysis of seven additional species within the Stachys genus. The findings demonstrate that the CP genomes of these species exhibit quadripartite structures, with lengths ranging from 14,523 to 150,599 bp. Overall, the genome structure remains relatively conserved, hosting 131 annotated genes, including 87 protein coding genes, 36 tRNA genes, and 8 rRNA genes. Additionally, 78 to 98 SSRs and long repeat sequences were detected , and notably, 6 highly variable regions were identified as potential molecular markers in the CP genome through sequence alignment. Phylogenetic analysis based on Bayesian inference and maximum likelihood methods strongly supported the phylogenetic position of the genus Stachys as a member of Stachydeae tribe. Overall, this comprehensive bioinformatics study of Stachys CP genomes lays the groundwork for phylogenetic classification, plant identification, genetic engineering, evolutionary studies, and breeding research concerning medicinal plants within the Stachys genus.

Re-Sequencing the Mitochondrial Genome Unveils a Novel Isomeric Form of NWB CMS Line in Radish and Functional Verification of Its Candidate Sterile Gene

Radish (Raphanus sativus L.) is a globally significant vegetable and relies on cytoplasmic male sterile (CMS) lines for hybrid seed production. The NWB CMS type is favored over Ogura CMS for its ease in maintainer screening. Despite its varied mitochondrial configurations and unvalidated sterile gene, we re-sequenced the mitochondrial genome of NWB CMS Tibet A and verified the function of the sterility gene via genetic transformation of Arabidopsis thaliana. The mitochondrial genomes of Tibet A could be assembled into circular DNA molecules, with a mitochondrial genome size of 239,184 bp. Our analysis indicated that the specific orf463a was the CMS-associated gene in Tibet A, sharing sequence consistency with the CMS gene in DCGMS and NWB CMS YB-A. Collinearity analysis showed that the mitochondrial genomes of NWB CMS Tibet A, DCGMS, and NWB CMS YB-A share the same mitotype, with structural variations due to recombination via 9731 bp long repeat sequences and 508 bp short repeat sequences. Driven by the Ap3 promoter, transgenic Arabidopsis with orf463a exhibited male sterility, confirming the gene’s potential role in CMS. In this study, we assembled a new isomeric form of NWB CMS mitochondrial genome and proved the function of the candidate sterile gene.

Telomere-to-telomere genome assembly of Oldenlandia diffusa.

We report the complete telomere-to-telomere genome assembly of Oldenlandia diffusa which renowned in traditional Chinese medicine, comprising 16 chromosomes and spanning 499.7Mb. The assembly showcases 28 telomeres and minimal gaps, with a total of only five. Repeat sequences constitute 46.41% of the genome, and 49,701 potential protein-coding genes have been predicted. Compared with O. corymbosa, O. diffusa exhibits chromosome duplication and fusion events, diverging 20.34 million years ago. Additionally, a total of 11 clusters of terpene synthase have been identified. The comprehensive genome sequence, gene catalog, and terpene synthase clusters of O. diffusa detailed in this study will significantly contribute to advancing research in this species' genetic, genomic, and pharmacological aspects.

Mitochondrial genome variation and intergenomic sequence transfers in Hevea species.

Among the Hevea species, rubber tree (Hevea brasiliensis) is the most important source of natural rubber. In previous studies, we sequenced the complete nuclear and chloroplast genomes of Hevea species, providing an invaluable resource for studying their phylogeny, disease resistance, and breeding. However, given that plant mitochondrial genomes are more complex and more difficult to assemble than that of the other organelles, little is known about their mitochondrial genome, which limits the comprehensive understanding of Hevea genomic evolution. In this study, we sequenced and assembled the mitochondrial genomes of four Hevea species. The four mitochondrial genomes had consistent GC contents, codon usages and AT skews. However, there were significant differences in the genome lengths and sequence repeats. Specifically, the circular mitochondrial genomes of the four Hevea species ranged from 935,732 to 1,402,206 bp, with 34-35 unique protein-coding genes, 35-38 tRNA genes, and 6-13 rRNA genes. In addition, there were 17,294-46,552 bp intergenomic transfer fragments between the chloroplast and mitochondrial genomes, consisting of eight intact genes (psaA, rrn16S, tRNA-Val, rrn5S, rrn4.5S, tRNA-Arg, tRNA-Asp, and tRNA-Asn), intergenic spacer regions and partial gene sequences. The evolutionary position of Hevea species, crucial for understanding its adaptive strategies and relation to other species, was verified by phylogenetic analysis based on the protein-coding genes in the mitochondrial genomes of 21 Malpighiales species. The findings from this study not only provide valuable insights into the structure and evolution of the Hevea mitochondrial genome but also lay the foundation for further molecular, evolutionary studies, and genomic breeding studies on rubber tree and other Hevea species, thereby potentially informing conservation and utilization strategies.