Subtelomeric Regions Research Articles

Characterization of Host-Specific Genes from Pine- and Grass-Associated Species of the Fusarium fujikuroi Species Complex.

The Fusarium fujikuroi species complex (FFSC) includes socioeconomically important pathogens that cause disease for numerous crops and synthesize a variety of secondary metabolites that can contaminate feedstocks and food. Here, we used comparative genomics to elucidate processes underlying the ability of pine-associated and grass-associated FFSC species to colonize tissues of their respective plant hosts. We characterized the identity, possible functions, evolutionary origins, and chromosomal positions of the host-range-associated genes encoded by the two groups of fungi. The 72 and 47 genes identified as unique to the respective genome groups were potentially involved in diverse processes, ranging from transcription, regulation, and substrate transport through to virulence/pathogenicity. Most genes arose early during the evolution of Fusarium/FFSC and were only subsequently retained in some lineages, while some had origins outside Fusarium. Although differences in the densities of these genes were especially noticeable on the conditionally dispensable chromosome of F. temperatum (representing the grass-associates) and F. circinatum (representing the pine-associates), the host-range-associated genes tended to be located towards the subtelomeric regions of chromosomes. Taken together, these results demonstrate that multiple mechanisms drive the emergence of genes in the grass- and pine-associated FFSC taxa examined. It also highlighted the diversity of the molecular processes potentially underlying niche-specificity in these and other Fusarium species.

Characterization of a Wheat-Dasypyrum breviaristatum Chromosome Addition and Its Derived Progenies Carrying Novel Dasypyrum-Specific Gliadin Genes

The construction of the 28-chromosome karyotype of Dasypyrum breviaristatum was undertaken using multicolor non-denaturing fluorescent in situ hybridization (ND-FISH) and Oligo-FISH painting protocols. A novel wheat-D. breviaristatum line D2138 contained 44 chromosomes including a pair of D. breviaristatum 6VbS.2VbL translocation chromosomes. Individual F2 and F3 progenies of a cross between D2138 with wheat lines CM62, MY11 and JM22, respectively, were characterized using ND-FISH and molecular markers. A relatively high chromosome alteration rate within wheat and D. breviaristatum 6VbS and 2VbL was observed in the three progeny populations, suggesting that chromosome 6VbS.2VbL has a gametocidal-like gene. The different types of translocation and deletion lines allowed localization of D. breviaristatum-specific gliadin coding genes on sub-telomeric regions of 6VbS by PCR and acid polyacrylamide gel electrophoresis analysis. The positive effect of the D. breviaristatum 6VbS on agronomic and quality characters was also demonstrated. The new wheat-D. breviaristatum 6VbS and 2VbL translocation lines will be useful as novel germplasm for breeding purposes.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton.

SummaryInterspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056–1255 ENAs showed transgressive expression in fibres at 5–15 dpa (days post‐anthesis) of some BILs, 520 of which were located in FL‐quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL‐QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR‐cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Identification and Mapping of bs8, a Novel Locus Conferring Resistance to Bacterial Spot Caused by Xanthomonas gardneri.

Although cultivars possessing recessive resistance alleles provide effective control of bacterial spot of pepper (Capsicum annuum), the deployed resistance gene, bs5, is ineffective against Xanthomonas gardneri, one of the pathogenic species. Resistance against X. gardneri was identified in C. annuum accession PI 163192, and this study sought to characterize this novel resistance and to map the resistance gene(s) to the pepper genome. We crossed PI 163192 with the susceptible cultivar Early Calwonder (ECW) to develop resistant near-isogenic lines (NILs) of ECW, designated ECW80R. The novel resistance in ECW80R was determined to be quantitative, recessively inherited, and non-hypersensitive-response causing, and inhibits lesion expansion and chlorosis. Presence of the resistance in NILs decreased the in planta bacterial population by ninefold compared with ECW. Bulked segregant analysis of resistant and susceptible individuals from an F2 population using whole genome single nucleotide polymorphisms identified a major resistance locus within an approximate 6-Mbp interval on the subtelomeric region of chromosome 11. We developed markers spanning this region and used these to genotype backcross F2 populations, which further delimited the resistance locus within a 2.3-Mbp interval. The novel resistance locus has been designated bs8. ECW80R and the linked markers developed in this study should prove useful for breeders seeking to advance this resistance into commercially relevant germplasm and for pyramiding bs8 with other resistance alleles such as bs5 and bs6. The allele bs8 will help prolong the durability of bacterial spot resistance in pepper and improve resistance to multiple species of Xanthomonas.

