Repeat Array Research Articles

Structural variation and DNA methylation shape the centromere-proximal meiotic crossover landscape in Arabidopsis.

Centromeres load kinetochore complexes onto chromosomes, which mediate spindle attachment and allow segregation during cell division. Although centromeres perform a conserved cellular function, their underlying DNA sequences are highly divergent within and between species. Despite variability in DNA sequence, centromeres are also universally suppressed for meiotic crossover recombination, across eukaryotes. However, the genetic and epigenetic factors responsible for suppression of centromeric crossovers remain to be completely defined. To explore the centromere-proximal meiotic recombination landscape, we map 14,397 crossovers against fully assembled Arabidopsis thaliana (A. thaliana) genomes. A. thaliana centromeres comprise megabase satellite repeat arrays that load nucleosomes containing the CENH3 histone variant. Each chromosome contains a structurally polymorphic region of ~3-4 megabases, which lack crossovers and include the satellite arrays. This polymorphic region is flanked by ~1-2 megabase low-recombination zones. These recombination-suppressed regions are enriched for Gypsy/Ty3 retrotransposons, and additionally contain expressed genes with high genetic diversity that initiate meiotic recombination, yet do not crossover. We map crossovers at high-resolution in proximity to CEN3, which resolves punctate centromere-proximal hotspots that overlap gene islands embedded in heterochromatin. Centromeres are densely DNA methylated and the recombination landscape is remodelled in DNA methylation mutants. We observe that the centromeric low-recombining zones decrease and increase crossovers in CG (met1) and non-CG (cmt3) mutants, respectively, whereas the core non-recombining zones remain suppressed. Our work relates the genetic and epigenetic organization of A. thaliana centromeres and flanking pericentromeric heterochromatin to the zones of crossover suppression that surround the CENH3-occupied satellite repeat arrays.

Genomic Insights of a Methicillin-Resistant Biofilm-Producing Staphylococcus aureus Strain Isolated From Food Handlers.

Methicillin-resistant Staphylococcus aureus (MRSA) is an important zoonotic pathogen associated with a wide range of infections in humans and animals. Thus, the emergence of MRSA clones poses an important threat to human and animal health. This study is aimed at elucidating the genomics insights of a strong biofilm-producing and multidrug-resistant (MDR) S. aureus MTR_BAU_H1 strain through whole-genome sequencing (WGS). The S. aureus MTR_BAU_H1 strain was isolated from food handlers' hand swabs in Bangladesh and phenotypically assessed for antimicrobial susceptibility and biofilm production assays. The isolate was further undergone to high throughput WGS and analysed using different bioinformatics tools to elucidate the genetic diversity, molecular epidemiology, sequence type (ST), antimicrobial resistance, and virulence gene distribution. Phenotypic analyses revealed that the S. aureus MTR_BAU_H1 strain is a strong biofilm-former and carries both antimicrobial resistance (e.g., methicillin resistance; mecA, beta-lactam resistance; blaZ and tetracycline resistance; tetC) and virulence (e.g., sea, tsst, and PVL) genes. The genome of the S. aureus MTR_BAU_H1 belonged to ST1930 that possessed three plasmid replicons (e.g., rep16, rep7c, and rep19), seven prophages, and two clustered regularly interspaced short palindromic repeat (CRISPR) arrays of varying sizes. Phylogenetic analysis showed a close evolutionary relationship between the MTR_BAU_H1 genome and other MRSA clones of diverse hosts and demographics. The MTR_BAU_H1 genome harbours 42 antimicrobial resistance genes (ARGs), 128 virulence genes, and 273 SEED subsystems coding for the metabolism of amino acids, carbohydrates, proteins, cofactors, vitamins, minerals, and lipids. This is the first-ever WGS-based study of a strong biofilm-producing and MDR S. aureus strain isolated from human hand swabs in Bangladesh that unveils new information on the resistomes (ARGs and correlated mechanisms) and virulence potentials that might be linked to staphylococcal pathogenesis in both humans and animals.

Whole genome sequencing of Lacticaseibacillus casei KACC92338 strain with strong antioxidant activity, reveals genes and gene clusters of probiotic and antimicrobial potential.

