Subtelomeric Regions Research Articles

Genomic localization bias of secondary metabolite gene clusters and association with histone modifications in Aspergillus

Abstract Fungi are well-known producers of bioactive secondary metabolites (SMs), which have been exploited for decades by humankind for various medical applications like therapeutics and antibiotics. SMs are synthesized by biosynthetic gene clusters (BGCs) – physically co-localized and co-regulated genes. Because BGCs are often regulated by histone post-translational modifications (PTMs), it was suggested that their chromosomal location is important for their expression. Studies in a few fungal species indicated an enrichment of BGCs in sub-telomeric regions; however, there is no evidence that BGCs with distinct genomic localization are regulated by different histone PTMs. Here, we used 174 Aspergillus species covering 22 sections to determine the correlation between BGC genomic localization, gene expression and histone PTMs. We found a high abundance and diversity of SM backbone genes across the Aspergillus genus, with notable unique genes within sections. Being unique or conserved in many species, BGCs showed a strong bias for being localized in low-synteny regions, regardless of their position in chromosomes. Using chromosome-level assemblies, we also confirmed a significantly biased localization in sub-telomeric regions. Notably, SM backbone genes in sub-telomeric regions and about half of those in low-synteny regions exhibit higher gene expression variability, likely due to the similar higher variability in H3K4me3 and H3K36me3 histone PTMs; while variations in histone H3 acetylation and H3K9me3 are not correlated to genomic localization and expression variation, as analyzed in two Aspergillus species. Expression variability across four Aspergillus species further supports that BGCs tend to be located in low-synteny regions and that regulation of expression in those regions likely involves different histone PTMs than the most commonly studied modifications.

A Plasmodium falciparum MORC protein complex modulates epigenetic control of gene expression through interaction with heterochromatin

Dynamic control of gene expression is critical for blood stage development of malaria parasites. Here, we used multi-omic analyses to investigate transcriptional regulation by the chromatin-associated microrchidia protein, MORC, during asexual blood stage development of the human malaria parasite Plasmodium falciparum. We show that PfMORC (PF3D7_1468100) interacts with a suite of nuclear proteins, including APETALA2 (ApiAP2) transcription factors (PfAP2-G5, PfAP2-O5, PfAP2-I, PF3D7_0420300, PF3D7_0613800, PF3D7_1107800, and PF3D7_1239200), a DNA helicase DS60 (PF3D7_1227100), and other chromatin remodelers (PfCHD1 and PfEELM2). Transcriptomic analysis of PfMORCHA-glmS knockdown parasites revealed 163 differentially expressed genes belonging to hypervariable multigene families, along with upregulation of genes mostly involved in host cell invasion. In vivo genome-wide chromatin occupancy analysis during both trophozoite and schizont stages of development demonstrates that PfMORC is recruited to repressed, multigene families, including the var genes in subtelomeric chromosomal regions. Collectively, we find that PfMORC is found in chromatin complexes that play a role in the epigenetic control of asexual blood stage transcriptional regulation and chromatin organization.

A Plasmodium falciparum MORC protein complex modulates epigenetic control of gene expression through interaction with heterochromatin.

Dynamic control of gene expression is critical for blood stage development of malaria parasites. Here, we used multi-omic analyses to investigate transcriptional regulation by the chromatin-associated microrchidia protein, MORC, during asexual blood stage development of the human malaria parasite Plasmodium falciparum. We show that PfMORC (PF3D7_1468100) interacts with a suite of nuclear proteins, including APETALA2 (ApiAP2) transcription factors (PfAP2-G5, PfAP2-O5, PfAP2-I, PF3D7_0420300, PF3D7_0613800, PF3D7_1107800, and PF3D7_1239200), a DNA helicase DS60 (PF3D7_1227100), and other chromatin remodelers (PfCHD1 and PfEELM2). Transcriptomic analysis of PfMORCHA-glmS knockdown parasites revealed 163 differentially expressed genes belonging to hypervariable multigene families, along with upregulation of genes mostly involved in host cell invasion. In vivo genome-wide chromatin occupancy analysis during both trophozoite and schizont stages of development demonstrates that PfMORC is recruited to repressed, multigene families, including the var genes in subtelomeric chromosomal regions. Collectively, we find that PfMORC is found in chromatin complexes that play a role in the epigenetic control of asexual blood stage transcriptional regulation and chromatin organization.