Genomic Diversity across Candida auris Clinical Isolates Shapes Rapid Development of Antifungal Resistance In Vitro and In Vivo.

ABSTRACTAntifungal drug resistance and tolerance pose a serious threat to global public health. In the human fungal pathogen, Candida auris, resistance to triazole, polyene, and echinocandin antifungals is rising, resulting in multidrug resistant isolates. Here, we use genome analysis and in vitro evolution of 17 new clinical isolates of C. auris from clades I and IV to determine how quickly resistance mutations arise, the stability of resistance in the absence of drug, and the impact of genetic background on evolutionary trajectories. We evolved each isolate in the absence of drug as well as in low and high concentrations of fluconazole. In just three passages, we observed genomic and phenotypic changes including karyotype alterations, aneuploidy, acquisition of point mutations, and increases in MIC values within the populations. Fluconazole resistance was stable in the absence of drug, indicating little to no fitness cost associated with resistance. Importantly, two isolates substantially increased resistance to ≥256 μg/mL fluconazole. Multiple evolutionary pathways and mutations associated with increased fluconazole resistance occurred simultaneously within the same population. Strikingly, the subtelomeric regions of C. auris were highly dynamic as deletion of multiple genes near the subtelomeres occurred during the three passages in several populations. Finally, we discovered a mutator phenotype in a clinical isolate of C. auris. This isolate had elevated mutation rates compared to other isolates and acquired substantial resistance during evolution in vitro and in vivo supporting that the genetic background of clinical isolates can have a significant effect on evolutionary potential.

Intra-Species Genomic Variation in the Pine Pathogen Fusarium circinatum.

Fusarium circinatum is an important global pathogen of pine trees. Genome plasticity has been observed in different isolates of the fungus, but no genome comparisons are available. To address this gap, we sequenced and assembled to chromosome level five isolates of F. circinatum. These genomes were analysed together with previously published genomes of F. circinatum isolates, FSP34 and KS17. Multi-sample variant calling identified a total of 461,683 micro variants (SNPs and small indels) and a total of 1828 macro structural variants of which 1717 were copy number variants and 111 were inversions. The variant density was higher on the sub-telomeric regions of chromosomes. Variant annotation revealed that genes involved in transcription, transport, metabolism and transmembrane proteins were overrepresented in gene sets that were affected by high impact variants. A core genome representing genomic elements that were conserved in all the isolates and a non-redundant pangenome representing all genomic elements is presented. Whole genome alignments showed that an average of 93% of the genomic elements were present in all isolates. The results of this study reveal that some genomic elements are not conserved within the isolates and some variants are high impact. The described genome-scale variations will help to inform novel disease management strategies against the pathogen.

Chromatin structure can introduce systematic biases in genome-wide analyses of Plasmodium falciparum

Background: The maintenance, regulation, and dynamics of heterochromatin in the human malaria parasite, Plasmodium falciparum, has drawn increasing attention due to its regulatory role in mutually exclusive virulence gene expression and the silencing of key developmental regulators. The advent of genome-wide analyses such as chromatin-immunoprecipitation followed by sequencing (ChIP-seq) has been instrumental in understanding chromatin composition; however, even in model organisms, ChIP-seq experiments are susceptible to intrinsic experimental biases arising from underlying chromatin structure. Methods: We performed a control ChIP-seq experiment, re-analyzed previously published ChIP-seq datasets and compared different analysis approaches to characterize biases of genome-wide analyses in P. falciparum. Results: We found that heterochromatic regions in input control samples used for ChIP-seq normalization are systematically underrepresented in regard to sequencing coverage across the P. falciparum genome. This underrepresentation, in combination with a non-specific or inefficient immunoprecipitation, can lead to the identification of false enrichment and peaks across these regions. We observed that such biases can also be seen at background levels in specific and efficient ChIP-seq experiments. We further report on how different read mapping approaches can also skew sequencing coverage within highly similar subtelomeric regions and virulence gene families. To ameliorate these issues, we discuss orthogonal methods that can be used to characterize bona fide chromatin-associated proteins. Conclusions: Our results highlight the impact of chromatin structure on genome-wide analyses in the parasite and the need for caution when characterizing chromatin-associated proteins and features.

Recombination landscape divergence between populations is marked by larger low-recombining regions in domesticated rye.