Lacticaseibacillus casei KACC92338 was originally isolated from Korean raw milk. The antioxidant activities and protective effect in vitro of this strain were evaluated extensively. The results showed that KACC92338 can tolerate hydrogen peroxide up to 2 mM and cell-free supernatant (CFS) had higher scavenging rates for DPPH, hydroxyl radical, reducing power, and iron chelating activities with 95.61 ± 1.59%, 34.10 ± 1.93%, 2.220 ± 0.82 and 81.06 ± 1.06%, respectively. Meanwhile, the CFS showed a protective effect on yeast cells against 10 mM hydrogen peroxide with a survival rate of 76.05 ± 5.65%. To explore the probiotic potential of KACC92338, whole genome assembly and gene clusters with probiotic properties were further analyzed. The genome size was 3,050,901 bp with a 47.96% GC ratio, and 63 contigs. The genome contains 3,048 genes composed of 2,981 coding sequences and 67 RNAs (including 57 tRNAs +9 rRNAs +1 tmRNA). Average Nucleotide Identity and genome-based taxonomy showed that the KACC92338 genome had close similarity with L. casei strains with 96% ANI. Functional annotation by EggNOG and KEGG revealed the presence of numerous genes putatively involved in carbohydrate- and amino acid-transport and metabolism, genetic information processing, and signaling and cellular processes. Additionally, several genes conferring probiotic characteristics such as tolerance to stress, heat, cold, acid, bile salts, oxidative stress, immunomodulation, and adhesion to intestinal epithelium were identified. Notably absent were acquired antibiotic resistance genes, virulence, and pathogenic factors, that prove KACC92338 is a safe strain. Besides, the defense mechanisms of KACC92338 include six prophage regions and three clustered regularly interspaced short palindromic repeat (CRISPR) arrays as acquired immune systems against mobile elements. Further, the BAGEL4 database determined antimicrobial bacteriocin clusters of class IIb: sakacin-P, Enterolysin_A, sactipeptides, and Enterocin X, which suggests the strain could exhibit a wide range of antimicrobial functions. Together, these findings show that the L. casei KACC92338 strain can be a potential probiotic candidate in producing functional fermented foods-, health care- and skin care products- with antioxidant properties. However, a few more mechanistic studies are necessary on the safety assurance and potential application of the strain as a probiotic agent.

Expression, purification, and biochemical analysis of the RNA-DNA hybrid helicase Sen1/SETX from Chaetomium thermophilum.

Yeast Sen1 and its vertebrate ortholog Senataxin (also known as SETX) are RNA-DNA resolving helicases. Sen1 and SETX are implicated in multiple critical nuclear functions not limited to but including DNA replication and repair, RNA processing, and transcription. These>200kDa helicases have a two-domain architecture with an N-terminal regulatory helical repeat array linked to an SF1b helicase motor core via a variable sized central linker of low complexity sequence. Given the size of these proteins, production of milligram quantities of protein that is suitable for biochemical, biophysical, and protein structural analysis has been challenging. To overcome these limitations, we developed a robust selectable high-yield YFP-fusion protein expression method for Sen1 production in mammalian cells, followed by purification on a high-affinity YFP-binding camelid nanobody support. Herein, we detail methods and protocols for the expression and purification of recombinant Sen1 from the thermophilic fungus Chaetomium thermophilum, and the quantitative characterization of its RNA-DNA duplex resolution activity.

Rex1BD and the 14-3-3 protein control heterochromatin organization at tandem repeats by linking RNAi and HDAC.

Tandem DNA repeats are often organized into heterochromatin that is crucial for genome organization and stability. Recent studies revealed that individual repeats within tandem DNA repeats can behave very differently. How DNA repeats are assembled into distinct heterochromatin structures remains poorly understood. Here, we developed a genome-wide genetic screen using a reporter gene at different units in a repeat array. This screen led to identification of a conserved protein Rex1BD required for heterochromatin silencing. Our structural analysis revealed that Rex1BD forms a four-helix bundle structure with a distinct charged electrostatic surface. Mechanistically, Rex1BD facilitates the recruitment of Clr6 histone deacetylase (HDAC) by interacting with histones. Interestingly, Rex1BD also interacts with the 14-3-3 protein Rad25, which is responsible for recruiting the RITS (RNA-induced transcriptional silencing) complex to DNA repeats. Our results suggest that coordinated action of Rex1BD and Rad25 mediates formation of distinct heterochromatin structure at DNA repeats via linking RNAi and HDAC pathways.