Units containing telomeric repeats are prevalent in subtelomeric regions of a Mesorhabditis isolate collected from the Republic of Korea.

Mesorhabditis is known for its somatic genome being only a small portion of the germline genome due to programmed DNA elimination. This phenotype may be associated with the maintenance of telomeres at the ends of fragmented somatic chromosomes. To comprehensively investigate the telomeric regions of Mesorhabditis nematodes at the sequence level, we endeavored to collect a Mesorhabditis nematode in the Republic of Korea and acquire its highly contiguous genome sequences. We isolated a Mesorhabditis nematode and assembled its 108-Mb draft genome using both 6.3Gb (53 ×) of short-read and 3.0Gb (25 × , N50 = 5.7kb) of nanopore-based long-read sequencing data. Our genome assembly exhibits comparable quality to the public genome of Mesorhabditis belari in terms of contiguity and evolutionary conserved genes. Unexpectedly, our Mesorhabditis genome has many more interstitial telomeric sequences (ITSs), specifically subtelomeric ones, compared to the genomes of Caenorhabditis elegans and M. belari. Moreover, several subtelomeric sequences containing ITSs had 4-26 homologous sequences, implying they are highly repetitive. Based on this highly repetitive nature, we hypothesize that subtelomeric ITSs might have accumulated through the action of transposable elements containing ITSs. It still remains elusive whether these ITS-containing units are associated with programmed DNA elimination, but they may facilitate new telomere formation after DNA elimination. Our genomic resources for Mesorhabditis can aid in understanding how its distinct phenotypes have evolved.

Abf1 negatively regulates the expression of EPA1 and affects adhesion in Candida glabrata.

Introduction. Adherence is a major virulence trait in Candida glabrata that, in many strains, depends on the EPA (epithelial adhesin) genes, which confer the ability to adhere to epithelial and endothelial cells of the host. The EPA genes are generally found at subtelomeric regions, which makes them subject to subtelomeric silencing. In C. glabrata, subtelomeric silencing depends on different protein complexes, such as silent information regulator and yKu complexes, and other proteins, such as Repressor/activator protein 1 (Rap1) and Abf1. At the EPA1 locus, which encodes the main adhesin Epa1, we previously found at least two cis-acting elements, the protosilencer Sil2126 and the negative element, that contribute to the propagation of silencing from the telomere to the subtelomeric region.Hypothesis. Abf1 binds to the regulatory regions of EPA1 and other regions at the telomere E-R, thereby negatively regulating EPA1 transcription.Aim. To determine whether Abf1 and Rap1 silencing proteins bind to previously identified cis-acting elements on the right telomere of chromosome E (E-R subtelomeric region), resulting in negative regulation of EPA1 transcription and infer Abf1 and Rap1 recognition sites in C. glabrata.Methodology. We used chromatin immunoprecipitation (ChIP) followed by quantitative PCR to determine the binding sites for Abf1 and Rap1 in the intergenic regions between EPA1 and EPA2 and HYR1 and EPA1, and mutants were used to determine the silencing level of the EPA1 promoter region.Results. We found that Abf1 predominantly binds to the EPA1 promoter region, leading to negative regulation of EPA1 expression. Furthermore, the mutant abf1-43, which lacks the last 43 amino acids at its C-terminal end and is defective for subtelomeric silencing, exhibits hyperadherence to epithelial cells in vitro compared to the parental strain, suggesting that EPA1 is derepressed. We also determined the motif-binding sequences for Abf1 and Rap1 in C. glabrata using data from the ChIP assays.Conclusion. Together these data indicate that Abf1 negatively regulates EPA1 expression, leading to decreased adhesion of C. glabrata to epithelial cells.

Highly divergent satellitomes of two barley species of agronomic importance, Hordeum chilense and H. vulgare