The genomic landscape of recombination plays an essential role in evolution. Patterns of recombination are highly variable along chromosomes, between sexes, individuals, populations, and species. In many eukaryotes, recombination rates are elevated in sub-telomeric regions and drastically reduced near centromeres, resulting in large low-recombining (LR) regions. The processes of recombination are influenced by genetic factors, such as different alleles of genes involved in meiosis and chromatin structure, as well as external environmental stimuli like temperature and overall stress. In this work, we focused on the genomic landscapes of recombination in a collection of 916 rye (Secale cereale) individuals. By analyzing population structure among individuals of different domestication status and geographic origin, we detected high levels of admixture, reflecting the reproductive biology of a self-incompatible, wind-pollinating grass species. We then analyzed patterns of recombination in overlapping subpopulations, which revealed substantial variation in the physical size of LR regions, with a tendency for larger LR regions in domesticated subpopulations. Genome-wide association scans (GWAS) for LR region size revealed a major quantitative-trait-locus (QTL) at which, among 18 annotated genes, an ortholog of histone H4 acetyltransferase ESA1 was located. Rye individuals belonging to domesticated subpopulations showed increased synaptonemal complex length, but no difference in crossover frequency, indicating that only the recombination landscape is different. Furthermore, the genomic region harboring rye ScESA1 showed moderate patterns of selection in domesticated subpopulations, suggesting that larger LR regions were indirectly selected during domestication to achieve more homogeneous populations for agricultural use.

Plasmodium falciparum PfRUVBL proteins bind at the TARE region and var gene promoter located in the subtelomeric region.

In order to survive and establish infection, the Plasmodium parasite employs various strategies to evade the host immune response. The var genes family, a repertoire of 60 genes, expresses parasite-specific protein PfEMP1, a variable surface antigen, on the membrane of infected erythrocytes, and by continuously switching the variants of PfEMP1, help the parasite to avoid detection and destruction by the host immune system during the intra-erythrocytic developmental cycle. Although chromatin modifications are recognised to be a prominent phenomenon in regulation of mono-allelic expression of these var genes, the precise histone codes and molecular players and mechanisms guiding these modifications have yet to be unravelled in depth. In this study, we have functionally characterised RUVBL proteins of Plasmodium falciparum and shown that PfMYST (an essential lysine acetyl transferase) and PfRUVBL protein complex occupy the TARE region and var gene promoter in the ring stage of the parasite. Further, we have demonstrated that the PfMYST/PfRUVBL complex interacts with core histones, H3 and H4. Overall the findings of this study add further information by identifying the potential role of epigenetic regulators, PfMYST and PfRUVBL, in the regulation of monoallelic expression of var genes in the malaria parasite.

The nature and genomic landscape of repetitive DNA classes in Chrysanthemum nankingense shows recent genomic changes.

Tandemly repeated DNA and transposable elements represent most of the DNA in higher plant genomes. High-throughput sequencing allows a survey of the DNA in a genome, but whole-genome assembly can miss a substantial fraction of highly repeated sequence motifs. Chrysanthemum nankingense (2n = 2x = 18; genome size = 3.07 Gb; Asteraceae), a diploid reference for the many auto- and allopolyploids in the genus, was considered as an ancestral species and serves as an ornamental plant and high-value food. We aimed to characterize the major repetitive DNA motifs, understand their structure and identify key features that are shaped by genome and sequence evolution. Graph-based clustering with RepeatExplorer was used to identify and classify repetitive motifs in 2.14 millions of 250-bp paired-end Illumina reads from total genomic DNA of C. nankingense. Independently, the frequency of all canonical motifs k-bases long was counted in the raw read data and abundant k-mers (16, 21, 32, 64 and 128) were extracted and assembled to generate longer contigs for repetitive motif identification. For comparison, long terminal repeat retrotransposons were checked in the published C. nankingense reference genome. Fluorescent in situ hybridization was performed to show the chromosomal distribution of the main types of repetitive motifs. Apart from rDNA (0.86 % of the total genome), a few microsatellites (0.16 %), and telomeric sequences, no highly abundant tandem repeats were identified. There were many transposable elements: 40 % of the genome had sequences with recognizable domains related to transposable elements. Long terminal repeat retrotransposons showed widespread distribution over chromosomes, although different sequence families had characteristic features such as abundance at or exclusion from centromeric or subtelomeric regions. Another group of very abundant repetitive motifs, including those most identified as low-complexity sequences (9.07 %) in the genome, showed no similarity to known sequence motifs or tandemly repeated elements. The Chrysanthemum genome has an unusual structure with a very low proportion of tandemly repeated sequences (~1.02 %) in the genome, and a high proportion of low-complexity sequences, most likely degenerated remains of transposable elements. Identifying the presence, nature and genomic organization of major genome fractions enables inference of the evolutionary history of sequences, including degeneration and loss, critical to understanding biodiversity and diversification processes in the genomes of diploid and polyploid Chrysanthemum, Asteraceae and plants more widely.