Array CGH + SNP in ALL Enhances Detection of Complex Cytogenetic Profiles and Suggests That PRSS1 Loss May Play a Role in Leukemogenesis

Introduction Cytogenetic alterations are essential prognostic biomarkers in patients (pts) with acute lymphoblastic leukemia/lymphoma (ALL), but are limited by the propensity of leukemic blasts to undergo apoptosis in culture, poor-quality chromosomes, and overgrowth by normal cells. Array CGH + SNP has 10,000-fold higher resolution than conventional cytogenetics and can identify submicroscopic gains, losses and copy number neutral loss of heterozygosity (cnLOH) irrespective of mitotic rate. Few studies outside of clinical trials have investigated microarray findings in adults with ALL or correlated microarray findings with clinical outcomes. Methods We retrospectively reviewed pts with ALL and mixed phenotype acute leukemia (MPAL) aged ≥18 years at diagnosis at our institution between 2017 and 2022. Array analysis was performed on diagnostic bone marrow or peripheral blood samples using high resolution CGH + SNP Genetisure cancer array platform (Agilent). Statistical analysis was done using R software (R core team, 2022). Results Baseline characteristics of included pts (n = 60) are shown in Table 1. CGH + SNP analysis identified 93.3% (56/60) of pts with genomic changes with a mean of 7 (range 0 - 28) aberrations per pt. Among pts with a normal karyotype, 81.8% (9/11) had new genomic abnormalities revealed by CGH + SNP, but significantly fewer aberrations were detected per pt compared to those with abnormal karyotypes (mean 3.09 v. 7.67, p = .008). Pts with Ph-like B-ALL/MPAL had significantly more aberrations on array than pts without Ph-like phenotype (mean 12.0 v. 5.85, p.= .002). Among the 78.3% of pts without complex cytogenetics on chromosome analysis, 44.7% (21/47) were found to have a complex genomic profile (≥5 aberrations) with the addition of CGH + SNP. However, there was no difference in CR rates or OS in pts with <5 or ≥5 microarray aberrations (p = .515 and p = .970, respectively). Quantity of copy number alterations (CNAs) as detected on array was not predictive of achieving end of induction CR/CRi/CRh (mean 6.08 in pts who achieved CR/CRi/CRh v. 7.02 in pts with primary refractory disease, p = .310). Hispanic pts had the same number of CNAs as non-Hispanic pts (mean 6.69 v. 6.94, p = .435). On linear regression analysis, there was no correlation between number of CNAs and age (r 2 = -.01, F = 0.1941, p = .661) or BMI (r 2 = 5.21 x10 -5, F = 0.0029, p = .957). The most frequent focal aberrations were detected in CDKN2A/B (n = 19), IKZF1 (n = 17), PRSS1 (n = 11), PAX5 (n = 8), BTG1 (n = 8) and RB1 (n = 8) (Figure 1). RB1 loss was enriched in pts with T-ALL (45.5% v. 6%, p = .003),and occurred at significantly higher rates than previously reported in T-ALL (Schraders, 2009). Loss of PRSS1 was detected in 18.3% (11/60) of the cohort. Homozygous PRSS1 deletions were identified in 3/11 (27.3%) pts. We identified the minimal overlapping deleted region on 7q34 in all 11 pts, which included the PRSS1 gene within the 124 kb region. Five of 6 pts who had repeat arrays performed at the end of induction had resolution of this aberration. The one pt who retained PRSS1 loss after induction had primary refractory disease, but after further treatment and attainment of CR1, this pt cleared the PRSS1 loss. Asparaginase associated pancreatitis, which has been linked to PRSS1 variations, occurred in only 1 pt with PRSS1 loss. Pts harboring PRSS1 loss were significantly less likely to achieve MRD negativity after induction (50% vs. 82%, P = .041) and had numerically more relapses (37.5% vs. 19%, p = .300), but no OS difference was observed (p = .290). Conclusions In our institutional cohort, the addition of CGH + SNP array reclassified almost half (44.7%) of pts without complex cytogenetics to harbor complex genomic changes (≥5 aberrations). Pts with Ph-like ALL had significantly more aberrations on array than those with Ph-positive or -negative B-ALL. Focal PRSS1 loss was associated with reduced rates of MRD clearance at the end of induction. This aberration resolved upon disease remission, raising the possibility that both homozygous and/or heterozygous deletions of 124 kb from the 7q34 chromosomal region of the PRSS1 gene may potentially have a role in ALL pathogenesis. Further studies with larger cohorts utilizing CGH + SNP analyses at additional time points during treatment are needed to explore the potential integration of microarray data for refined ALL risk stratification as well as to elucidate the functional role of PRSS1 in leukemogenesis.