In this paper, we have performed an in-depth study of the complete set of the satellite DNA (satDNA) families (i.e. the satellitomes) in the genome of two barley species of agronomic value in a breeding framework, H. chilense (H1 and H7 accessions) and H. vulgare (H106 accession), which can be useful tools for studying chromosome associations during meiosis. The study has led to the analysis of a total of 18 satDNA families in H. vulgare, 25 satDNA families in H. chilense (accession H1) and 27 satDNA families in H. chilense (accession H7) that constitute 46 different satDNA families forming 36 homology groups. Our study highlights different important contributions of evolutionary and applied interests. Thus, both barley species show very divergent satDNA profiles, which could be partly explained by the differential effects of domestication versus wildlife. Divergence derives from the differential amplification of different common ancestral satellites and the emergence of new satellites in H. chilense, usually from pre-existing ones but also random sequences. There are also differences between the two H. chilense accessions, which support genetically distinct groups. The fluorescence in situ hybridization (FISH) patterns of some satDNAs yield distinctive genetic markers for the identification of specific H. chilense or H. vulgare chromosomes. Some of the satellites have peculiar structures or are related to transposable elements which provide information about their origin and expansion. Among these, we discuss the existence of different (peri)centromeric satellites that supply this region with some plasticity important for centromere evolution. These peri(centromeric) satDNAs and the set of subtelomeric satDNAs (a total of 38 different families) are analyzed in the framework of breeding as the high diversity found in the subtelomeric regions might support their putative implication in chromosome recognition and pairing during meiosis, a key point in the production of addition/substitution lines and hybrids.

Haplotype-resolved de novo assembly revealed unique characteristics of alternative lengthening of telomeres in mouse embryonic stem cells.

Telomeres protect chromosome ends from DNA damage responses, and their dysfunction triggers genomic alterations like chromosome fusion and rearrangement, which can lead to cellular death. Certain cells, including specific cancer cells, adopt alternative lengthening of telomere (ALT) to counteract dysfunctional telomeres and proliferate indefinitely. While telomere instability and ALT activity are likely major sources of genomic alteration, the patterns and consequences of such changes at the nucleotide level in ALT cells remainunexplored. Here we generated haplotype-resolved genome assemblies for type I ALT mouse embryonic stem cells, facilitated by highly accurateor ultra-long reads and Hi-C reads. High-quality genome revealed ALT-specific complex chromosome end structures and various genomic alterations including over 1000 structural variants (SVs). The unique sequence (mTALT) used as a template for type I ALT telomeres showed traces of being recruited into the genome, with mTALT being replicated with remarkably high accuracy. Subtelomeric regions exhibited distinct characteristics: resistance to the accumulation of SVs and small variants. We genotyped SVs at allele resolution, identifying genes (Rgs6, Dpf3and Tacc2) crucial for maintaining ALT telomere stability. Our genome assembly-based approach elucidated the unique characteristics of ALT genome, offering insights into the genome evolution of cells surviving telomere-derived crisis.

Linked regulation of genome integrity and senescence-associated inflammation by p53.

Genomic instability and inflammation are distinct hallmarks of aging, but the connection between them is poorly understood. Understanding their interrelationship will help unravel new mechanisms and therapeutic targets of aging and age-associated diseases. Here we report a novel mechanism directly linking genomic instability and inflammation in senescent cells through a mitochondria-regulated molecular circuit driven by p53 and cytoplasmic chromatin fragments (CCF). We show, through activation or inactivation of p53 by genetic and pharmacologic approaches, that p53 suppresses CCF accumulation and the downstream inflammatory senescence-associated secretory phenotype (SASP), without affecting cell cycle arrest. p53 activation suppressed CCF formation by promoting DNA repair, and this is reflected in maintenance of genomic integrity, particularly in subtelomeric regions, as shown by single cell genome resequencing. Activation of p53 in aged mice by pharmacological inhibition of MDM2 reversed signatures of aging, including age- and senescence-associated transcriptomic signatures of inflammation and age-associated accumulation of monocytes and macrophages in liver. Remarkably, mitochondria in senescent cells suppressed p53 activity by promoting CCF formation and thereby restricting ATM-dependent nuclear DNA damage signaling. These data provide evidence for a mitochondria-regulated p53 signaling circuit in senescent cells that controls DNA repair, genome integrity, and senescence- and age-associated inflammation. This pathway is immunomodulatory in mice and a potential target for healthy aging interventions by small molecules already shown to activate p53.