De Novo Assembly of Plasmodium knowlesi Genomes From Clinical Samples Explains the Counterintuitive Intrachromosomal Organization of Variant SICAvar and kir Multiple Gene Family Members.

Plasmodium knowlesi, a malaria parasite of Old World macaque monkeys, is used extensively to model Plasmodium biology. Recently, P. knowlesi was found in the human population of Southeast Asia, particularly Malaysia. P. knowlesi causes uncomplicated to severe and fatal malaria in the human host with features in common with the more prevalent and virulent malaria caused by Plasmodium falciparum. As such, P. knowlesi presents a unique opportunity to develop experimental translational model systems for malaria pathophysiology informed by clinical data from same-species human infections. Experimental lines of P. knowlesi represent well-characterized genetically stable parasites, and to maximize their utility as a backdrop for understanding malaria pathophysiology, genetically diverse contemporary clinical isolates, essentially wild-type, require comparable characterization. The Oxford Nanopore PCR-free long-read sequencing platform was used to sequence and de novo assemble P. knowlesi genomes from frozen clinical samples. The sequencing platform and assembly pipelines were designed to facilitate capturing data and describing, for the first time, P. knowlesi schizont-infected cell agglutination (SICA) var and Knowlesi-Interspersed Repeats (kir) multiple gene families in parasites acquired from nature. The SICAvar gene family members code for antigenically variant proteins analogous to the virulence-associated P. falciparum erythrocyte membrane protein (PfEMP1) multiple var gene family. Evidence presented here suggests that the SICAvar family members have arisen through a process of gene duplication, selection pressure, and variation. Highly evolving genes including PfEMP1family members tend to be restricted to relatively unstable sub-telomeric regions that drive change with core genes protected in genetically stable intrachromosomal locations. The comparable SICAvar and kir gene family members are counter-intuitively located across chromosomes. Here, we demonstrate that, in contrast to conserved core genes, SICAvar and kir genes occupy otherwise gene-sparse chromosomal locations that accommodate rapid evolution and change. The novel methods presented here offer the malaria research community not only new tools to generate comprehensive genome sequence data from small clinical samples but also new insight into the complexity of clinically important real-world parasites.

The Shape of Silence: A Role for the Nuclear Envelope in Transcriptional Regulation

The lipid drug edelfosine integrates into the nuclear envelope (NE) of S. cerevisiae (budding yeast) and results in its morphological deformation. As biologists we are guided by a central tenant that form and function are inextricably linked. This is true for the thumb on a hand and a protein produced in a cell. We hypothesized that this will also be true for the organelles inside a eukaryotic cell, wherein changes in nuclear shape would lead to changes in nuclear functions, like alterations in gene expression and protein localization. Using live cell imaging combined with yeast genetics and biochemical approaches, we found that edelfosine treatment led to dispersion of the membrane‐associated SIR complex and triggered processing of the membrane‐sensing transcription factors, Mga2 and Spt23. By a complementary approach of RNA‐seq followed by transcriptomic analysis, we identified that targets of Mga2 as well as subtelomeric regions were upregulated in response to edelfosine, while ribosomal protein targets of Rap1 were downregulated. Our work suggests that disruption of the nuclear membrane is sufficient to trigger changes in membrane‐associated transcription factors and chromatin remodellers. Our data indicate that nuclear membrane integrity is linked to transcriptional regulation. The NE could be a novel target for chemotherapeutics as lipid drugs like edelfosine do not cause DNA damage and would likely not pose mutagenic threats.

Transcription of ncRNAs promotes repair of UV induced DNA lesions in Saccharomyces cerevisiae subtelomeres.

