TRDB--the Tandem Repeats Database

The role of centromeric repeats and transcripts in kinetochore assembly and function

Centromeres are the chromosomal domains, where the kinetochore protein complex is formed, mediating proper segregation of chromosomes during cell division. Although the function of centromeres has remained conserved during evolution, centromeric DNA is highly variable, even in closely related species. In addition, the composition of the kinetochore complexes varies among organisms. Therefore, it is assumed that the centromeric position is determined epigenetically, and the centromeric histone H3 (CENH3) serves as an epigenetic marker. The loading of CENH3 onto centromeres depends on centromere-licensing factors, chaperones, and transcription of centromeric repeats. Several proteins that regulate CENH3 loading and kinetochore assembly interact with the centromeric transcripts and DNA in a sequence-independent manner. However, the functional aspects of these interactions are not fully understood. This review discusses the variability of centromeric sequences in different organisms and the regulation of their transcription through the RNA Pol II and RNAi machinery. The data suggest that the interaction of proteins involved in CENH3 loading and kinetochore assembly with centromeric DNA and transcripts plays a role in centromere, and possibly neocentromere, formation in a sequence-independent manner.

Analysis and benchmarking of small and large genomic variants across tandem repeats

Tandem repeats (TRs) are highly polymorphic in the human genome, have thousands of associated molecular traits and are linked to over 60 disease phenotypes. However, they are often excluded from at-scale studies because of challenges with variant calling and representation, as well as a lack of a genome-wide standard. Here, to promote the development of TR methods, we created a catalog of TR regions and explored TR properties across 86 haplotype-resolved long-read human assemblies. We curated variants from the Genome in a Bottle (GIAB) HG002 individual to create a TR dataset to benchmark existing and future TR analysis methods. We also present an improved variant comparison method that handles variants greater than 4 bp in length and varying allelic representation. The 8.1% of the genome covered by the TR catalog holds ~24.9% of variants per individual, including 124,728 small and 17,988 large variants for the GIAB HG002 'truth-set' TR benchmark. We demonstrate the utility of this pipeline across short-read and long-read technologies.

The analysis of complete genome sequence and comparative genomics of Vibrio parahaemolyticus LF1113 in Hainan

Vibrio parahaemolyticus is a Gram-negative, halophilic and polymorphic coccobacillus. It is world-widely distributed and has resulted in great economic losses since its first appearance. In this study, a pathogenic strain was isolated from diseased pearl gentian grouper and identified as V. parahaemolyticus based on the sequencing results of 16S rDNA gene. In order to gain a comprehensive understanding of this isolation, the whole genome sequencing was conducted. Phylogenetic analysis of the complete genomes of 16 Vibrio species showed that LF1113, ATCC17802, ATCC33787, 2210633, FORC 004, and 160807 were the most closely related. Animal experiments demonstrated that the isolated LF1113 strain was pathogenic in a fish model. This study is the first study to describe the complete genome sequence of a V. parahaemolyticus isolate, which infected pearl gentian grouper from an outbreak in a fish factory farm in Hainan. The results will expand our understanding of genetic characteristics, pathogenesis, diagnostics and disease prevention of V. parahaemolyticus, and lay the foundation for further study.

Bacterial host adaptation through sequence and structural variations of a single type III effector gene

Molecular mechanisms underlying quantitative variations of pathogenicity remain elusive. Here, we identified the Xanthomonas campestris XopJ6 effector that triggers disease resistance in cauliflower and Arabidopsis thaliana. XopJ6 is a close homolog of the Ralstoniapseudosolanacearum PopP2 YopJ family acetyltransferase. XopJ6 is recognized by the RRS1-R/RPS4 NLR pair that integrates a WRKY decoy domain mimicking effector targets. We identified a XopJ6 natural variant carrying a single residue substitution in XopJ6 WRKY-binding site that disrupts interaction with WRKY proteins. This mutation allows XopJ6 to evade immune perception while retaining some XopJ6 virulence functions. Interestingly, xopJ6 resides in a Tn3-family transposon likely contributing to xopJ6 copy number variation (CNV). Using synthetic biology, we demonstrate that xopJ6 CNV tunes pathogen virulence on Arabidopsis through gene dosage-mediated modulation of xopJ6 expression. Together, our findings highlight how sequence and structural genetic variations restricted at a particular effector gene contribute to bacterial host adaptation.

Benchmarking of small and large variants across tandem repeats

Tandem repeats (TRs) are highly polymorphic in the human genome, have thousands of associated molecular traits, and are linked to over 60 disease phenotypes. However, their complexity often excludes them from at-scale studies due to challenges with variant calling, representation, and lack of a genome-wide standard. To promote TR methods development, we create a comprehensive catalog of TR regions and explore its properties across 86 samples. We then curate variants from the GIAB HG002 individual to create a tandem repeat benchmark. We also present a variant comparison method that handles small and large alleles and varying allelic representation. The 8.1% of the genome covered by the TR catalog holds ∼24.9% of variants per individual, including 124,728 small and 17,988 large variants for the GIAB HG002 TR benchmark. We work with the GIAB community to demonstrate the utility of this benchmark across short and long read technologies.

Polymeric structure of the Cannabis sativa L. mitochondrial genome identified with an assembly graph model

Cannabis sativa L. belongs to the family Cannabaceae in Rosales. It has been widely used as medicines, building materials, and textiles. Elucidating its genome is critical for molecular breeding and synthetic biology study. Many studies have shown that the mitochondrial genomes (mitogenomes) and even chloroplast genomes (plastomes) had complex polymeric structures. Using the Nanopore sequencing platform, we sequenced, assembled, and analyzed its mitogenome and plastome. The resulting unitig graph suggested that the mitogenome had a complex polymeric structure. However, a gap-free, circular sequence was further assembled from the unitig graph. In contrast, a circular sequence representing the plastome was obtained. The mitogenome major conformation was 415,837 bp long, and the plastome was 153,927 bp long. To test if the repeat sequences promote recombination, which corresponds to the branch points in the structure, we tested the sequences around repeats by long-read mapping. Among 208 pairs of predicted repeats, the mapping results supported the presence of cross-over around 25 pairs of repeats. Subsequent PCR amplification confirmed the presence of cross-over around 15 of the 25 repeats. By comparing the mitogenome and plastome sequences, we identified 19 mitochondria plastid DNAs, including seven complete genes (trnW-CCA, trnP-UGG, psbJ, trnN-GUU, trnD-GUC, trnH-GUG, trnM-CAU) and nine gene fragments. Furthermore, the selective pressure analysis results showed that five genes (atp1, ccmB, ccmC, cox1, nad7) had 19 positively selected sites. Lastly, we predicted 28 RNA editing sites. A total of 8 RNA editing sites located in the coding regions were successfully validated by PCR amplification and Sanger sequencing, of which four were synonymous, and four were nonsynonymous. In particular, the RNA editing events appeared to be tissue-specific in C. sativa mitogenome. In summary, we have confirmed the major confirmation of C. sativa mitogenome and characterized its structural features in detail. These results provide critical information for future variety breeding and resource development for C. sativa.

ScrepYard: An online resource for disulfide-stabilized tandem repeat peptides

Receptor avidity through multivalency is a highly sought-after property of ligands. While readily available in nature in the form of bivalent antibodies, this property remains challenging to engineer in synthetic molecules. The discovery of several bivalent venom peptides containing two homologous and independently folded domains (in a tandem repeat arrangement) has provided a unique opportunity to better understand the underpinning design of multivalency in multimeric biomolecules, as well as how naturally occurring multivalent ligands can be identified. In previous work, we classified these molecules as a larger class termed secreted cysteine-rich repeat-proteins (SCREPs). Here, we present an online resource; ScrepYard, designed to assist researchers in identification of SCREP sequences of interest and to aid in characterizing this emerging class of biomolecules. Analysis of sequences within the ScrepYard reveals that two-domain tandem repeats constitute the most abundant SCREP domain architecture, while the interdomain "linker" regions connecting the functional domains are found to be abundant in amino acids with short or polar sidechains and contain an unusually high abundance of proline residues. Finally, we demonstrate the utility of ScrepYard as a virtual screening tool for discovery of putatively multivalent peptides, by using it as a resource to identify a previously uncharacterized serine protease inhibitor and confirm its predicted activity using an enzyme assay.

Mapping separase-mediated cleavage in situ