Novel centromeric repetitive DNA elements reveal karyotype dynamics in polyploid sainfoin (Onobrychis viciifolia).

Polyploidy is a common feature in eukaryotes with one of paramount consequences leading to better environmental adaptation. Heterochromatin is often located at telomeres and centromeres and contains repetitive DNA sequences. Sainfoin (Onobrychis viciifolia) is an important perennial forage legume for sustainable agriculture. However, there are only a few studies on the sainfoin genome and chromosomes. In this study, novel tandem repetitive DNA sequences of the sainfoin genome (OnVi180, OnVi169, OnVi176 and OnVidimer) were characterized using bioinformatics, molecular and cytogenetic approaches. The OnVi180 and OnVi169 elements colocalized within functional centromeres. The OnVi176 and OnVidimer elements were localized in centromeric, subtelomeric and interstitial regions. We constructed a sainfoin karyotype that distinguishes the seven basic chromosome groups. Our study provides the first detailed description of heterochromatin and chromosome structure of sainfoin and proposes an origin of heterozygous ancestral genomes, possibly from the same ancestral diploid species, not necessarily from different species, or for chromosome rearrangements after polyploidy. Overall, we discuss our novel and complementary findings in a polyploid crop with unique and complex chromosomal features.

The highly dynamic satellitomes of cultivated wheat species.

Durum wheat, Triticum turgidum, and bread wheat, Triticum aestivum, are two allopolyploid species of very recent origin that have been subjected to intense selection programs during the thousands of years they have been cultivated. In this paper, we study the durum wheat satellitome and establish a comparative analysis with the previously published bread wheat satellitome. We revealed the durum wheat satellitome using the satMiner protocol which is based on consecutive rounds of clustering of Illumina reads by RepeatExplorer2, and estimated abundance and variation for each identified satDNA with RepeatMasker v4.0.5. We have also performed a deep satDNA families characterization including chromosomal location by Fluorescence In Situ Hybridization (FISH) in durum wheat and its comparison with FISH patterns in bread wheat. Basic Local Alignment Search Tool (BLAST®) was used for trailing each satDNA in the assembly of durum wheat genome through NCBI's Genome Data Viewer (GDW) and the genome assemblies of both species were compared. Sequence divergence and consensus turnover rate (CTR) between homologous satDNA families of durum and bread wheat were estimated using MEGA11. This study reveals that in an exceedingly short period, significant qualitative and quantitative changes have occurred in the set of satellite DNAs (satDNAs) of both species, with expansions/contractions of the number of repeats and the loci per satellite, different in each species, and a high rate of sequence change for most of these satellites, in addition to the emergence/loss of satDNAs not shared between the two species analysed. These evolutionary changes in satDNA are common between species but what is truly remarkable and novel about this study is that these processes have taken place in less than the last ~8000 years separating the two species, indicating an accelerated evolution of their satDNAs. These results, together with the relationship of many of these satellites with transposable elements and the polymorphisms they generate at the level of centromeres and subtelomeric regions of their chromosomes, are analysed and discussed in the context of the evolutionary origin of these species and the selection pressure exerted by man throughout the history of their cultivation.

PfMORC protein regulates chromatin accessibility and transcriptional repression in the human malaria parasite, Plasmodium falciparum.

The environmental challenges the human malaria parasite, Plasmodium falciparum, faces during its progression into its various lifecycle stages warrant the use of effective and highly regulated access to chromatin for transcriptional regulation. Microrchidia (MORC) proteins have been implicated in DNA compaction and gene silencing across plant and animal kingdoms. Accumulating evidence has shed light into the role MORC protein plays as a transcriptional switch in apicomplexan parasites. In this study, using CRISPR/Cas9 genome editing tool along with complementary molecular and genomics approaches, we demonstrate that PfMORC not only modulates chromatin structure and heterochromatin formation throughout the parasite erythrocytic cycle, but is also essential to the parasite survival. Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) experiments suggest that PfMORC binds to not only sub-telomeric regions and genes involved in antigenic variation but may also play a role in modulating stage transition. Protein knockdown experiments followed by chromatin conformation capture (Hi-C) studies indicate that downregulation of PfMORC impairs key histone marks and induces the collapse of the parasite heterochromatin structure leading to its death. All together these findings confirm that PfMORC plays a crucial role in chromatin structure and gene regulation, validating this factor as a strong candidate for novel antimalarial strategies.