Ultraviolet light causes DNA lesions that are removed by nucleotide excision repair (NER). The efficiency of NER is conditional to transcription and chromatin structure. UV induced photoproducts are repaired faster in the gene transcribed strands than in the non-transcribed strands or in transcriptionally inactive regions of the genome. This specificity of NER is known as transcription-coupled repair (TCR). The discovery of pervasive non-coding RNA transcription (ncRNA) advocates for ubiquitous contribution of TCR to the repair of UV photoproducts, beyond the repair of active gene-transcribed strands. Chromatin rules transcription, and telomeres form a complex structure of proteins that silences nearby engineered ectopic genes. The essential protective function of telomeres also includes preventing unwanted repair of double-strand breaks. Thus, telomeres were thought to be transcriptionally inert, but more recently, ncRNA transcription was found to initiate in subtelomeric regions. On the other hand, induced DNA lesions like the UV photoproducts must be recognized and repaired also at the ends of chromosomes. In this study, repair of UV induced DNA lesions was analyzed in the subtelomeric regions of budding yeast. The T4-endonuclease V nicking-activity at cyclobutene pyrimidine dimer (CPD) sites was exploited to monitor CPD formation and repair. The presence of two photoproducts, CPDs and pyrimidine (6,4)-pyrimidones (6-4PPs), was verified by the effective and precise blockage of Taq DNA polymerase at these sites. The results indicate that UV photoproducts in silenced heterochromatin are slowly repaired, but that ncRNA transcription enhances NER throughout one subtelomeric element, called Y', and in distinct short segments of the second, more conserved element, called X. Therefore, ncRNA-transcription dependent TCR assists global genome repair to remove CPDs and 6-4PPs from subtelomeric DNA.

Gestational exposure to NO2 aggravates placental senescence

Placental senescence is a normal physiological process of placenta, while premature placental senescence has been confirmed to be associated with some adverse pregnancy complications. Epidemiological studies indicate that NO2 exposure can aggravate placental senescence which is represented by fibrosis and abnormal telomere homeostasis, etc. In this study, pregnant C57BL/6 mice were exposed to NO2 (2.5 ppm, 5 h/day) daily in a dynamic exposure chamber throughout the gestation period, and were sacrificed at embryonic day 13.5 (E13.5), E15.5 and E18.5. Placenta were harvested and conducted for histopathological examination and telomere evaluation. Our results showed that gestational NO2 exposure significantly aggravated placental fibrosis and calcification, and up-regulated the related bio-markers (connective tissue growth factor (Ctgf) and transforming growth factor-β1 (Tgf-β1)) at E18.5. In addition, gestational exposure to NO2 also activated senescence related pathway (p53/p21) at E18.5. Furthermore, gestational NO2 exposure significantly shortened telomere length at E18.5, and the expression of telomere homeostasis regulation genes telomeric repeat binding factor 1 (Trf1), protection of telomeres 1a (Pot1a) and Pot1b were significantly increased while telomerase reverse transcriptase (Tert) was suppressed after NO2 exposure at E13.5 or E18.5, respectively. Importantly, DNA methylation status of the 22nd at E13.5 and 32nd at E18.5 site in sub-telomeric region of chromosome 1 was significantly altered. Based on the above results, our present study indicated that gestational NO2 exposure could lead to premature placental senescence during the late trimester of pregnancy via aggravation of fibrosis and telomere length shortening regulated by telomere regulatory enzyme and DNA methylation.

Full-Length Genome of an Ogataea polymorpha Strain CBS4732 ura3Δ Reveals Large Duplicated Segments in Subtelomeric Regions