Separase is a protease that performs critical functions in the maintenance of genetic homeostasis. Among them, the cleavage of the meiotic cohesin during meiosis is a key step in producing gametes in eukaryotes. However, the exact chromosomal localization of this proteolytic cleavage was not addressed due to the lack of experimental tools. To this end, we developed a method based on monoclonal antibodies capable of recognizing the predicted neo-epitopes produced by separase-mediated proteolysis in the RAD21 and REC8 cohesin subunits. To validate the epigenomic strategy of mapping cohesin proteolysis, anti-RAD21 neo-epitopes antibodies were used in ChIP-On-ChEPseq analysis of human cells undergoing mitotic anaphase. Second, a similar analysis applied for mapping of REC8 cleavage in germline cells in Macaque showed a correlation with a subset of alpha-satellites and other repeats, directly demonstrating that the site-specific mei-cohesin proteolysis hotspots are coincident but not identical with centromeres. The sequences for the corresponding immunoglobulin genes show a convergence of antibodies with close specificity. This approach could be potentially used to investigate cohesin ring opening events in other chromosomal locations, if applied to single cells.

Characterization of the mitochondrial genomes of Bradysia hygida, Phytosciara flavipes and Trichosia splendens (Diptera: Sciaridae) and novel insights on the control region of sciarid mitogenomes

Sciarids, also called "fungus gnats" are small, almost entirely dark-coloured insects. Sciarid larvae feed on different substrates and can infest agricultural crops and mushroom nurseries, causing economic losses. Of the 2174 Diptera mitogenome sequences currently available in GenBank, only eight are from the Sciaridae family, none of which are complete circular molecules. Here we describe the mitogenome sequences of three sciarid species: Phytosciara flavipes, Trichosia splendens and Bradysia hygida and provide novel insights on the control region of sciarid mitogenomes. The assembled mitogenomes range from 16,062 bp in P. flavipes to 17,095 bp in B. hygida. All 13 protein coding genes, 22 tRNAs and 2 rRNAs characteristic of insect mitogenomes were identified, but the sequence of the control region could not be determined. Experimental results suggest that the B. hygida control region is about 21 kb long resulting in a 37 kb long mitogenome which constitutes the largest insect mitochondrial genome described so far. Phylogenetic analysis using all Bibionomorpha mitogenome sequences available in GenBank strongly supports the Sciaridae monophyly and led to the identification of species and subfamily specific gene rearrangements. Our study extends the knowledge of this large and diverse insect family that includes agricultural pest species.

Genome Analyses of Two Blueberry Pathogens: Diaporthe amygdali CAA958 and Diaporthe eres CBS 160.32

The genus Diaporthe includes pathogenic species distributed worldwide and affecting a wide variety of hosts. Diaporthe amygdali and Diaporthe eres have been found to cause cankers, dieback, or twig blights on economically important crops such as soybean, almond, grapevine, and blueberry. Despite their importance as plant pathogens, the strategies of species of Diaporthe to infect host plants are poorly explored. To provide a genomic basis of pathogenicity, the genomes of D. amygdali CAA958 and D. eres CBS 160.32 were sequenced and analyzed. Cellular transporters involved in the transport of toxins, ions, sugars, effectors, and genes implicated in pathogenicity were detected in both genomes. Hydrolases and oxidoreductases were the most prevalent carbohydrate-active enzymes (CAZymes). However, analyses of the secreted proteins revealed that the secretome of D. eres CBS 160.32 is represented by 5.4% of CAZymes, whereas the secreted CAZymes repertoire of D. amygdali CAA958 represents 29.1% of all secretomes. Biosynthetic gene clusters (BGCs) encoding compounds related to phytotoxins and mycotoxins were detected in D. eres and D. amygdali genomes. The core gene clusters of the phytotoxin Fusicoccin A in D. amygdali are reported here through a genome-scale assembly. Comparative analyses of the genomes from 11 Diaporthe species revealed an average of 874 CAZymes, 101 secondary metabolite BGCs, 1640 secreted proteins per species, and genome sizes ranging from 51.5 to 63.6 Mbp. This study offers insights into the overall features and characteristics of Diaporthe genomes. Our findings enrich the knowledge about D. eres and D. amygdali, which will facilitate further research into the pathogenicity mechanisms of these species.

Ribodysgenesis: sudden genome instability in the yeast Saccharomyces cerevisiae arising from RNase H2 cleavage at genomic-embedded ribonucleotides

Ribonucleotides can be incorporated into DNA during replication by the replicative DNA polymerases. These aberrant DNA subunits are efficiently recognized and removed by Ribonucleotide Excision Repair, which is initiated by the heterotrimeric enzyme RNase H2. While RNase H2 is essential in higher eukaryotes, the yeast Saccharomyces cerevisiae can survive without RNase H2 enzyme, although the genome undergoes mutation, recombination and other genome instability events at an increased rate. Although RNase H2 can be considered as a protector of the genome from the deleterious events that can ensue from recognition and removal of embedded ribonucleotides, under conditions of high ribonucleotide incorporation and retention in the genome in a RNase H2-negative strain, sudden introduction of active RNase H2 causes massive DNA breaks and genome instability in a condition which we term 'ribodysgenesis'. The DNA breaks and genome instability arise solely from RNase H2 cleavage directed to the ribonucleotide-containing genome. Survivors of ribodysgenesis have massive loss of heterozygosity events stemming from recombinogenic lesions on the ribonucleotide-containing DNA, with increases of over 1000X from wild-type. DNA breaks are produced over one to two divisions and subsequently cells adapt to RNase H2 and ribonucleotides in the genome and grow with normal levels of genome instability.

Insights into the inhibition of type I-F CRISPR-Cas system by a multifunctional anti-CRISPR protein AcrIF24

CRISPR-Cas systems are prokaryotic adaptive immune systems and phages use anti-CRISPR proteins (Acrs) to counteract these systems. Here, we report the structures of AcrIF24 and its complex with the crRNA-guided surveillance (Csy) complex. The HTH motif of AcrIF24 can bind the Acr promoter region and repress its transcription, suggesting its role as an Aca gene in self-regulation. AcrIF24 forms a homodimer and further induces dimerization of the Csy complex. Apart from blocking the hybridization of target DNA to the crRNA, AcrIF24 also induces the binding of non-sequence-specific dsDNA to the Csy complex, similar to AcrIF9, although this binding seems to play a minor role in AcrIF24 inhibitory capacity. Further structural and biochemical studies of the Csy-AcrIF24-dsDNA complexes and of AcrIF24 mutants reveal that the HTH motif of AcrIF24 and the PAM recognition loop of the Csy complex are structural elements essential for this non-specific dsDNA binding. Moreover, AcrIF24 and AcrIF9 display distinct characteristics in inducing non-specific DNA binding. Together, our findings highlight a multifunctional Acr and suggest potential wide distribution of Acr-induced non-specific DNA binding.

Phages use anti-CRISPR proteins (Acrs) to counteract the bacterial CRISPR-Cas systems. Here, the authors characterize AcrIF24, which functions as an Aca (Acr-associated) to repress and regulate its own transcription, dimerizes the Csy complex, blocks the hybridization of target DNA, and tethers non-sequence-specific DNA to the Csy complex.

Bacterial N4-methylcytosine as an epigenetic mark in eukaryotic DNA

DNA modifications are used to regulate gene expression and defend against invading genetic elements. In eukaryotes, modifications predominantly involve C5-methylcytosine (5mC) and occasionally N6-methyladenine (6mA), while bacteria frequently use N4-methylcytosine (4mC) in addition to 5mC and 6mA. Here we report that 4mC can serve as an epigenetic mark in eukaryotes. Bdelloid rotifers, tiny freshwater invertebrates with transposon-poor genomes rich in foreign genes, lack canonical eukaryotic C5-methyltransferases for 5mC addition, but encode an amino-methyltransferase, N4CMT, captured from bacteria >60 Mya. N4CMT deposits 4mC at active transposons and certain tandem repeats, and fusion to a chromodomain shapes its “histone-read-DNA-write” architecture recognizing silent chromatin marks. Furthermore, amplification of SETDB1 H3K9me3 histone methyltransferases yields variants preferentially binding 4mC-DNA, suggesting “DNA-read-histone-write” partnership to maintain chromatin-based silencing. Our results show how non-native DNA methyl groups can reshape epigenetic systems to silence transposons and demonstrate the potential of horizontal gene transfer to drive regulatory innovation in eukaryotes.