Mitogenome sequences of domestic cats demonstrate lineage expansions and dynamic mutation processes in a mitochondrial minisatellite

BackgroundAs a population genetic tool, mitochondrial DNA is commonly divided into the ~ 1-kb control region (CR), in which single nucleotide variant (SNV) diversity is relatively high, and the coding region, in which selective constraint is greater and diversity lower, but which provides an informative phylogeny. In some species, the CR contains variable tandemly repeated sequences that are understudied due to heteroplasmy. Domestic cats (Felis catus) have a recent origin and therefore traditional CR-based analysis of populations yields only a small number of haplotypes.ResultsTo increase resolution we used Nanopore sequencing to analyse 119 cat mitogenomes via a long-amplicon approach. This greatly improves discrimination (from 15 to 87 distinct haplotypes in our dataset) and defines a phylogeny showing similar starlike topologies within all major clades (haplogroups), likely reflecting post-domestication expansion. We sequenced RS2, a CR tandem array of 80-bp repeat units, placing RS2 array structures within the phylogeny and increasing overall haplotype diversity. Repeat number varies between 3 and 12 (median: 4) with over 30 different repeat unit types differing largely by SNVs. Five SNVs show evidence of independent recurrence within the phylogeny, and seven are involved in at least 11 instances of rapid spread along repeat arrays within haplogroups.ConclusionsIn defining mitogenome variation our study provides key information for the forensic genetic analysis of cat hair evidence, and for the first time a phylogenetically informed picture of tandem repeat variation that reveals remarkably dynamic mutation processes at work in the mitochondrion.

Flavones provide resistance to DUX4-induced toxicity via an mTor-independent mechanism

Facioscapulohumeral muscular dystrophy (FSHD) is among the most common of the muscular dystrophies, affecting nearly 1 in 8000 individuals, and is a cause of profound disability. Genetically, FSHD is linked to the contraction and/or epigenetic de-repression of the D4Z4 repeat array on chromosome 4, thereby allowing expression of the DUX4 gene in skeletal muscle. If the DUX4 transcript incorporates a stabilizing polyadenylation site the myotoxic DUX4 protein will be synthesized, resulting in muscle wasting. The mechanism of toxicity remains unclear, as many DUX4-induced cytopathologies have been described, however cell death does primarily occur through caspase 3/7-dependent apoptosis. To date, most FSHD therapeutic development has focused on molecular methods targeting DUX4 expression or the DUX4 transcript, while therapies targeting processes downstream of DUX4 activity have received less attention. Several studies have demonstrated that inhibition of multiple signal transduction pathways can ameliorate DUX4-induced toxicity, and thus compounds targeting these pathways have the potential to be developed into FSHD therapeutics. To this end, we have screened a group of small molecules curated based on their reported activity in relevant pathways and/or structural relationships with known toxicity-modulating molecules. We have identified a panel of five compounds that function downstream of DUX4 activity to inhibit DUX4-induced toxicity. Unexpectedly, this effect was mediated through an mTor-independent mechanism that preserved expression of ULK1 and correlated with an increase in a marker of active cellular autophagy. This identifies these flavones as compounds of interest for therapeutic development, and potentially identifies the autophagy pathway as a target for therapeutics.

Spt5 C-terminal repeat domain phosphorylation and length negatively regulate heterochromatin through distinct mechanisms.