BackgroundCurrently, methylotrophic yeasts (e.g., Pichia pastoris, Ogataea polymorpha, and Candida boindii) are subjects of intense genomics studies in basic research and industrial applications. In the genus Ogataea, most research is focused on three basic O. polymorpha strains-CBS4732, NCYC495, and DL-1. However, the relationship between CBS4732, NCYC495, and DL-1 remains unclear, as the genomic differences between them have not be exactly determined without their high-quality complete genomes. As a nutritionally deficient mutant derived from CBS4732, the O. polymorpha strain CBS4732 ura3Δ (named HU-11) is being used for high-yield production of several important proteins or peptides. HU-11 has the same reference genome as CBS4732 (noted as HU-11/CBS4732), because the only genomic difference between them is a 5-bp insertion.ResultsIn the present study, we have assembled the full-length genome of O. polymorpha HU-11/CBS4732 using high-depth PacBio and Illumina data. Long terminal repeat retrotransposons (LTR-rts), rDNA, 5′ and 3′ telomeric, subtelomeric, low complexity and other repeat regions were exactly determined to improve the genome quality. In brief, the main findings include complete rDNAs, complete LTR-rts, three large duplicated segments in subtelomeric regions and three structural variations between the HU-11/CBS4732 and NCYC495 genomes. These findings are very important for the assembly of full-length genomes of yeast and the correction of assembly errors in the published genomes of Ogataea spp. HU-11/CBS4732 is so phylogenetically close to NCYC495 that the syntenic regions cover nearly 100% of their genomes. Moreover, HU-11/CBS4732 and NCYC495 share a nucleotide identity of 99.5% through their whole genomes. CBS4732 and NCYC495 can be regarded as the same strain in basic research and industrial applications.ConclusionThe present study preliminarily revealed the relationship between CBS4732, NCYC495, and DL-1. Our findings provide new opportunities for in-depth understanding of genome evolution in methylotrophic yeasts and lay the foundations for the industrial applications of O. polymorpha CBS4732, NCYC495, DL-1, and their derivative strains. The full-length genome of O. polymorpha HU-11/CBS4732 should be included into the NCBI RefSeq database for future studies of Ogataea spp.

De Novo Assembly of 20 Chicken Genomes Reveals the Undetectable Phenomenon for Thousands of Core Genes on Microchromosomes and Subtelomeric Regions.

The gene numbers and evolutionary rates of birds were assumed to be much lower than those of mammals, which is in sharp contrast to the huge species number and morphological diversity of birds. It is, therefore, necessary to construct a complete avian genome and analyze its evolution. We constructed a chicken pan-genome from 20 de novo assembled genomes with high sequencing depth, and identified 1,335 protein-coding genes and 3,011 long noncoding RNAs not found in GRCg6a. The majority of these novel genes were detected across most individuals of the examined transcriptomes but were seldomly measured in each of the DNA sequencing data regardless of Illumina or PacBio technology. Furthermore, different from previous pan-genome models, most of these novel genes were overrepresented on chromosomal subtelomeric regions and microchromosomes, surrounded by extremely high proportions of tandem repeats, which strongly blocks DNA sequencing. These hidden genes were proved to be shared by all chicken genomes, included many housekeeping genes, and enriched in immune pathways. Comparative genomics revealed the novel genes had 3-fold elevated substitution rates than known ones, updating the knowledge about evolutionary rates in birds. Our study provides a framework for constructing a better chicken genome, which will contribute toward the understanding of avian evolution and the improvement of poultry breeding.

Copy Number Variation Analysis of Euploid Pregnancy Loss.

Objectives: Copy number variant (CNV) is believed to be the potential genetic cause of pregnancy loss. However, CNVs less than 3 Mb in euploid products of conceptions (POCs) remain largely unexplored. The aim of this study was to investigate the features of CNVs less than 3 Mb in POCs and their potential clinical significance in pregnancy loss/fetal death. Methods: CNV data were extracted from a cohort in our institution and 19 peer-reviewed publications, and only those CNVs less than 3 Mb detected in euploid pregnancy loss/fetal death were included. We conducted a CNV map to analyze the distribution of CNVs in chromosomes using R packages karyoploteR_1.10.5. Gene names and annotated gene types covered by those CNVs were mined from the human Release 19 reference genome file and GENECODE database. We assessed the expression patterns and the consequences of murine knock-out of those genes using TiGER and Mouse Genome Informatics (MGI) databases. Functional enrichment and pathway analysis for genes in CNVs were performed using clusterProfiler V3.12.0. Result: Breakpoints of 564 CNVs less than 3 Mb were obtained from 442 euploid POCs, with 349 gains and 185 losses. The CNV map showed that CNVs were distributed in all chromosomes, with the highest frequency detected in chromosome 22 and the lowest frequency in chromosome Y, and CNVs showed a higher density in the pericentromeric and sub-telomeric regions. A total of 5,414 genes mined from the CNV regions (CNVRs), Gene Ontology (GO), and pathway analysis showed that the genes were significantly enriched in multiple terms, especially in sensory perception, membrane region, and tight junction. A total of 995 protein-coding genes have been reported to present mammalian phenotypes in MGI, and 276 of them lead to embryonic lethality or abnormal embryo/placenta in knock-out mouse models. CNV located at 19p13.3 was the most common CNV of all POCs. Conclusion: CNVs less than 3 Mb in euploid POCs distribute unevenly in all chromosomes, and a higher density was seen in the pericentromeric and sub-telomeric regions. The genes in those CNVRs are significantly enriched in biological processes and pathways that are important to embryonic/fetal development. CNV in 19p13.3 and the variations of ARID3A and FSTL3 might contribute to pregnancy loss.