Eukaryotic DNA can be methylated as 5-methylcytosine and N6-methyladenine, but whether other forms of DNA methylation occur has been controversial. Here the authors show that a bacterial DNA methyltransferase was acquired >60 Mya in bdelloid rotifers that catalyzes N4-methylcytosine addition and is involved in suppression of transposon proliferation.

The Genome of Rhyzopertha dominica (Fab.) (Coleoptera: Bostrichidae): Adaptation for Success

The lesser grain borer, Rhyzopertha dominica (F.) (Coleoptera: Bostrichidae), is a major global pest of cereal grains. Infestations are difficult to control as larvae feed inside grain kernels, and many populations are resistant to both contact insecticides and fumigants. We sequenced the genome of R. dominica to identify genes responsible for important biological functions and develop more targeted and efficacious management strategies. The genome was assembled from long read sequencing and long-range scaffolding technologies. The genome assembly is 479.1 Mb, close to the predicted genome size of 480.4 Mb by flow cytometry. This assembly is among the most contiguous beetle assemblies published to date, with 139 scaffolds, an N₅₀ of 53.6 Mb, and L₅₀ of 4, indicating chromosome-scale scaffolds. Predicted genes from biologically relevant groups were manually annotated using transcriptome data from adults and different larval tissues to guide annotation. The expansion of carbohydrase and serine peptidase genes suggest that they combine to enable efficient digestion of cereal proteins. A reduction in the copy number of several detoxification gene families relative to other coleopterans may reflect the low selective pressure on these genes in an insect that spends most of its life feeding internally. Chemoreceptor genes contain elevated numbers of pseudogenes for odorant receptors that also may be related to the recent ontogenetic shift of R. dominica to a diet consisting primarily of stored grains. Analysis of repetitive sequences will further define the evolution of bostrichid beetles compared to other species. The data overall contribute significantly to coleopteran genetic research.

Genome-scale patterns in the loss of heterozygosity incidence in Saccharomyces cerevisiae

Former studies have established that loss of heterozygosity can be a key driver of sequence evolution in unicellular eukaryotes and tissues of metazoans. However, little is known about whether the distribution of loss of heterozygosity events is largely random or forms discernible patterns across genomes. To initiate our experiments, we introduced selectable markers to both arms of all chromosomes of the budding yeast. Subsequent extensive assays, repeated over several genetic backgrounds and environments, provided a wealth of information on the genetic and environmental determinants of loss of heterozygosity. Three findings stand out. First, the number of loss of heterozygosity events per unit time was more than 25 times higher for growing than starving cells. Second, loss of heterozygosity was most frequent when regions of homology around a recombination site were identical, about a half-% sequence divergence was sufficient to reduce its incidence. Finally, the density of loss of heterozygosity events was highly dependent on the genome's physical architecture. It was several-fold higher on short chromosomal arms than on long ones. Comparably large differences were seen within a single arm where regions close to a centromere were visibly less affected than regions close, though usually not strictly adjacent, to a telomere. We suggest that the observed uneven distribution of loss of heterozygosity events could have been caused not only by an uneven density of initial DNA damages. Location-depended differences in the mode of DNA repair, or its effect on fitness, were likely to operate as well.

Comparative genomic analysis of Escherichia coli strains obtained from continuous imipenem stress evolution

The carbapenem-resistant Escherichia coli (E. coli) has aroused increasing attention worldwide, especially in terms of imipenem (IMP) resistance. The molecular mechanism of IMP resistance remains unclear. This study aimed to explore the resistance mechanisms of IMP in E. coli. Susceptible Sx181-0-1 strain was induced into resistance strains by adaptive laboratory evolution. The drug resistance spectrum was measured using the disk diffusion and microbroth dilution methods. Whole-genome sequencing and resequencing were used to analyze the non-synonymous single-nucleotide polymorphisms (nsSNPs) between the primary susceptible strain and resistant strains. The expression levels of these genes with nsSNPs were identified by real-time quantitative PCR (RT-qPCR). Resistance phenotype appeared in the induced 15th generation (induction time = 183 h). Sx181-32 and Sx181-256, which had minimum inhibitory concentrations of IMP of 8 and 64 µg mL-1, were isolated during continuous subculture exposed to increasing concentrations of IMP, respectively. Nineteen nsSNPs were observed both in Sx181-32 and Sx181-256, including rpsU, sdaC, zwf, ttuC, araJ, dacC, mrdA, secF, dacD, lpxD, mrcB, ftsI, envZ, and two unknown function genes (orf01892 and orf01933). Among these 15 genes, five genes (dacC, mrdA, lpxD, mrcB, and ftsI) were mainly involved in cell wall synthesis. The mrdA (V338A, L378P, and M574I) and mrcB (P784L, A736V, and T708A) had three amino acid substitutions, respectively. The expression levels of rpsU, ttuC and orf01933 were elevated in both Sx181-32 and Sx181-256 compared to Sx181-0-1. The expression levels of these genes were elevated in Sx181-256, except for araJ. Bacteria developed resistance to antimicrobials by regulating various biological processes, among which the most involved is the cell wall synthesis (dacC, mrdA, lpxD, mrcB, and ftsI). The combination mutations of mrdA, envZ, and ftsI genes may increase the resistance to IMP. Our study could improve the understanding of the molecular mechanism underlying the IMP resistance of E. coli.

A CTCF-Binding Element and Histone Deacetylation Cooperatively Maintain Chromatin Loops, Linking to Long-Range Gene Regulation in Cancer Genomes

Background

Genes spanning long chromosomal domains are coordinately regulated in human genome, which contribute to global gene dysregulation and carcinogenesis in cancer. It has been noticed that epigenetic modification and chromatin architecture may participate in the regulation process. However, the regulation patterns and functional elements of long-range gene regulation are unclear.

Methods

Based on the clinical transcriptome data from different tumor sets, a novel expressional correlation analysis pipeline was performed to classify the co-regulated regions and subsets of intercorrelated regions. The GLAM2 program was used to predict conserved DNA elements that enriched in regions. Two conserved elements were selected to delete in Ishikawa and HeLa cells by CRISPR-Cas9. SAHA treatment and HDAC knockdown were used to change the histone acetylation status. Using qPCR, MTT, and scratch healing assay, we evaluate the effect on gene expression and cancer cell phenotype. By DNA pull-down and ChIP, the element-binding proteins were testified. 3C and 3D-FISH were performed to depict the alteration in chromatin architecture.

Results

In multiple cancer genomes, we classified subsets of coordinately regulated regions (sub-CRRs) that possibly shared the same regulatory mechanisms and exhibited similar expression patterns. A new conserved DNA element (CRE30) was enriched in sub-CRRs and associated with cancer patient survival. CRE30 could restrict gene regulation in sub-CRRs and affect cancer cell phenotypes. DNA pull-down showed that multiple proteins including CTCF were recruited on the CRE30 locus, and ChIP assay confirmed the CTCF-binding signals. Subsequent results uncovered that as an essential element, CRE30 maintained chromatin loops and mediated a compact chromatin architecture. Moreover, we found that blocking global histone deacetylation induced chromatin loop disruption and CTCF dropping in the region containing CRE30, linked to promoted gene regulation. Additionally, similar effects were observed with CRE30 deletion in another locus of chromosome 8.

Conclusions

Our research clarified a new functional element that recruits CTCF and collaborates with histone deacetylation to maintain high-order chromatin organizations, linking to long-range gene regulation in cancer genomes. The findings highlight a close relationship among conserved DNA element, epigenetic modification, and chromatin architecture in long-range gene regulation process.

Genomic and Metabolomic Analyses of the Marine Fungus Emericellopsis cladophorae: Insights into Saltwater Adaptability Mechanisms and Its Biosynthetic Potential