Heterochromatin is a condensed chromatin structure that represses transcription of repetitive DNA elements and developmental genes, and is required for genome stability. Paradoxically, transcription of heterochromatic sequences is required for establishment of heterochromatin in diverse eukaryotic species. As such, components of the transcriptional machinery can play important roles in establishing heterochromatin. How these factors coordinate with heterochromatin proteins at nascent heterochromatic transcripts remains poorly understood. In the model eukaryote Schizosaccharomyces pombe (S. pombe), heterochromatin nucleation can be coupled to processing of nascent transcripts by the RNA interference (RNAi) pathway, or to other post-transcriptional mechanisms that are RNAi-independent. Here we show that the RNA polymerase II processivity factor Spt5 negatively regulates heterochromatin in S. pombe through its C-terminal domain (CTD). The Spt5 CTD is analogous to the CTD of the RNA polymerase II large subunit, and is comprised of multiple repeats of an amino acid motif that is phosphorylated by Cdk9. We provide evidence that genetic ablation of Spt5 CTD phosphorylation results in aberrant RNAi-dependent nucleation of heterochromatin at an ectopic location, as well as inappropriate spread of heterochromatin proximal to centromeres. In contrast, truncation of Spt5 CTD repeat number enhanced RNAi-independent heterochromatin formation and bypassed the requirement for RNAi. We relate these phenotypes to the known Spt5 CTD-binding factor Prf1/Rtf1. This separation of function argues that Spt5 CTD phosphorylation and CTD length restrict heterochromatin through unique mechanisms. More broadly, our findings argue that length and phosphorylation of the Spt5 CTD repeat array have distinct regulatory effects on transcription.

Optical Genome Mapping for the Molecular Diagnosis of Facioscapulohumeral Muscular Dystrophy: Advancement and Challenges.

Facioscapulohumeral muscular dystrophy (FSHD) is the second most common muscular dystrophy in adults, and it is associated with local D4Z4 chromatin relaxation, mostly via the contraction of the D4Z4 macrosatellite repeat array on chromosome 4q35. In this study, we aimed to investigate the use of Optical Genome Mapping (OGM) as a diagnostic tool for testing FSHD cases from the UK and India and to compare OGM performance with that of traditional techniques such as linear gel (LGE) and Pulsed-field gel electrophoresis (PFGE) Southern blotting (SB). A total of 6 confirmed and 19 suspected FSHD samples were processed with LGE and PFGE, respectively. The same samples were run using a Saphyr Genome-Imaging Instrument (1-color), and the data were analysed using custom EnFocus FSHD analysis. OGM was able to confirm the diagnosis of FSHD1 in all FSHD1 cases positive for SB (n = 17), and D4Z4 sizing highly correlated with PFGE-SB (p < 0.001). OGM correctly identified cases with mosaicism for the repeat array contraction (n = 2) and with a duplication of the D4Z4 repeat array. OGM is a promising new technology able to unravel structural variants in the genome and seems to be a valid tool for diagnosing FSHD1.

Characterization of integrated Marek's disease virus genomes supports a model of integration by homology-directed recombination and telomere-loop-driven excision.

Marek's disease virus (MDV) is a ubiquitous chicken pathogen that inflicts a large economic burden on the poultry industry, despite worldwide vaccination programs. MDV is only partially controlled by available vaccines, and the virus retains the ability to replicate and spread between vaccinated birds. Following an initial infection, MDV enters a latent state and integrates into host telomeres and this may be a prerequisite for malignant transformation, which is usually fatal. To understand the mechanism that underlies the dynamic relationship between integrated-latent and reactivated MDV, we have characterized integrated MDV (iMDV) genomes and their associated telomeres. This revealed a single orientation among iMDV genomes and the loss of some terminal sequences that is consistent with integration by homology-directed recombination and excision via a telomere-loop-mediated process.

Expanding the repertoire of human tandem repeat RNA-binding proteins.

Protein regions consisting of arrays of tandem repeats are known to bind other molecular partners, including nucleic acid molecules. Although the interactions between repeat proteins and DNA are already widely explored, studies characterising tandem repeat RNA-binding proteins are lacking. We performed a large-scale analysis of human proteins devoted to expanding the knowledge about tandem repeat proteins experimentally reported as RNA-binding molecules. This work is timely because of the release of a full set of accurate structural models for the human proteome amenable to repeat detection using structural methods. The main goal of our analysis was to build a comprehensive set of human RNA-binding proteins that contain repeats at the sequence or structure level. Our results showed that the combination of sequence and structural methods finds significantly more tandem repeat proteins than either method alone. We identified 219 tandem repeat proteins that bind RNA molecules and characterised the overlap between repeat regions and RNA-binding regions as a first step towards assessing their functional relationship. We observed differences in the characteristics of repeat regions predicted by sequence-based or structure-based methods in terms of their sequence composition, their functions and their protein domains.