A Conserved Hypothetical Gene Is Required but Not Sufficient for Ptr ToxC Production in Pyrenophora tritici-repentis.

The fungus Pyrenophora tritici-repentis causes tan spot, an important foliar disease of wheat worldwide. The fungal pathogen produces three necrotrophic effectors, namely Ptr ToxA, Ptr ToxB, and Ptr ToxC to induce necrosis or chlorosis in wheat. Both Ptr ToxA and Ptr ToxB are proteins, and their encoding genes have been cloned. Ptr ToxC was characterized as a low-molecular weight molecule 20 years ago but the one or more genes controlling its production in P. tritici-repentis are unknown. Here, we report the genetic mapping, molecular cloning, and functional analysis of a fungal gene that is required for Ptr ToxC production. The genetic locus controlling the production of Ptr ToxC, termed ToxC, was mapped to a subtelomeric region using segregating biparental populations, genome sequencing, and association analysis. Additional marker analysis further delimited ToxC to a 173-kb region. The predicted genes in the region were examined for presence/absence polymorphism in different races and isolates leading to the identification of a single candidate gene. Functional validation showed that this gene was required but not sufficient for Ptr ToxC production, thus it is designated as ToxC1. ToxC1 encoded a conserved hypothetical protein likely located on the vacuole membrane. The gene was highly expressed during infection, and only one haplotype was identified among 120 isolates sequenced. Our work suggests that Ptr ToxC is not a protein and is likely produced through a cascade of biosynthetic pathway. The identification of ToxC1 is a major step toward revealing the Ptr ToxC biosynthetic pathway and studying its molecular interactions with host factors.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Comparative Genome Analysis Reveals Cis-Regulatory Elements on Gene-Sized Chromosomes of Ciliated Protists.

Gene-sized chromosomes are a distinct feature of the macronuclear genome in ciliated protists known as spirotrichs. These nanochromosomes are often only several kilobase pairs long and contain a coding region for a single gene. However, the ways in which transcription is regulated on nanochromosomes is still largely unknown. Here, we generated macronuclear genome assemblies for two species of Pseudokeronopsis ciliates to better understand transcription regulation on gene-sized chromosomes. We searched within the short subtelomeric regions for potential cis-regulatory elements and identified distinct AT-rich sequences conserved in both species, at both the 5’ and 3’ end of each gene. We further acquired transcriptomic data for these species, which showed the 5’ cis-regulatory element is associated with active gene expression. Gene family evolution analysis suggests nanochromosomes in spirotrichs may originated approximately 900 million years ago. Together our comparative genomic analyses reveal novel insights into the biological roles of cis-regulatory elements on gene-sized chromosomes.

TOR complex 2 contributes to regulation of gene expression via inhibiting Gcn5 recruitment to subtelomeric and DNA replication stress genes.

The fission yeast TOR complex 2 (TORC2) is required for gene silencing at subtelomeric regions and for the induction of gene transcription in response to DNA replication stress. Thus, TORC2 affects transcription regulation both negatively and positively. Whether these two TORC2-dependent functions share a common molecular mechanism is currently unknown. Here, we show that Gad8 physically interacts with proteins that regulate transcription, including subunits of the Spt-Ada-Gcn5-acetyltransferase (SAGA) complex and the BET bromodomain protein Bdf2. We demonstrate that in the absence of TORC2, Gcn5, the histone acetyltransferase subunit of SAGA, accumulates at subtelomeric genes and at non-induced promoters of DNA replication genes. Remarkably, the loss of Gcn5 in TORC2 mutant cells restores gene silencing as well as transcriptional induction in response to DNA replication stress. Loss of Bdf2 alleviates excess of Gcn5 binding in TORC2 mutant cells and also rescues the aberrant regulation of transcription in these cells. Furthermore, the loss of either SAGA or Bdf2 suppresses the sensitivity of TORC2 mutant cells to a variety of stresses, including DNA replication, DNA damage, temperature and nutrient stresses. We suggest a role of TORC2 in transcriptional regulation that is critical for gene silencing and gene induction in response to stress and involves the binding of Gcn5 to the chromatin.