The genus Emericellopsis is found in terrestrial, but mainly in marine, environments with a worldwide distribution. Although Emericellopsis has been recognized as an important source of bioactive compounds, the range of metabolites expressed by the species of this genus, as well as the genes involved in their production are still poorly known. Untargeted metabolomics, using UPLC- QToF–MS/MS, and genome sequencing (Illumina HiSeq) was performed to unlock E. cladophorae MUM 19.33 chemical diversity. The genome of E. cladophorae is 26.9 Mb and encodes 8572 genes. A large set of genes encoding carbohydrate-active enzymes (CAZymes), secreted proteins, transporters, and secondary metabolite biosynthetic gene clusters were identified. Our analysis also revealed genomic signatures that may reflect a certain fungal adaptability to the marine environment, such as genes encoding for (1) the high-osmolarity glycerol pathway; (2) osmolytes’ biosynthetic processes; (3) ion transport systems, and (4) CAZymes classes allowing the utilization of marine polysaccharides. The fungal crude extract library constructed revealed a promising source of antifungal (e.g., 9,12,13-Trihydroxyoctadec-10-enoic acid, hymeglusin), antibacterial (e.g., NovobiocinA), anticancer (e.g., daunomycinone, isoreserpin, flavopiridol), and anti-inflammatory (e.g., 2’-O-Galloylhyperin) metabolites. We also detected unknown compounds with no structural match in the databases used. The metabolites’ profiles of E. cladophorae MUM 19.33 fermentations were salt dependent. The results of this study contribute to unravel aspects of the biology and ecology of this marine fungus. The genome and metabolome data are relevant for future biotechnological exploitation of the species.

Genome and Metabolome MS-Based Mining of a Marine Strain of Aspergillus affinis

Aspergillus section Circumdati encompasses several species that express both beneficial (e.g., biochemical transformation of steroids and alkaloids, enzymes and metabolites) and harmful compounds (e.g., production of ochratoxin A (OTA)). Given their relevance, it is important to analyze the genetic and metabolic diversity of the species of this section. We sequenced the genome of Aspergillus affinis CMG 70, isolated from sea water, and compared it with the genomes of species from section Circumdati, including A. affinis's strain type. The A. affinis genome was characterized considering secondary metabolites biosynthetic gene clusters (BGCs), carbohydrate-active enzymes (CAZymes), and transporters. To uncover the biosynthetic potential of A. affinis CMG 70, an untargeted metabolomics (LC-MS/MS) approach was used. Cultivating the fungus in the presence and absence of sea salt showed that A. affinis CMG 70 metabolite profiles are salt dependent. Analyses of the methanolic crude extract revealed the presence of both unknown and well-known Aspergillus compounds, such as ochratoxin A, anti-viral (e.g., 3,5-Di-tert-butyl-4-hydroxybenzoic acid and epigallocatechin), anti-bacterial (e.g., 3-Hydroxybenzyl alcohol, l-pyroglutamic acid, lecanoric acid), antifungal (e.g., lpyroglutamic acid, 9,12,13-Trihydroxyoctadec-10-enoic acid, hydroxyferulic acid), and chemotherapeutic (e.g., daunomycinone, mitoxantrone) related metabolites. Comparative analysis of 17 genomes from 16 Aspergillus species revealed abundant CAZymes (568 per species), secondary metabolite BGCs (73 per species), and transporters (1359 per species). Some BGCs are highly conserved in this section (e.g., pyranonigrin E and UNII-YC2Q1O94PT (ACR toxin I)), while others are incomplete or completely lost among species (e.g., bikaverin and chaetoglobosins were found exclusively in series Sclerotiorum, while asperlactone seemed completely lost). The results of this study, including genome analysis and metabolome characterization, emphasize the molecular diversity of A. affinis CMG 70, as well as of other species in the section Circumdati.

MicroSEC filters sequence errors for formalin-fixed and paraffin-embedded samples

The clinical sequencing of tumors is usually performed on formalin-fixed, paraffin-embedded samples and results in many sequencing errors. We identified that most of these errors are detected in chimeric reads caused by single-strand DNA molecules with microhomology. During the end-repair step of library preparation, mutations are introduced by the mis-annealing of two single-strand DNA molecules comprising homologous sequences. The mutated bases are distributed unevenly near the ends in the individual reads. Our filtering pipeline, MicroSEC, focuses on the uneven distribution of mutations in each read and removes the sequencing errors in formalin-fixed, paraffin-embedded samples without over-eliminating the mutations detected also in fresh frozen samples. Amplicon-based sequencing using 97 mutations confirmed that the sensitivity and specificity of MicroSEC were 97% (95% confidence interval: 82–100%) and 96% (95% confidence interval: 88–99%), respectively. Our pipeline will increase the reliability of the clinical sequencing and advance the cancer research using formalin-fixed, paraffin-embedded samples.

Masachika Ikegami and Shinji Kohsaka et al. develop MicroSEC, a computational pipeline to filter sequencing artifacts from archival formalin-fixed and paraffin-embedded samples. Given that archival FFPE tissue is of great interest for genomic analysis, but difficult to reliably analyze, this tool may improve the ability of researchers to probe sequencing data from these samples.

Sequence and phylogenetic analysis of the thrombospondin-related adhesive protein gene of Babesia gibsoni isolates in dogs in South India

Babesia gibsoni, the causative agent of canine piroplasmosis, is a tick-borne intraerythrocytic protozoan parasite predominantly reported in Asian countries. The present study aimed at genotypic characterization of B. gibsoni isolates prevalent in dogs in Kerala, a southern state of India. Blood samples were collected from 272 dogs in Kerala and B. gibsoni infection was detected by microscopy and polymerase chain reaction (PCR). Molecular confirmation of B. gibsoni parasites was carried out by 18S rRNA nested-PCR, followed by sequencing. Nested-PCR detected a higher percentage of dogs (40.44%) positive for B. gibsoni infection than microscopy where 15.81% dogs were detected positive for infection. Genetic characterization of B. gibsoni isolates (n = 11) prevalent in dogs in the state of Kerala was carried out by PCR amplification and sequencing of the 855 bp thrombospondin-related adhesive protein (TRAP) gene fragment. Phylogenetic analysis of the B. gibsoni TRAP (BgTRAP) gene revealed that B. gibsoni isolates from Kerala formed a distinct cluster with the isolates from north India and Bangladesh, away from other East Asian isolates. Nucleotide analysis of the tandem repeats of BgTRAP gene showed considerable genetic variation among Indian isolates that was shared by B. gibsoni isolates of Bangladesh but not by the isolates of East Asian countries. The results of the present study further confirmed that B. gibsoni parasites in a distinct genetic clade are endemic in dogs in India and Bangladesh. However, elaborate studies are required for better understanding of the genetic diversity of B. gibsoni.

Genome-wide characterization of human minisatellite VNTRs: population-specific alleles and gene expression differences

Variable Number Tandem Repeats (VNTRs) are tandem repeat (TR) loci that vary in copy number across a population. Using our program, VNTRseek, we analyzed human whole genome sequencing datasets from 2770 individuals in order to detect minisatellite VNTRs, i.e., those with pattern sizes ≥7 bp. We detected 35 638 VNTR loci and classified 5676 as commonly polymorphic (i.e. with non-reference alleles occurring in >5% of the population). Commonly polymorphic VNTR loci were found to be enriched in genomic regions with regulatory function, i.e. transcription start sites and enhancers. Investigation of the commonly polymorphic VNTRs in the context of population ancestry revealed that 1096 loci contained population-specific alleles and that those could be used to classify individuals into super-populations with near-perfect accuracy. Search for quantitative trait loci (eQTLs), among the VNTRs proximal to genes, indicated that in 187 genes expression differences correlated with VNTR genotype. We validated our predictions in several ways, including experimentally, through the identification of predicted alleles in long reads, and by comparisons showing consistency between sequencing platforms. This study is the most comprehensive analysis of minisatellite VNTRs in the human population to date.

Different types of disease-causing noncoding variants revealed by genomic and gene expression analyses in families with X-linked intellectual disability

The pioneering discovery research of X-linked intellectual disability (XLID) genes has benefitted thousands of individuals worldwide; however, approximately 30% of XLID families still remain unresolved. We postulated that noncoding variants that affect gene regulation or splicing may account for the lack of a genetic diagnosis in some cases. Detecting pathogenic, gene-regulatory variants with the same sensitivity and specificity as structural and coding variants is a major challenge for Mendelian disorders. Here, we describe three pedigrees with suggestive XLID where distinctive phenotypes associated with known genes guided the identification of three different noncoding variants. We used comprehensive structural, single-nucleotide, and repeat expansion analyses of genome sequencing. RNA-Seq from patient-derived cell lines, reverse-transcription polymerase chain reactions, Western blots, and reporter gene assays were used to confirm the functional effect of three fundamentally different classes of pathogenic noncoding variants: a retrotransposon insertion, a novel intronic splice donor, and a canonical splice variant of an untranslated exon. In one family, we excluded a rare coding variant in ARX, a known XLID gene, in favor of a regulatory noncoding variant in OFD1 that correlated with the clinical phenotype. Our results underscore the value of genomic research on unresolved XLID families to aid novel, pathogenic noncoding variant discovery.