Autosomal dominant in cis D4Z4 repeat array duplication alleles in facioscapulohumeral dystrophy.

Facioscapulohumeral dystrophy (FSHD) has a unique genetic aetiology resulting in partial chromatin relaxation of the D4Z4 macrosatellite repeat array on 4qter. This D4Z4 chromatin relaxation facilitates inappropriate expression of the transcription factor DUX4 in skeletal muscle. DUX4 is encoded by a retrogene that is embedded within the distal region of the D4Z4 repeat array. In the European population, the D4Z4 repeat array is usually organized in a single array that ranges between 8 and 100 units. D4Z4 chromatin relaxation and DUX4 derepression in FSHD is most often caused by repeat array contraction to 1-10 units (FSHD1) or by a digenic mechanism requiring pathogenic variants in a D4Z4 chromatin repressor like SMCHD1, combined with a repeat array between 8 and 20 units (FSHD2). With a prevalence of 1.5% in the European population, in cis duplications of the D4Z4 repeat array, where two adjacent D4Z4 arrays are interrupted by a spacer sequence, are relatively common but their relationship to FSHD is not well understood. In cis duplication alleles were shown to be pathogenic in FSHD2 patients; however, there is inconsistent evidence for the necessity of an SMCHD1 mutation for disease development. To explore the pathogenic nature of these alleles we compared in cis duplication alleles in FSHD patients with or without pathogenic SMCHD1 variant. For both groups we showed duplication-allele-specific DUX4 expression. We studied these alleles in detail using pulsed-field gel electrophoresis-based Southern blotting and molecular combing, emphasizing the challenges in the characterization of these rearrangements. Nanopore sequencing was instrumental to study the composition and methylation of the duplicated D4Z4 repeat arrays and to identify the breakpoints and the spacer sequence between the arrays. By comparing the composition of the D4Z4 repeat array of in cis duplication alleles in both groups, we found that specific combinations of proximal and distal repeat array sizes determine their pathogenicity. Supported by our algorithm to predict pathogenicity, diagnostic laboratories should now be furnished to accurately interpret these in cis D4Z4 repeat array duplications, alleles that can easily be missed in routine settings.

Tandem NBPF 3mer HORs (Olduvai triplets) in Neanderthal and two novel HOR tandem arrays in human chromosome 1 T2T-CHM13 assembly

It is known that the ~ 1.6 kb Neuroblastoma BreakPoint Family (NBPF) repeats are human specific and contributing to cognitive capabilities, with increasing frequency in higher order repeat 3mer HORs (Olduvai triplets). From chimpanzee to modern human there is a discontinuous jump from 0 to ~ 50 tandemly organized 3mer HORs. Here we investigate the structure of NBPF 3mer HORs in the Neanderthal genome assembly of Pääbo et al., comparing it to the results obtained for human hg38.p14 chromosome 1. Our findings reveal corresponding NBPF 3mer HOR arrays in Neanderthals with slightly different monomer structures and numbers of HOR copies compared to humans. Additionally, we compute the NBPF 3mer HOR pattern for the complete telomere-to-telomere human genome assembly (T2T-CHM13) by Miga et al., identifying two novel tandem arrays of NBPF 3mer HOR repeats with 5 and 9 NBPF 3mer HOR copies. We hypothesize that these arrays correspond to novel NBPF genes (here referred to as NBPFA1 and NBPFA2). Further improving the quality of the Neanderthal genome using T2T-CHM13 as a reference would be of great interest in determining the presence of such distant novel NBPF genes in the Neanderthal genome and enhancing our understanding of human evolution.

Measuring open chromatin and DNA methylation in repeat arrays.

Measuring open chromatin and DNA methylation in repeat arrays.

Assembly of 43 human Y chromosomes reveals extensive complexity and variation.