Transposon-sequencing (Tn-seq) of the Candida glabrata reference strain CBS138 reveals epigenetic plasticity, structural variation, and intrinsic mechanisms of resistance to micafungin.

Candida glabrata (also called Nakaseomyces glabratus) is an opportunistic pathogen that can resist common antifungals and rapidly acquire multidrug resistance. A large amount of genetic variation exists between isolates, which complicates generalizations. Portable transposon-sequencing (Tn-seq) methods can efficiently provide genome-wide information on strain differences and genetic mechanisms. Using the Hermes transposon, the CBS138 reference strain and a commonly studied derivative termed 2001 were subjected to Tn-seq in control conditions and after exposure to varying doses of the clinical antifungal micafungin. The approach revealed large differences between these strains, including a 131-kb tandem duplication and a variety of fitness differences. Additionally, both strains exhibited up to 1,000-fold increased transposon accessibility in subtelomeric regions relative to the BG2 strain, indicative of open subtelomeric chromatin in these isolates and large epigenetic variation within the species. Unexpectedly, the Pdr1 transcription factor conferred resistance to micafungin through targets other than CDR1. Other micafungin resistance pathways were also revealed including mannosyltransferase activity and biosynthesis of the lipid precursor sphingosine, the inhibition of which by SDZ 90-215 and myriocin enhanced the potency of micafungin in vitro. These findings provide insights into the complexity of the C. glabrata species as well as strategies for improving antifungal efficacy.

Male-biased recombination at chromosome ends in a songbird revealed by precisely mapping crossover positions.

Recombination plays a crucial role in evolution by generating novel haplotypes and disrupting linkage between genes, thereby enhancing the efficiency of selection. Here, we analyze the genomes of 12 great reed warblers (Acrocephalus arundinaceus) in a 3-generation pedigree to identify precise crossover positions along the chromosomes. We located more than 200 crossovers and found that these were highly concentrated toward the telomeric ends of the chromosomes. Apart from this major pattern in the recombination landscape, we found significantly higher frequencies of crossovers in genic compared with intergenic regions, and in exons compared with introns. Moreover, while the number of recombination events was similar between the sexes, the crossovers were located significantly closer to the ends of paternal compared with maternal chromosomes. In conclusion, our study of the great reed warbler revealed substantial variation in crossover frequencies within chromosomes, with a distinct bias toward the sub-telomeric regions, particularly on the paternal side. These findings emphasize the importance of thoroughly screening the entire length of chromosomes to characterize the recombination landscape and uncover potential sex-biases in recombination.

French National Protocol for diagnosis and care of facioscapulohumeral muscular dystrophy (FSHD).

Facioscapulohumeral muscular dystrophy (FSHD) is one of the most common genetically inherited myopathies in adults. It is characterized by incomplete penetrance and variable expressivity. Typically, FSHD patients display asymmetric weakness of facial, scapular, and humeral muscles that may progress to other muscle groups, particularly the abdominal and lower limb muscles. Early-onset patients display more severe muscle weakness and atrophy, resulting in a higher frequency of associated skeletal abnormalities. In these patients, multisystem involvement, including respiratory, ocular, and auditory, is more frequent and severe and may include the central nervous system. Adult-onset FSHD patients may also display some degree of multisystem involvement which mainly remains subclinical. In 95% of cases, FSHD patients carry a pathogenic contraction of the D4Z4 repeat units (RUs) in the subtelomeric region of chromosome 4 (4q35), which leads to the expression of DUX4 retrogene, toxic for muscles (FSHD1). Five percent of patients display the same clinical phenotype in association with a mutation in the SMCHD1 gene located in chromosome 18, inducing epigenetic modifications of the 4q D4Z4 repeated region and expression of DUX4 retrogene. This review highlights the complexities and challenges of diagnosing and managing FSHD, underscoring the importance of standardized approaches for optimal patient outcomes. It emphasizes the critical role of multidisciplinary care in addressing the diverse manifestations of FSHD across different age groups, from skeletal abnormalities in early-onset cases to the often-subclinical multisystem involvement in adults. With no current cure, the focus on alleviating symptoms and slowing disease progression through coordinated care is paramount.