A chromosome-level genome assembly for the Pacific oyster Crassostrea gigas

Background

The Pacific oyster (Crassostrea gigas) is a bivalve mollusc with vital roles in coastal ecosystems and aquaculture globally. While extensive genomic tools are available for C. gigas, highly contiguous reference genomes are required to support both fundamental and applied research. Herein we report the creation and annotation of a chromosome-level assembly for C. gigas.

Findings

High-coverage long- and short-read sequence data generated on Pacific Biosciences and Illumina platforms were used to generate an initial assembly, which was then scaffolded into 10 pseudo-chromosomes using both Hi-C sequencing and a high-density linkage map. The assembly has a scaffold N50 of 58.4 Mb and a contig N50 of 1.8 Mb, representing a step advance on the previously published C. gigas assembly. Annotation based on Pacific Biosciences Iso-Seq and Illumina RNA-Seq resulted in identification of ∼30,000 putative protein-coding genes. Annotation of putative repeat elements highlighted an enrichment of Helitron rolling-circle transposable elements, suggesting their potential role in shaping the evolution of the C. gigas genome.

Conclusions

This new chromosome-level assembly will be an enabling resource for genetics and genomics studies to support fundamental insight into bivalve biology, as well as for selective breeding of C. gigas in aquaculture.

A mosquito small RNA genomics resource reveals dynamic evolution and host responses to viruses and transposons

Although mosquitoes are major transmission vectors for pathogenic arboviruses, viral infection has little impact on mosquito health. This immunity is caused in part by mosquito RNA interference (RNAi) pathways that generate antiviral small interfering RNAs (siRNAs) and Piwi-interacting RNAs (piRNAs). RNAi also maintains genome integrity by potently repressing mosquito transposon activity in the germline and soma. However, viral and transposon small RNA regulatory pathways have not been systematically examined together in mosquitoes. Therefore, we developed an integrated mosquito small RNA genomics (MSRG) resource that analyzes the transposon and virus small RNA profiles in mosquito cell cultures and somatic and gonadal tissues across four medically important mosquito species. Our resource captures both somatic and gonadal small RNA expression profiles within mosquito cell cultures, and we report the evolutionary dynamics of a novel Mosquito-Conserved piRNA Cluster Locus (MCpiRCL) made up of satellite DNA repeats. In the larger culicine mosquito genomes we detected highly regular periodicity in piRNA biogenesis patterns coinciding with the expansion of Piwi pathway genes. Finally, our resource enables detection of cross talk between piRNA and siRNA populations in mosquito cells during a response to virus infection. The MSRG resource will aid efforts to dissect and combat the capacity of mosquitoes to tolerate and spread arboviruses.

Dark Matter of Primate Genomes: Satellite DNA Repeats and Their Evolutionary Dynamics

A substantial portion of the primate genome is composed of non-coding regions, so-called "dark matter", which includes an abundance of tandemly repeated sequences called satellite DNA. Collectively known as the satellitome, this genomic component offers exciting evolutionary insights into aspects of primate genome biology that raise new questions and challenge existing paradigms. A complete human reference genome was recently reported with telomere-to-telomere human X chromosome assembly that resolved hundreds of dark regions, encompassing a 3.1 Mb centromeric satellite array that had not been identified previously. With the recent exponential increase in the availability of primate genomes, and the development of modern genomic and bioinformatics tools, extensive growth in our knowledge concerning the structure, function, and evolution of satellite elements is expected. The current state of knowledge on this topic is summarized, highlighting various types of primate-specific satellite repeats to compare their proportions across diverse lineages. Inter- and intraspecific variation of satellite repeats in the primate genome are reviewed. The functional significance of these sequences is discussed by describing how the transcriptional activity of satellite repeats can affect gene expression during different cellular processes. Sex-linked satellites are outlined, together with their respective genomic organization. Mechanisms are proposed whereby satellite repeats might have emerged as novel sequences during different evolutionary phases. Finally, the main challenges that hinder the detection of satellite DNA are outlined and an overview of the latest methodologies to address technological limitations is presented.

Genetic Diversity of the Noncoding Control Region of the Novel Human Polyomaviruses

The genomes of polyomaviruses are characterized by their tripartite organization with an early region, a late region and a noncoding control region (NCCR). The early region encodes proteins involved in replication and transcription of the viral genome, while expression of the late region generates the capsid proteins. Transcription regulatory sequences for expression of the early and late genes, as well as the origin of replication are encompassed in the NCCR. Cell tropism of polyomaviruses not only depends on the appropriate receptors on the host cell, but cell-specific expression of the viral genes is also governed by the NCCR. Thus far, 15 polyomaviruses have been isolated from humans, though it remains to be established whether all of them are genuine human polyomaviruses (HPyVs). The sequences of the NCCR of these HPyVs show high genetic variability and have been best studied in the human polyomaviruses BK and JC. Rearranged NCCRs in BKPyV and JCPyV, the first HPyVs to be discovered approximately 30 years ago, have been associated with the pathogenic properties of these viruses in nephropathy and progressive multifocal leukoencephalopathy, respectively. Since 2007, thirteen novel PyVs have been isolated from humans: KIPyV, WUPyV, MCPyV, HPyV6, HPyV7, TSPyV, HPyV9, HPyV10, STLPyV, HPyV12, NJPyV, LIPyV and QPyV. This review describes all NCCR variants of the new HPyVs that have been reported in the literature and discusses the possible consequences of NCCR diversity in terms of promoter strength, putative transcription factor binding sites and possible association with diseases.

Complete genome analysis of Vibrio mimicus strain SCCF01, a highly virulent isolate from the freshwater catfish

Vibrio mimicus is a foodborne pathogen, which is widely distributed in the aquatic environment. Moreover, it is often involved in aquatic animal diseases. In recent years, V. mimicus is an emerging pathogen in some species of Siluriformes. The strain SCCF01 was isolated from yellow catfish (Pelteobagrus fulvidraco). In this study, we aimed to perform genomic analysis of V. mimicus strain SCCF01 to identify genetic features and evolutionary relationships. Information on gene function and classification was obtained by functional annotation, and circular graph of strain SCCF01 genome, which was created by Circos v0.64. Information on virulence genes (adhesion, flagellum system, exotoxin, and secretory system, etc.) was obtained by virulence genes annotation. Genome element prediction showed that most of the mobile elements were distributed in chromosome I. Therefore, chromosome I of SCCF01 genome has more plasticity than chromosome II and might be larger in size. Genomic linear relationship between the strain of V. mimicus and strain SCCF01 was analyzed by linear pairwise comparison but was unable to determine the relationship. Gene family analysis predicted that the evolutionary direction of strain SCCF01 was: clinical strain → environmental strain → SCCF01 strain. Phylogenetic analysis showed that the strain SCCF01 was more closely related to environmental strains. According to gene family analysis and phylogenetic analysis, we speculated that strain SCCF01 has probably diverged from environmental strains.

Satellite DNA-like repeats are dispersed throughout the genome of the Pacific oyster Crassostrea gigas carried by Helentron non-autonomous mobile elements

Satellite DNAs (satDNAs) are long arrays of tandem repeats typically located in heterochromatin and span the centromeres of eukaryotic chromosomes. Despite the wealth of knowledge about satDNAs, little is known about a fraction of short, satDNA-like arrays dispersed throughout the genome. Our survey of the Pacific oyster Crassostrea gigas sequenced genome revealed genome assembly replete with satDNA-like tandem repeats. We focused on the most abundant arrays, grouped according to sequence similarity into 13 clusters, and explored their flanking sequences. Structural analysis showed that arrays of all 13 clusters represent central repeats of 11 non-autonomous elements named Cg_HINE, which are classified into the Helentron superfamily of DNA transposons. Each of the described elements is formed by a unique combination of flanking sequences and satDNA-like central repeats, coming from one, exceptionally two clusters in a consecutive order. While some of the detected Cg_HINE elements are related according to sequence similarities in flanking and repetitive modules, others evidently arose in independent events. In addition, some of the Cg_HINE's central repeats are related to the classical C. gigas satDNA, interconnecting mobile elements and satDNAs. Genome-wide distribution of Cg_HINE implies non-autonomous Helentrons as a dynamic system prone to efficiently propagate tandem repeats in the C. gigas genome.