The prevalence of highly repetitive sequences within the human Y chromosome has prevented its complete assembly to date1 and led to its systematic omission from genomic analyses. Here we present de novo assemblies of 43 Y chromosomes spanning 182,900 years of human evolution and report considerable diversity in size and structure. Half of the male-specific euchromatic region is subject to large inversions with a greater than twofold higher recurrence rate compared with all other chromosomes2. Ampliconic sequences associated with these inversions show differing mutation rates that are sequence context dependent, and some ampliconic genes exhibit evidence for concerted evolution with the acquisition and purging of lineage-specific pseudogenes. The largest heterochromatic region in the human genome, Yq12, is composed of alternating repeat arrays that show extensive variation in the number, size and distribution, but retain a 1:1 copy-number ratio. Finally, our data suggest that the boundary between the recombining pseudoautosomal region 1 and the non-recombining portions of the X and Y chromosomes lies 500 kb away from the currently established1 boundary. The availability of fully sequence-resolved Y chromosomes from multiple individuals provides a unique opportunity for identifying new associations of traits with specific Y-chromosomal variants and garnering insights into the evolution and function of complex regions of the human genome.

DYW cytidine deaminase domains have a long-range impact on RNA recognition by the PPR array of chimeric plant C-to-U RNA editing factors and strongly affect target selection.

The protein factors for the specific C-to-U RNA editing events in plant mitochondria and chloroplasts possess unique arrays of RNA-binding pentatricopeptide repeats (PPRs) linked to carboxy-terminal cytidine deaminase DYW domains via the extension motifs E1 and E2. The E1 and E2 motifs have distant similarities to tetratricopeptide repeats known to mediate protein-protein interactions but their precise function is unclear. Here, we investigate the tolerance of PPR56 and PPR65, two functionally characterized RNA editing factors of the moss Physcomitrium patens, for the creation of chimeras by variably replacing their C-terminal protein regions. Making use of a heterologous RNA editing assay system in Escherichia coli we find that heterologous DYW domains can strongly restrict or widen the spectrum of off-targets in the bacterial transcriptome for PPR56. Surprisingly, our data suggest that these changes are not only caused by the preference of a given heterologous DYW domain for the immediate sequence environment of the cytidine to be edited but also by a long-range impact on the nucleotide selectivity of the upstream PPRs.

Emergence of high-level tigecycline resistance due to the amplification of a tet(A) gene variant in clinical carbapenem-resistant Klebsiella pneumoniae

ObjectiveTo investigate the prevalence of a tet(A) gene variant and its role in developing high-level tigecycline resistance among carbapenem-resistant Klebsiella pneumoniae (CRKP) clinical isolates. MethodsThe mechanism of high-level tigecycline resistance in CRKP mediated by a tet(A) variant was explored by induction experiments, antimicrobial susceptibility testing, whole-genome sequencing and bioinformatics analysis. The amplification and overexpression of the tet(A) variant were measured by the determination of sequencing depth, gene copy numbers, and qRT-PCR. ResultsA high rate (62.1%, 998/1607) of tet(A) variant carriage was observed among 1607 CRKP clinical isolates from Henan Province, China. High-level tigecycline resistance could rapidly develop by the amplification of the tet(A) variant in these isolates. The analysis of the raw sequencing data and the plasmid mapping depth revealed that the ΔtnpA homologous sequence of Tn1721 supports the amplification of the region that harbours the tet(A) variant by forming a large number of repeat arrays through translocatable units (TUs). Moreover, the epidemiological analysis of tet(A) variant-carrying structures among 1607 clinical CRKPs showed that the TU structure is widely present. ConclusionThe presence of a tigecycline resistance-mediating tet(A) variant in CRKP clinical isolates represents a greater health concern than initially thought and should be monitored consistently.

Cycles of satellite and transposon evolution in Arabidopsis centromeres.

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

TRASH: Tandem Repeat Annotation and Structural Hierarchy.

The advent of long-read DNA sequencing is allowing complete assembly of highly repetitive genomic regions for the first time, including the megabase-scale satellite repeat arrays found in many eukaryotic centromeres. The assembly of such repetitive regions creates a need for their de novo annotation, including patterns of higher order repetition. To annotate tandem repeats, methods are required that can be widely applied to diverse genome sequences, without prior knowledge of monomer sequences. Tandem Repeat Annotation and Structural Hierarchy (TRASH) is a tool that identifies and maps tandem repeats in nucleotide sequence, without prior knowledge of repeat composition. TRASH analyses a fasta assembly file, identifies regions occupied by repeats and then precisely maps them and their higher order structures. To demonstrate the applicability and scalability of TRASH for centromere research, we apply our method to the recently published Col-CEN genome of Arabidopsis thaliana and the complete human CHM13 genome. TRASH is freely available at:https://github.com/vlothec/TRASH and supported on Linux.