Cell-type-specific effects of age and sex on human cortical neurons

Altered transcriptional and epigenetic regulation of brain cell types may contribute to cognitive changes withadvanced age. Using single-nucleus multi-omic DNA methylation and transcriptome sequencing (snmCT-seq) in frontal cortex from young adult and aged donors, we found widespread age- and sex-related variation in specific neuron types. The proportion of inhibitory SST- and VIP-expressing neurons was reduced in aged donors. Excitatory neurons had more profound age-related changes in their gene expression and DNA methylation than inhibitory cells. Hundreds of genes involved in synaptic activity, including EGR1, were less expressed in aged adults. Genes located in subtelomeric regions increased their expression with age and correlated with reduced telomere length. We further mapped cell-type-specific sex differences in gene expression and X-inactivation escape genes. Multi-omic single-nucleus epigenomes and transcriptomes provide new insight into the effects of age and sex on human neurons.

Complete telomere-to-telomere genomes uncover virulence evolution conferred by chromosome fusion in oomycete plant pathogens

Variations in chromosome number are occasionally observed among oomycetes, a group that includes many plant pathogens, but the emergence of such variations and their effects on genome and virulence evolution remain ambiguous. We generated complete telomere-to-telomere genome assemblies for Phytophthora sojae, Globisporangium ultimum, Pythium oligandrum, and G. spinosum. Reconstructing the karyotype of the most recent common ancestor in Peronosporales revealed that frequent chromosome fusion and fission drove changes in chromosome number. Centromeres enriched with Copia-like transposons may contribute to chromosome fusion and fission events. Chromosome fusion facilitated the emergence of pathogenicity genes and their adaptive evolution. Effectors tended to duplicate in the sub-telomere regions of fused chromosomes, which exhibited evolutionary features distinct to the non-fused chromosomes. By integrating ancestral genomic dynamics and structural predictions, we have identified secreted Ankyrin repeat-containing proteins (ANKs) as a novel class of effectors in P. sojae. Phylogenetic analysis and experiments further revealed that ANK is a specifically expanded effector family in oomycetes. These results revealed chromosome dynamics in oomycete plant pathogens, and provided novel insights into karyotype and effector evolution.

Exchange of subtelomeric regions between chromosomes 4q and 10q reverts the FSHD genotype and phenotype.

The most common form of facioscapulohumeral dystrophy (FSHD1) is caused by a partial loss of the D4Z4 macrosatellite repeat array in the subtelomeric region of chromosome 4. Patients with FSHD1 typically carry 1 to 10 D4Z4 repeats, whereas nonaffected individuals have 11 to 150 repeats. The ~150-kilobyte subtelomeric region of the chromosome 10q exhibits a ~99% sequence identity to the 4q, including the D4Z4 array. Nevertheless, contractions of the chr10 array do not cause FSHD or any known disease, as in most people D4Z4 array on chr10 is flanked by the nonfunctional polyadenylation signal, not permitting the DUX4 expression. Here, we attempted to correct the FSHD genotype by a CRISPR-Cas9-induced exchange of the chr4 and chr10 subtelomeric regions. We demonstrated that the induced t(4;10) translocation can generate recombinant genotypes translated into improved FSHD phenotype. FSHD myoblasts with the t(4;10) exhibited reduced expression of the DUX4 targets, restored PAX7 target expression, reduced sensitivity to oxidative stress, and improved differentiation capacity.

Annotations of four high-quality indigenous chicken genomes identify more than one thousand missing genes in subtelomeric regions and micro-chromosomes with high G/C contents.