Genomics Reveals the Metabolic Potential and Functions in the Redistribution of Dissolved Organic Matter in Marine Environments of the Genus Thalassotalea

Members of the bacterial genus Thalassotalea have been isolated recently from various marine environments, including marine invertebrates. A metagenomic study of the Deepwater Horizon oil plume has identified genes involved in aromatic hydrocarbon degradation in the Thalassotalea genome, shedding light on its potential role in the degradation of crude oils. However, the genomic traits of the genus are not well-characterized, despite the ability of the species to degrade complex natural compounds, such as agar, gelatin, chitin, or starch. Here, we obtained a complete genome of a new member of the genus, designated PS06, isolated from marine sediments containing dead marine benthic macroalgae. Unexpectedly, strain PS06 was unable to grow using most carbohydrates as sole carbon sources, which is consistent with the finding of few ABC transporters in the PS06 genome. A comparative analysis of 12 Thalassotalea genomes provided insights into their metabolic potential (e.g., microaerobic respiration and carbohydrate utilization) and evolutionary stability [including a low abundance of clustered regularly interspaced short palindromic repeats (CRISPR) loci and prophages]. The diversity and frequency of genes encoding extracellular enzymes for carbohydrate metabolism in the 12 genomes suggest that members of Thalassotalea contribute to nutrient cycling by the redistribution of dissolved organic matter in marine environments. Our study improves our understanding of the ecological and genomic properties of the genus Thalassotalea.

Molecular detection and genetic characteristics of Babesia gibsoni in dogs in Shaanxi Province, China

Background

Several members of genus Babesia are important pathogens causing babesiosis in dogs. In China, at least five Babesia species have been described in dogs or ticks. This study sought to determine the prevalence and molecular characteristics of various Babesia spp. in dogs in cities in Shaanxi Province in China, including Xi’an and Hanzhong.

Methods

A total of 371 blood samples were collected from pet dogs presenting to veterinary clinics in the cities of Xi’an and Hanzhong in Shaanxi, China. Babesia spp. DNA was detected via amplification of partial 18S rRNA genes by semi-nested PCR. Almost full-length 18S rRNA, ITS, partial TRAP and complete cytb genes were recovered for analysis of the genetic characteristics and relationships with known isolates.

Results

A single species, Babesia gibsoni, was identified in dogs in Xi’an and Hanzhong. Consistently, B. gibsoni was also detected in 14 ticks collected from positive dogs. Sequence similarities and phylogenetic analysis suggested that the isolates identified herein showed a closer genetic relationship with isolates from East Asian countries rather than India, Bangladesh, or the USA. Sequence analysis based on tandem repeat analysis of the TRAP gene further revealed that specific haplotypes were circulating in both Xi’an and Hanzhong, with no specific regionality. In addition, 10.9% of all isolates with atovaquone (ATV)-resistance were identified because of M121I mutation in the deduced cytb protein.

Conclusions

This study revealed a high prevalence rate of Babesia infection. Babesia gibsoni was the only Babesia species identified in cases of canine babesiosis in the cities of Xi’an and Hanzhong cities in Shaanxi, China. In addition, the TRAP gene presented high genetic diversity across isolates. Such information is useful for elucidating the epidemiological characteristics of canine babesiosis, as well as the overall genetic diversity of Babesia spp. circulating in dog populations in Shaanxi Province.

MSDB: a comprehensive, annotated database of microsatellites

Microsatellites are short tandem repeats of 1–6 nucleotide motifs, studied for their utility as genome markers and in forensics. Recent evidence points to the role of microsatellites in important regulatory functions, and their length polymorphisms at coding regions are linked to various neurodegenerative disorders in humans. Microsatellites show a taxon-specific enrichment in eukaryotic genomes, and their evolution remains poorly understood. Though other databases of microsatellites exist, they fall short on several fronts. MSDB (MicroSatellite DataBase) is a collection of >4 billion microsatellites from 37 680 genomes presented in a user-friendly web portal for easy, interactive analysis and visualization. This is by far the most comprehensive, annotated, updated database to access and analyze microsatellite data of multiple species. The features of MSDB enable users to explore the data as tables that can be filtered and exported, and also as interactive charts to view and compare the data of multiple species simultaneously. Its modularity and architecture permit seamless updates with new data, making it a powerful tool and useful resource to researchers working on this important class of DNA elements, particularly in context of their evolution and emerging roles in genome organization and gene regulation.

The inflated mitochondrial genomes of siphonous green algae reflect processes driving expansion of noncoding DNA and proliferation of introns

Within the siphonous green algal order Bryopsidales, the size and gene arrangement of chloroplast genomes has been examined extensively, while mitochondrial genomes have been mostly overlooked. The recently published mitochondrial genome of Caulerpa lentillifera is large with expanded noncoding DNA, but it remains unclear if this is characteristic of the entire order. Our study aims to evaluate the evolutionary forces shaping organelle genome dynamics in the Bryopsidales based on the C. lentillifera and Ostreobium quekettii mitochondrial genomes. In this study, the mitochondrial genome of O. quekettii was characterised using a combination of long and short read sequencing, and bioinformatic tools for annotation and sequence analyses. We compared the mitochondrial and chloroplast genomes of O. quekettii and C. lentillifera to examine hypotheses related to genome evolution. The O. quekettii mitochondrial genome is the largest green algal mitochondrial genome sequenced (241,739 bp), considerably larger than its chloroplast genome. As with the mtDNA of C. lentillifera, most of this excess size is from the expansion of intergenic DNA and proliferation of introns. Inflated mitochondrial genomes in the Bryopsidales suggest effective population size, recombination and/or mutation rate, influenced by nuclear-encoded proteins, differ between the genomes of mitochondria and chloroplasts, reducing the strength of selection to influence evolution of their mitochondrial genomes.

Tandem repeats lead to sequence assembly errors and impose multi-level challenges for genome and protein databases

The widespread occurrence of repetitive stretches of DNA in genomes of organisms across the tree of life imposes fundamental challenges for sequencing, genome assembly, and automated annotation of genes and proteins. This multi-level problem can lead to errors in genome and protein databases that are often not recognized or acknowledged. As a consequence, end users working with sequences with repetitive regions are faced with 'ready-to-use' deposited data whose trustworthiness is difficult to determine, let alone to quantify. Here, we provide a review of the problems associated with tandem repeat sequences that originate from different stages during the sequencing-assembly-annotation-deposition workflow, and that may proliferate in public database repositories affecting all downstream analyses. As a case study, we provide examples of the Atlantic cod genome, whose sequencing and assembly were hindered by a particularly high prevalence of tandem repeats. We complement this case study with examples from other species, where mis-annotations and sequencing errors have propagated into protein databases. With this review, we aim to raise the awareness level within the community of database users, and alert scientists working in the underlying workflow of database creation that the data they omit or improperly assemble may well contain important biological information valuable to others.

K-mer-Based Motif Analysis in Insect Species across Anopheles, Drosophila, and Glossina Genera and Its Application to Species Classification

Short k-mer sequences from DNA are both conserved and diverged across species owing to their functional significance in speciation, which enables their use in many species classification algorithms. In the present study, we developed a methodology to analyze the DNA k-mers of whole genome, 5' UTR, intron, and 3' UTR regions from 58 insect species belonging to three genera of Diptera that include Anopheles, Drosophila, and Glossina. We developed an improved algorithm to predict and score k-mers based on a scheme that normalizes k-mer scores in different genomic subregions. This algorithm takes advantage of the information content of the whole genome as opposed to other algorithms or studies that analyze only a small group of genes. Our algorithm uses k-mers of lengths 7-9 bp for the whole genome, 5' and 3' UTR regions as well as the intronic regions. Taxonomical relationships based on the whole-genome k-mer signatures showed that species of the three genera clustered together quite visibly. We also improved the scoring and filtering of these k-mers for accurate species identification. The whole-genome k-mer content correlation algorithm showed that species within a single genus correlated tightly with each other as compared to other genera. The genomes of two Aedes and one Culex species were also analyzed to demonstrate how newly sequenced species can be classified using the algorithm. Furthermore, working with several dozen species has enabled us to assign a whole-genome k-mer signature for each of the 58 Dipteran species by making all-to-all pairwise comparison of the k-mer content. These signatures were used to compare the similarity between species and to identify clusters of species displaying similar signatures.