Although multiple chicken genomes have been assembled and annotated, the numbers of protein-coding genes in chicken genomes and their variation among breeds are still uncertain due to the low quality of these genome assemblies and limited resources used in their gene annotations. To fill these gaps, we recently assembled genomes of four indigenous chicken breeds with distinct traits at chromosome-level. In this study, we annotated genes in each of these assembled genomes using a combination of RNA-seq- and homology-based approaches. We identified varying numbers (17,497-17,718) of protein-coding genes in the four indigenous chicken genomes, while recovering 51 of the 274 "missing" genes in birds in general, and 36 of the 174 "missing" genes in chickens in particular. Intriguingly, based on deeply sequenced RNA-seq data collected in multiple tissues in the four breeds, we found 571 ~ 627 protein-coding genes in each genome, which were missing in the annotations of the reference chicken genomes (GRCg6a and GRCg7b/w). After removing redundancy, we ended up with a total of 1,420 newly annotated genes (NAGs). The NAGs tend to be found in subtelomeric regions of macro-chromosomes (chr1 to chr5, plus chrZ) and middle chromosomes (chr6 to chr13, plus chrW), as well as in micro-chromosomes (chr14 to chr39) and unplaced contigs, where G/C contents are high. Moreover, the NAGs have elevated quadruplexes G frequencies, while both G/C contents and quadruplexes G frequencies in their surrounding regions are also high. The NAGs showed tissue-specific expression, and we were able to verify 39 (92.9%) of 42 randomly selected ones in various tissues of the four chicken breeds using RT-qPCR experiments. Most of the NAGs were also encoded in the reference chicken genomes, thus, these genomes might harbor more genes than previously thought. The NAGs are widely distributed in wild, indigenous and commercial chickens, and they might play critical roles in chicken physiology. Counting these new genes, chicken genomes harbor more genes than originally thought.

Genome assembly of a symbiotic balantidia (Balantidium ctenopharyngodoni) in fish hindgut

Balantidium ctenopharyngodoni is identified as the sole ciliate species that exclusively resides within the hindgut of grass carp with high prevalence and intensity. In this study, the successful cultivation of B. ctenopharyngodoni enabled us to collect enough cells for genome sequencing. Consequently, we acquired a high-quality genome assembly spanning 68.66 Mb, encompassing a total of 22,334 nanochromosomes. Furthermore, we predicted 29,348 protein-coding genes, and 95.5% of them was supported by the RNA-seq data. The trend of GC content in the subtelomeric regions of single-gene chromosomes was similar to other ciliates containing nanochromosomes. A large number of genes encoding carbohydrate-binding modules with affinities for starch and peptidoglycans was identified. The identification of mitochondrion-related organelles (MROs) within genome indicates its well-suited adaptation to the anaerobic conditions in the hindgut environment. In summary, our results will offer resources for understanding the genetic basis and molecular adaptations of balantidia to hindgut of herbivorous fish.

Cytogenomic Characterization of Transposable Elements and Satellite DNA in Passiflora L. Species.

The species Passiflora alata, P. cincinnata, and P. edulis have great economic value due to the use of their fruits for human consumption. In this study, we compared the repetitive genome fractions of these three species. The compositions of the repetitive DNA of these three species' genomes were analyzed using clustering and identification of the repetitive sequences with RepeatExplorer. It was found that repetitive DNA content represents 74.70%, 66.86%, and 62.24% of the genome of P. alata, P. edulis, and P. cincinnata, respectively. LTR Ty3/Gypsy retrotransposons represent the highest genome proportions in P. alata and P. edulis, while Ty1/Copia comprises the largest proportion of P. cincinnata genome. Chromosomal mapping by Fluorescent In Situ Hybridization (FISH) showed that LTR retrotransposons have a dispersed distribution along chromosomes. The subtelomeric region of chromosomes is where 145 bp satellite DNA is located, suggesting that these elements may play important roles in genome structure and organization in these species. In this work, we obtained the first global characterization of the composition of repetitive DNA in Passiflora, showing that an increase in genome size is related to an increase in repetitive DNA, which represents an important evolutionary route for these species.