Fluorescence spectroscopic detection and measurement of single telomere molecules

Telomeres are the end-caps of chromosomes that serve to protect the integrity of the genome. Below certain critical lengths, the telomeres can no longer fulfill their protective function, and chromosomal instability ensues. Telomeres shorten during normal cell division due to the end replication problem and are implicated in the development of various aging-associated diseases, including cancer. Telomere length has the potential to serve as a useful biomarker in the field of aging and cancer. However, existing methods of telomere measurement are either too laborious, unable to provide absolute measurement of individual telomere lengths, or limited to certain chromosomes or cell types. Here, we describe an easy single-molecule, fluorescence spectroscopic method for measuring the length of telomeres that permits the profiling of absolute telomere lengths in any DNA sample. We have demonstrated the accurate detection of telomeres as short as 100 bp using cloned telomere standards, and have profiled telomere lengths in human cancer cell lines and primary cells. Since this method allows direct comparison between samples, it could greatly improve the clinical utility of telomere biomarkers.

RB1 Deletion in Retinoblastoma Protein Pathway-Disrupted Cells Results in DNA Damage and Cancer Progression

Proliferative control in cancer cells is frequently disrupted by mutations in the retinoblastoma protein (RB) pathway. Intriguingly, RB1 mutations can arise late in tumorigenesis in cancer cells whose RB pathway is already compromised by another mutation. In this study, we present evidence for increased DNA damage and instability in cancer cells with RB pathway defects when RB1 mutations are induced. We generated isogenic RB1 mutant genotypes with CRISPR/Cas9 in a number of cell lines. Cells with even one mutant copy of RB1 have increased basal levels of DNA damage and increased mitotic errors. Elevated levels of reactive oxygen species as well as impaired homologous recombination repair underlie this DNA damage. When xenografted into immunocompromised mice, RB1 mutant cells exhibit an elevated propensity to seed new tumors in recipient lungs. This study offers evidence that late-arising RB1 mutations can facilitate genome instability and cancer progression that are beyond the preexisting proliferative control deficit.

Estimation of duplication history under a stochastic model for tandem repeats

BACKGROUND

Tandem repeat sequences are common in the genomes of many organisms and are known to cause important phenomena such as gene silencing and rapid morphological changes. Due to the presence of multiple copies of the same pattern in tandem repeats and their high variability, they contain a wealth of information about the mutations that have led to their formation. The ability to extract this information can enhance our understanding of evolutionary mechanisms.

RESULTS

We present a stochastic model for the formation of tandem repeats via tandem duplication and substitution mutations. Based on the analysis of this model, we develop a method for estimating the relative mutation rates of duplications and substitutions, as well as the total number of mutations, in the history of a tandem repeat sequence. We validate our estimation method via Monte Carlo simulation and show that it outperforms the state-of-the-art algorithm for discovering the duplication history. We also apply our method to tandem repeat sequences in the human genome, where it demonstrates the different behaviors of micro- and mini-satellites and can be used to compare mutation rates across chromosomes. It is observed that chromosomes that exhibit the highest mutation activity in tandem repeat regions are the same as those thought to have the highest overall mutation rates. However, unlike previous works that rely on comparing human and chimpanzee genomes to measure mutation rates, the proposed method allows us to find chromosomes with the highest mutation activity based on a single genome, in essence by comparing (approximate) copies of the pattern in tandem repeats.

CONCLUSION

The prevalence of tandem repeats in most organisms and the efficiency of the proposed method enable studying various aspects of the formation of tandem repeats and the surrounding sequences in a wide range of settings.

AVAILABILITY

The implementation of the estimation method is available at http://ips.lab.virginia.edu/smtr .

Tandem Repeats Contribute to Coding Sequence Variation in Bumblebees (Hymenoptera: Apidae)

Tandem repeats (TRs) are highly dynamic regions of the genome. Mutations at these loci represent a significant source of genetic variation and can facilitate rapid adaptation. Bumblebees are important pollinating insects occupying a wide range of habitats. However, to date, molecular mechanisms underlying the potential adaptation of bumblebees to diverse habitats are largely unknown. In the present study, we investigate how TRs contribute to genetic variation in bumblebees, thus potentially facilitating adaptation. We identified 26,595 TRs from the assembled 18 chromosome sequences of the buff-tailed bumblebee (Bombus terrestris), 66.7% of which reside in genic regions. We also compared TRs found in B. terrestris with those present in the assembled genome sequence of a congener, B. impatiens. We found that a total of 1,137 TRs were variable in length between the two sequenced bumblebee species, and further analysis reveals that 101 of them are located within coding regions. These 101 TRs are responsible for coding sequence variation and correspond to protein sequence length variation between the two bumblebee species. The variability of identified TRs in coding regions between bumblebees was confirmed by PCR amplification of a subset of loci. Functional classification of bumblebee genes where coding sequences include variable-length TRs suggests that a majority of genes (87%) that could be assigned to a protein class are related to transcriptional regulation. Our results show that TRs contribute to coding sequence variation in bumblebees, and thus may facilitate the adaptation of bumblebees through diversifying proteins involved in controlling gene expression.

Selective Cytotoxicity of the NAMPT Inhibitor FK866 Toward Gastric Cancer Cells With Markers of the Epithelial-Mesenchymal Transition, Due to Loss of NAPRT

BACKGROUND & AIMS

Markers of the epithelial-to-mesenchymal transition (EMT) in gastric tumor tissues are associated with poor patient outcomes. We performed a screen to identify pharmacologic compounds that kill gastric cancer cells with EMT-associated gene expression patterns and investigate their mechanisms.

METHODS

We identified 29 gastric cancer cell lines with a gene expression signature previously associated with an EMT subtype, based on data from RNA sequence analyses, and confirmed the mesenchymal phenotypes of 7 lines (Hs746T, SNU1750, MKN1, SK4, SNU484, SNU668, and YCC11), based on invasive activity and protein markers. We screened 1,345 compounds for their ability to kill cells with the EMT signature compared with cell lines without this pattern. We tested the effects of identified compounds in BALB/c nude mice bearing GA077 tumors; mice were given intraperitoneal injections of the compound or vehicle (control) twice daily for 24 days and tumor growth was monitored. Proteins associated with the toxicity of the compounds were overexpressed in MKN1 and SNU484 cells or knocked down in MKN45 and SNU719 using small interfering RNAs. We performed immunohistochemical analyses of 942 gastric cancer tissues and investigated associations between EMT markers and protein expression patterns.

RESULTS

The nicotinamide phosphoribosyltransferase inhibitor FK866 killed 6 of 7 gastric cancer cell lines with EMT-associated gene expression signatures but not gastric cancer cells without this signature. The 6 EMT-subtype gastric cell lines expressed significantly low levels of nicotinic acid phosphoribosyltransferase (NAPRT), which makes the cells hypersensitive to nicotinamide phosphoribosyltransferase inhibition. Gastric cell lines that expressed higher levels of NAPRT, regardless of EMT markers, were sensitized to FK866 after knockdown of NAPRT, whereas overexpression of NAPRT in deficient EMT cell lines protected them from FK866-mediated toxicity. Administration of FK866 to nude mice with tumors grown from GA077 cells (human gastric cancer tumors of the EMT subtype) led to tumor regression in 2 weeks; FK866 did not affect tumors grown from MKN45 cells without the EMT expression signature. Loss of NAPRT might promote the EMT, because it stabilizes β-catenin. We correlated the EMT gene expression signature with lower levels of NAPRT in 942 gastric tumors from patients; we also found lower levels of NAPRT mRNA in colorectal, pancreatic, and lung adenocarcinoma tissues with the EMT gene expression signature.

CONCLUSIONS

FK866 selectively kills gastric cancer cells with an EMT gene expression signature by inhibiting nicotinamide phosphoribosyltransferase in cells with NAPRT deficiency. Loss of NAPRT expression, frequently through promoter hypermethylation, is observed in many gastric tumors of the EMT subtype. FK866 might be used to treat patients with tumors of this subtype.

GtTR: Bayesian estimation of absolute tandem repeat copy number using sequence capture and high throughput sequencing

BACKGROUND

Tandem repeats comprise significant proportion of the human genome including coding and regulatory regions. They are highly prone to repeat number variation and nucleotide mutation due to their repetitive and unstable nature, making them a major source of genomic variation between individuals. Despite recent advances in high throughput sequencing, analysis of tandem repeats in the context of complex diseases is still hindered by technical limitations. We report a novel targeted sequencing approach, which allows simultaneous analysis of hundreds of repeats. We developed a Bayesian algorithm, namely - GtTR - which combines information from a reference long-read dataset with a short read counting approach to genotype tandem repeats at population scale. PCR sizing analysis was used for validation.

RESULTS

We used a PacBio long-read sequenced sample to generate a reference tandem repeat genotype dataset with on average 13% absolute deviation from PCR sizing results. Using this reference dataset GtTR generated estimates of VNTR copy number with accuracy within 95% high posterior density (HPD) intervals of 68 and 83% for capture sequence data and 200X WGS data respectively, improving to 87 and 94% with use of a PCR reference. We show that the genotype resolution increases as a function of depth, such that the median 95% HPD interval lies within 25, 14, 12 and 8% of the its midpoint copy number value for 30X, 200X WGS, 395X and 800X capture sequence data respectively. We validated nine targets by PCR sizing analysis and genotype estimates from sequencing results correlated well with PCR results.

CONCLUSIONS

The novel genotyping approach described here presents a new cost-effective method to explore previously unrecognized class of repeat variation in GWAS studies of complex diseases at the population level. Further improvements in accuracy can be obtained by improving accuracy of the reference dataset.

Who Needs This Junk, or Genomic Dark Matter

Centromeres (CEN), pericentromeric regions (periCEN), and subtelomeric regions (subTel) comprise the areas of constitutive heterochromatin (HChr). Tandem repeats (TRs or satellite DNA) are the main components of HChr forming no less than 10% of the mouse and human genome. HChr is assembled within distinct structures in the interphase nuclei of many species - chromocenters. In this review, the main classes of HChr repeat sequences are considered in the order of their number increase in the sequencing reads of the mouse chromocenters (ChrmC). TRs comprise ~70% of ChrmC occupying the first place. Non-LTR (-long terminal repeat) retroposons (mainly LINE, long interspersed nuclear element) are the next (~11%), and endogenous retroviruses (ERV; LTR-containing) are in the third position (~9%). HChr is not enriched with ERV in comparison with the whole genome, but there are differences in distribution of certain elements: while MaLR-like elements (ERV3) are dominant in the whole genome, intracisternal A-particles and corresponding LTR (ERV2) are prevalent in HChr. Most of LINE in ChrmC is represented by the 2-kb fragment at the end of the 2nd open reading frame and its flanking regions. Almost all tandem repeats classified as CEN or periCEN are contained in ChrmC. Our previous classification revealed 60 new mouse TR families with 29 of them being absent in ChrmC, which indicates their location on chromosome arms. TR transcription is necessary for maintenance of heterochromatic status of the HChr genome part. A burst of TR transcription is especially important in embryogenesis and other cases of radical changes in the cell program, including carcinogenesis. The recently discovered mechanism of epigenetic regulation with noncoding sequences transcripts, long noncoding RNA, and its role in embryogenesis and pluripotency maintenance is discussed.

Dynamics of sister chromatid resolution during cell cycle progression

Faithful genome transmission in dividing cells requires that the two copies of each chromosome's DNA package into separate but physically linked sister chromatids. The linkage between sister chromatids is mediated by cohesin, yet where sister chromatids are linked and how they resolve during cell cycle progression has remained unclear. In this study, we investigated sister chromatid organization in live human cells using dCas9-mEGFP labeling of endogenous genomic loci. We detected substantial sister locus separation during G2 phase irrespective of the proximity to cohesin enrichment sites. Almost all sister loci separated within a few hours after their respective replication and then rapidly equilibrated their average distances within dynamic chromatin polymers. Our findings explain why the topology of sister chromatid resolution in G2 largely reflects the DNA replication program. Furthermore, these data suggest that cohesin enrichment sites are not persistent cohesive sites in human cells. Rather, cohesion might occur at variable genomic positions within the cell population.

An intronic VNTR affects splicing of ABCA7 and increases risk of Alzheimer's disease

Mutations leading to premature termination codons in ATP-Binding Cassette Subfamily A Member 7 (ABCA7) are high penetrant risk factors of Alzheimer’s disease (AD). The influence of other genetic variants in ABCA7 and downstream functional mechanisms, however, is poorly understood. To address this knowledge gap, we investigated tandem repetitive regions in ABCA7 in a Belgian cohort of 1529 AD patients and control individuals and identified an intronic variable number tandem repeat (VNTR). We observed strong association between VNTR length and a genome-wide associated signal for AD in the ABCA7 locus. Expanded VNTR alleles were highly enriched in AD patients [odds ratio = 4.5 (1.3–24.2)], and VNTR length inversely correlated with amyloid β_1–42 in cerebrospinal fluid and ABCA7 expression. In addition, we identified three novel ABCA7 alternative splicing events. One isoform in particular—which is formed through exon 19 skipping—lacks the first nucleotide binding domain of ABCA7 and is abundant in brain tissue. We observed a tight correlation between exon 19 skipping and VNTR length. Our findings underline the importance of studying repetitive DNA in complex disorders and expand the contribution of genetic and transcript variation in ABCA7 to AD.

Analysis of microRNA profile of Anopheles sinensis by deep sequencing and bioinformatic approaches

BACKGROUND

microRNAs (miRNAs) are small non-coding RNAs widely identified in many mosquitoes. They are reported to play important roles in development, differentiation and innate immunity. However, miRNAs in Anopheles sinensis, one of the Chinese malaria mosquitoes, remain largely unknown.

METHODS

We investigated the global miRNA expression profile of An. sinensis using Illumina Hiseq 2000 sequencing. Meanwhile, we applied a bioinformatic approach to identify potential miRNAs in An. sinensis. The identified miRNA profiles were compared and analyzed by two approaches. The selected miRNAs from the sequencing result and the bioinformatic approach were confirmed with qRT-PCR. Moreover, target prediction, GO annotation and pathway analysis were carried out to understand the role of miRNAs in An. sinensis.

RESULTS

We identified 49 conserved miRNAs and 12 novel miRNAs by next-generation high-throughput sequencing technology. In contrast, 43 miRNAs were predicted by the bioinformatic approach, of which two were assigned as novel. Comparative analysis of miRNA profiles by two approaches showed that 21 miRNAs were shared between them. Twelve novel miRNAs did not match any known miRNAs of any organism, indicating that they are possibly species-specific. Forty miRNAs were found in many mosquito species, indicating that these miRNAs are evolutionally conserved and may have critical roles in the process of life. Both the selected known and novel miRNAs (asi-miR-281, asi-miR-184, asi-miR-14, asi-miR-nov5, asi-miR-nov4, asi-miR-9383, and asi-miR-2a) could be detected by quantitative real-time PCR (qRT-PCR) in the sequenced sample, and the expression patterns of these miRNAs measured by qRT-PCR were in concordance with the original miRNA sequencing data. The predicted targets for the known and the novel miRNAs covered many important biological roles and pathways indicating the diversity of miRNA functions. We also found 21 conserved miRNAs and eight counterparts of target immune pathway genes in An. sinensis based on the analysis of An. gambiae.

CONCLUSIONS

Our results provide the first lead to the elucidation of the miRNA profile in An. sinensis. Unveiling the roles of mosquito miRNAs will undoubtedly lead to a better understanding of mosquito biology and mosquito-pathogen interactions. This work lays the foundation for the further functional study of An. sinensis miRNAs and will facilitate their application in vector control.

Describing Genomic and Epigenomic Traits Underpinning Emerging Fungal Pathogens

An unprecedented number of pathogenic fungi are emerging and causing disease in animals and plants, putting the resilience of wild and managed ecosystems in jeopardy. While the past decades have seen an increase in the number of pathogenic fungi, they have also seen the birth of new big data technologies and analytical approaches to tackle these emerging pathogens. We review how the linked fields of genomics and epigenomics are transforming our ability to address the challenge of emerging fungal pathogens. We explore the methodologies and bioinformatic toolkits that currently exist to rapidly analyze the genomes of unknown fungi, then discuss how these data can be used to address key questions that shed light on their epidemiology. We show how genomic approaches are leading a revolution into our understanding of emerging fungal diseases and speculate on future approaches that will transform our ability to tackle this increasingly important class of emerging pathogens.