Sequence-based estimation of minisatellite and microsatellite repeat variability

Phenotypic tolerance for rDNA copy number variation within the natural range of C. elegans

The genes for ribosomal RNA (rRNA) are encoded by ribosomal DNA (rDNA), whose structure is notable for being present in arrays of tens to thousands of tandemly repeated copies in eukaryotic genomes. The exact number of rDNA copies per genome is highly variable within a species, with differences between individuals measuring in potentially hundreds of copies and megabases of DNA. The extent to which natural variation in rDNA copy number impacts whole-organism phenotypes such as fitness and lifespan is poorly understood, in part due to difficulties in manipulating such large and repetitive tracts of DNA even in model organisms. Here, we used the natural resource of copy number variation in C. elegans wild isolates to generate new tools and investigated the phenotypic consequences of this variation. Specifically, we generated a panel of recombinant inbred lines (RILs) using a laboratory strain derivative with ∼130 haploid rDNA copies and a wild isolate with ∼417 haploid rDNA copies, one of the highest validated C. elegans rDNA copy number arrays. We find that rDNA copy number is stable in the RILs, rejecting prior hypotheses that predicted copy number instability and copy number reversion. To isolate effects of rDNA copy number on phenotype, we produced a series of near isogenic lines (NILs) with rDNA copy numbers representing the high and low end of the rDNA copy number spectrum in C. elegans wild isolates. We find no correlation between rDNA copy number and phenotypes of rRNA abundance, competitive fitness, early life fertility, lifespan, or global transcriptome under standard laboratory conditions. These findings demonstrate a remarkable ability of C. elegans to tolerate substantial variation in a locus critical to fundamental cell function. Our study provides strain resources for future investigations into the boundaries of this tolerance.

Genome-wide profiling of polymorphic short tandem repeats and their influence on gene expression and trait variation in diverse rice populations

Short tandem repeats (STRs) modulate gene expression and contribute to trait variation. However, a systematic evaluation of the genomic characteristics of STRs has not been conducted, and their influence on gene expression in rice remains unclear. Here, we construct a map of 137,629 polymorphic STRs in the rice (Oryza sativa L.) genome using a population-scale resequencing dataset. A genome-wide survey encompassing 4726 accessions shows that the occurrence frequency, mutational patterns, chromosomal distribution, and functional properties of STRs are correlated with the sequences and lengths of repeat motifs. Leveraging a transcriptome dataset from 127 rice accessions, we identify 44,672 expression STRs (eSTRs) by modeling gene expression in response to the length variation of STRs. These eSTRs are notably enriched in the regulatory regions of genes with active transcriptional signatures. Population analysis identifies numerous STRs that have undergone genetic divergence among different rice groups and 1726 tagged STRs that may be associated with agronomic traits. By editing the (ACT)7 STR in OsFD1 promoter, we further experimentally validate its role in regulating gene expression and phenotype. Our study highlights the contribution of STRs to transcriptional regulation in plants and establishes the foundation for their potential use as alternative targets for genetic improvement.

The landcape of Helicobacter pylori-mediated DNA breaks links bacterial genotoxicity to its oncogenic potential

BACKGROUND

Helicobacter pylori (H. pylori) infection is a significant risk factor for gastric cancer (GC) development. A growing body of evidence suggests a causal link between infection with H. pylori and increased DNA breakage in the host cells. While several mechanisms have been proposed for this damage, their relative impact on the overall bacterial genotoxicity is unknown. Moreover, the link between the formation of DNA damage following infection and the emergence of cancerous structural variants (SV) in the genome of infected cells remained unexplored.

METHODS

We constructed a high-resolution map of genomic H. pylori-induced recurrent break sites using the END-seq method on AGS human gastric cells before and after infection. We next applied END-seq to cycling and arrested cells to identify the role of DNA replication on break formation. Recurrent H. pylori-mediated break sites were further characterized by analyzing published RNA-seq, DRIP-seq, and GRO-seq data at these sites. γH2AX staining and comet assay were used for DNA breakage quantification. Liquid chromatography-mass spectrometry (LC-MS) assay was used to quantify cellular concentrations of dNTPs.

RESULTS

Our data indicated that sites of recurrent H. pylori-mediated DNA breaks are ubiquitous across cell types, localized at replication-related fragile sites, and their breakage is dependent on replication. Consistent with that, we found that H. pylori inflicts nucleotide depletion, and that rescuing the cellular nucleotide pool largely reduced H. pylori-induced DNA breaks. Intriguingly, we found that this genotoxic mechanism operates independently of H. pylori cag pathogenicity island (CagPAI) that encodes for the bacterial type 4 secretion system (T4SS), and its virulence factor, CagA, which was previously implicated in increasing DNA damage by downregulating the DNA damage response. Finally, we show that sites of recurrent H. pylori-mediated breaks coincide with chromosomal deletions observed in patients with intestinal-type GC and that this link potentially elucidates the persistent transcriptional alterations observed in cancer driver genes.

CONCLUSIONS

Our findings indicate that dNTP depletion by H. pylori is a key component of its genotoxicity and suggest a link between H. pylori genotoxicity and its oncogenic potential.

Evolutionary feedback from the environment shapes mechanisms that generate genome variation

Darwin recognized that 'a grand and almost untrodden field of inquiry will be opened, on the causes and laws of variation.' However, because the Modern Synthesis assumes that the intrinsic probability of any individual mutation is unrelated to that mutation's potential adaptive value, attention has been focused on selection rather than on the intrinsic generation of variation. Yet many examples illustrate that the term 'random' mutation, as widely understood, is inaccurate. The probabilities of distinct classes of variation are neither evenly distributed across a genome nor invariant over time, nor unrelated to their potential adaptive value. Because selection acts upon variation, multiple biochemical mechanisms can and have evolved that increase the relative probability of adaptive mutations. In effect, the generation of heritable variation is in a feedback loop with selection, such that those mechanisms that tend to generate variants that survive recurring challenges in the environment would be captured by this survival and thus inherited and accumulated within lineages of genomes. Moreover, because genome variation is affected by a wide range of biochemical processes, genome variation can be regulated. Biochemical mechanisms that sense stress, from lack of nutrients to DNA damage, can increase the probability of specific classes of variation. A deeper understanding of evolution involves attention to the evolution of, and environmental influences upon, the intrinsic variation generated in gametes, in other words upon the biochemical mechanisms that generate variation across generations. These concepts have profound implications for the types of questions that can and should be asked, as omics databases become more comprehensive, detection methods more sensitive, and computation and experimental analyses even more high throughput and thus capable of revealing the intrinsic generation of variation in individual gametes. These concepts also have profound implications for evolutionary theory, which, upon reflection it will be argued, predicts that selection would increase the probability of generating adaptive mutations, in other words, predicts that the ability to evolve itself evolves.

Characterization of genome-wide STR variation in 6487 human genomes

Short tandem repeats (STRs) are abundant and highly mutagenic in the human genome. Many STR loci have been associated with a range of human genetic disorders. However, most population-scale studies on STR variation in humans have focused on European ancestry cohorts or are limited by sequencing depth. Here, we depicted a comprehensive map of 366,013 polymorphic STRs (pSTRs) constructed from 6487 deeply sequenced genomes, comprising 3983 Chinese samples (~31.5x, NyuWa) and 2504 samples from the 1000 Genomes Project (~33.3x, 1KGP). We found that STR mutations were affected by motif length, chromosome context and epigenetic features. We identified 3273 and 1117 pSTRs whose repeat numbers were associated with gene expression and 3'UTR alternative polyadenylation, respectively. We also implemented population analysis, investigated population differentiated signatures, and genotyped 60 known disease-causing STRs. Overall, this study further extends the scale of STR variation in humans and propels our understanding of the semantics of STRs.

The mutational dynamics of short tandem repeats in large, multigenerational families

BACKGROUND

Short tandem repeats (STRs) compose approximately 3% of the genome, and mutations at STR loci have been linked to dozens of human diseases including amyotrophic lateral sclerosis, Friedreich ataxia, Huntington disease, and fragile X syndrome. Improving our understanding of these mutations would increase our knowledge of the mutational dynamics of the genome and may uncover additional loci that contribute to disease. To estimate the genome-wide pattern of mutations at STR loci, we analyze blood-derived whole-genome sequencing data for 544 individuals from 29 three-generation CEPH pedigrees. These pedigrees contain both sets of grandparents, the parents, and an average of 9 grandchildren per family.

RESULTS

We use HipSTR to identify de novo STR mutations in the 2nd generation of these pedigrees and require transmission to the third generation for validation. Analyzing approximately 1.6 million STR loci, we estimate the empirical de novo STR mutation rate to be 5.24 × 10^-5 mutations per locus per generation. Perfect repeats mutate about 2 × more often than imperfect repeats. De novo STRs are significantly enriched in Alu elements.

CONCLUSIONS

Approximately 30% of new STR mutations occur within Alu elements, which compose only 11% of the genome, but only 10% are found in LINE-1 insertions, which compose 17% of the genome. Phasing these mutations to the parent of origin shows that parental transmission biases vary among families. We estimate the average number of de novo genome-wide STR mutations per individual to be approximately 85, which is similar to the average number of observed de novo single nucleotide variants.

Genome-wide survey of D/E repeats in human proteins uncovers their instability and aids in identifying their role in the chromatin regulator ATAD2

D/E repeats are stretches of aspartic and/or glutamic acid residues found in over 150 human proteins. We examined genomic stability of D/E repeats and functional characteristics of D/E repeat-containing proteins vis-à-vis the proteins with poly-Q or poly-A repeats, which are known to undergo pathologic expansions. Mining of tumor sequencing data revealed that D/E repeat-coding regions are similar to those coding poly-Qs and poly-As in increased incidence of trinucleotide insertions/deletions but differ in types and incidence of substitutions. D/E repeat-containing proteins preferentially function in chromatin metabolism and are the more likely to be nuclear and interact with core histones, the longer their repeats are. One of the longest D/E repeats of unknown function is in ATAD2, a bromodomain family ATPase frequently overexpressed in tumors. We demonstrate that D/E repeat deletion in ATAD2 suppresses its binding to nascent and mature chromatin and to the constitutive pericentromeric heterochromatin, where ATAD2 represses satellite transcription.

•

Many human proteins contain runs of aspartic/glutamic acid residues (D/E repeats)

•

D/E repeats show increased incidence of in-frame insertions/deletions in tumors

•

Nuclear and histone-interacting proteins often have long D/E repeats

•

D/E repeat of the oncogene ATAD2 controls its binding to pericentric chromatin

The Chloroplast Genome of the Lichen Photobiont Trebouxiophyceae sp. DW1 and Its Phylogenetic Implications

Lichens are symbiotic associations of algae and fungi. The genetic mechanism of the symbiosis of lichens and the influence of symbiosis on the size and composition of the genomes of symbiotic algae have always been intriguing scientific questions explored by lichenologists. However, there were limited data on lichen genomes. Therefore, we isolated and purified a lichen symbiotic alga to obtain a single strain (Trebouxiophyceae sp. DW1), and then obtained its chloroplast genome information by next-generation sequencing (NGS). The chloroplast genome is 129,447 bp in length, and the GC content is 35.2%. Repetitive sequences with the length of 30–35 bp account for 1.27% of the total chloroplast genome. The simple sequence repeats are all mononucleotide repeats. Codon usage analysis showed that the genome tended to use codon ending in A/U. By comparing the length of different regions of Trebouxiophyceae genomes, we found that the changes in the length of exons, introns, and intergenic sequences affect the size of genomes. Trebouxiophyceae had an unstable chloroplast genome structure, with IRs repeatedly losing during evolution. Phylogenetic analysis showed that Trebouxiophyceae is paraphyletic, and Trebouxiophyceae sp. DW1 is sister to the clade of Koliella longiseta and Pabia signiensis.

The effect of habitat and climatic on microsatellite diversity and allele length variation

Microsatellite loci have been shown to vary according to environment. We studied allelic length variation and diversity in eight microsatellite loci along a sharp climatic and habitat gradient in Israel, using the eastern spadefoot toad (Pelobates syriacus) as our model system. We found a gradual increase in allele lengths from north to south. We used a distance-based redundancy analysis (dbRDA) to associate between allele length and habitat and climatic measures and found that geography and annual climate explained the significant variation in allele length. We also used additional measurements pertaining to demography, heterozygosity and allelic diversity to explore four different hypotheses that might explain the variations in allele length. Our results suggest that the changes we observed in allele lengths may not be purely random but could be influenced by the differential mutation rate and/or local environmental conditions operating at the different locations.

Mononucleotide A-repeats may Play a Regulatory Role in Endothermic Housekeeping Genes

Background:

Coding and non-coding short tandem repeats (STRs) facilitate a great diversity of phenotypic traits. The imbalance of mononucleotide A-repeats around transcription start sites (TSSs) was found in 3 mammals: H. sapiens, M. musculus, and R. norvegicus.

Principal Findings:

We found that the imbalance pattern originated in some vertebrates. A similar pattern was observed in mammals and birds, but not in amphibians and reptiles. We proposed that the enriched A-repeats upstream of TSSs is a novel hallmark of endotherms or warm-blooded animals. Gene ontology analysis indicates that the primary function of upstream A-repeats involves metabolism, cellular transportation, and sensory perception (smell and chemical stimulus) through housekeeping genes.

Conclusions:

Upstream A-repeats may play a regulatory role in the metabolic process of endothermic animals.

Population-Scale Polymorphic Short Tandem Repeat Provides an Alternative Strategy for Allele Mining in Cotton

Short tandem repeats (STRs), which vary in size due to featuring variable numbers of repeat units, are present throughout most eukaryotic genomes. To date, few population-scale studies identifying STRs have been reported for crops. Here, we constructed a high-density polymorphic STR map by investigating polymorphic STRs from 911 Gossypium hirsutum accessions. In total, we identified 556,426 polymorphic STRs with an average length of 21.1 bp, of which 69.08% were biallelic. Moreover, 7,718 (1.39%) were identified in the exons of 6,021 genes, which were significantly enriched in transcription, ribosome biogenesis, and signal transduction. Only 5.88% of those exonic STRs altered open reading frames, of which 97.16% were trinucleotide. An alternative strategy STR-GWAS analysis revealed that 824 STRs were significantly associated with agronomic traits, including 491 novel alleles that undetectable by previous SNP-GWAS methods. For instance, a novel polymorphic STR consisting of GAACCA repeats was identified in GH_D06G1697, with its (GAACCA)₅ allele increasing fiber length by 1.96-4.83% relative to the (GAACCA)₄ allele. The database CottonSTRDB was further developed to facilitate use of STR datasets in breeding programs. Our study provides functional roles for STRs in influencing complex traits, an alternative strategy STR-GWAS for allele mining, and a database serving the cotton community as a valuable resource.

Post-harvest biology and recent advances of storage technologies in sugarcane

Sugarcane deteriorates at a quick rate, just like other perishable crops. The quick loss of sucrose content in sugarcane from the time it is harvested has a significant impact on sugar recovery. This problem of post-harvest sucrose losses in sugarcane is a serious concern in cane-producing countries, as it not only leads to low sugar recovery in mills, but also to poor sugar refining. Unreasonable delays in cane transportation from the fields to the mill are frequently linked to a number of problems related to primary or secondary sucrose losses, all of which contribute to a significant reduction in cane weight and sugar recovery. In sugar mills, the processing of damaged or stale canes also presents a number of challenges, including increased viscosity due to dextran generation, formation of acetic acid, and dextrans due to Leuconostoc spp. invasion, and so on. The combination of all of these variables results in low sugar quality, resulting in significant losses for sugar mills. The primary and secondary losses caused by post-harvest sucrose degradation in sugarcane are enlisted. The employment of physico-chemical technologies in farmers' fields and sugar mills to control and minimize these losses has also been demonstrated.

Short Tandem Repeats as a High-Resolution Marker for Capturing Recent Orangutan Population Evolution

The genus Pongo is ideal to study population genetics adaptation, given its remarkable phenotypic divergence and the highly contrasting environmental conditions it’s been exposed to. Studying its genetic variation bears the promise to reveal a motion picture of these great apes’ evolutionary and adaptive history, and also helps us expand our knowledge of the patterns of adaptation and evolution. In this work, we advance the understanding of the genetic variation among wild orangutans through a genome-wide study of short tandem repeats (STRs). Their elevated mutation rate makes STRs ideal markers for the study of recent evolution within a given population. Current technological and algorithmic advances have rendered their sequencing and discovery more accurate, therefore their potential can be finally leveraged in population genetics studies. To study patterns of population variation within the wild orangutan population, we genotyped the short tandem repeats in a population of 21 individuals spanning four Sumatran and Bornean (sub-) species and eight Southeast Asian regions. We studied the impact of sequencing depth on our ability to genotype STRs and found that the STR copy number changes function as a powerful marker, correctly capturing the demographic history of these populations, even the divergences as recent as 10 Kya. Moreover, gene ontology enrichments for genes close to STR variants are aligned with local adaptations in the two islands. Coupled with more advanced STR-compatible population models, and selection tests, genomic studies based on STRs will be able to reduce the gap caused by the missing heritability for species with recent adaptations.

Length variation in short tandem repeats affects gene expression in natural populations of Arabidopsis thaliana

The genetic basis for the fine-tuned regulation of gene expression is complex and ultimately influences the phenotype and thus the local adaptation of natural populations. Short tandem repeats (STRs) consisting of repetitive DNA motifs have been shown to regulate gene expression. STRs are variable in length within a population and serve as a heritable, but semi-reversible, reservoir of standing genetic variation. For sessile organisms, such as plants, STRs could be of major importance in fine-tuning gene expression as a response to a shifting local environment. Here, we used a transcriptome dataset from natural accessions of Arabidopsis thaliana to investigate population-wide gene expression patterns in light of genome-wide STR variation. We empirically modeled gene expression as a response to the STR length within and around the gene and demonstrated that an association between gene expression and STR length variation is unequivocally present in the sampled population. To support our model, we explored the promoter activity in a transcriptional regulator involved in root hair formation and provided experimentally determined causality between coding sequence length variation and promoter activity. Our results support a general link between gene expression variation and STR length variation in A. thaliana.

The Simple Biology of Flipons and Condensates Enhances the Evolution of Complexity

The classical genetic code maps nucleotide triplets to amino acids. The associated sequence composition is complex, representing many elaborations during evolution of form and function. Other genomic elements code for the expression and processing of RNA transcripts. However, over 50% of the human genome consists of widely dispersed repetitive sequences. Among these are simple sequence repeats (SSRs), representing a class of flipons, that under physiological conditions, form alternative nucleic acid conformations such as Z-DNA, G4 quartets, I-motifs, and triplexes. Proteins that bind in a structure-specific manner enable the seeding of condensates with the potential to regulate a wide range of biological processes. SSRs also encode the low complexity peptide repeats to patch condensates together, increasing the number of combinations possible. In situations where SSRs are transcribed, SSR-specific, single-stranded binding proteins may further impact condensate formation. Jointly, flipons and patches speed evolution by enhancing the functionality of condensates. Here, the focus is on the selection of SSR flipons and peptide patches that solve for survival under a wide range of environmental contexts, generating complexity with simple parts.

Evolutionary pressures and codon bias in low complexity regions of plasmodia

The biological meaning of low complexity regions in the proteins of Plasmodium species is a topic of discussion in evolutionary biology. There is a debate between selectionists and neutralists, who either attribute or do not attribute an effect of low-complexity regions on the fitness of these parasites, respectively. In this work, we comparatively study 22 Plasmodium species to understand whether their low complexity regions undergo a neutral or, rather, a selective and species-dependent evolution. The focus is on the connection between the codon repertoire of the genetic coding sequences and the occurrence of low complexity regions in the corresponding proteins. The first part of the work concerns the correlation between the length of plasmodial proteins and their propensity at embedding low complexity regions. Relative synonymous codon usage, entropy, and other indicators reveal that the incidence of low complexity regions and their codon bias is species-specific and subject to selective evolutionary pressure. We also observed that protein length, a relaxed selective pressure, and a broad repertoire of codons in proteins, are strongly correlated with the occurrence of low complexity regions. Overall, it seems plausible that the codon bias of low-complexity regions contributes to functional innovation and codon bias enhancement of proteins on which Plasmodium species rest as successful evolutionary parasites.

PolyQ length co-evolution in neural proteins

Intermolecular co-evolution optimizes physiological performance in functionally related proteins, ultimately increasing molecular co-adaptation and evolutionary fitness. Polyglutamine (polyQ) repeats, which are over-represented in nervous system-related proteins, are increasingly recognized as length-dependent regulators of protein function and interactions, and their length variation contributes to intraspecific phenotypic variability and interspecific divergence. However, it is unclear whether polyQ repeat lengths evolve independently in each protein or rather co-evolve across functionally related protein pairs and networks, as in an integrated regulatory system. To address this issue, we investigated here the length evolution and co-evolution of polyQ repeats in clusters of functionally related and physically interacting neural proteins in Primates. We observed function-/disease-related polyQ repeat enrichment and evolutionary hypervariability in specific neural protein clusters, particularly in the neurocognitive and neuropsychiatric domains. Notably, these analyses detected extensive patterns of intermolecular polyQ length co-evolution in pairs and clusters of functionally related, physically interacting proteins. Moreover, they revealed both direct and inverse polyQ length co-variation in protein pairs, together with complex patterns of coordinated repeat variation in entire polyQ protein sets. These findings uncover a whole system of co-evolving polyQ repeats in neural proteins with direct implications for understanding polyQ-dependent phenotypic variability, neurocognitive evolution and neuropsychiatric disease pathogenesis.

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.

Comprehensive identification and characterization of simple sequence repeats based on the whole-genome sequences of 14 forest and fruit trees

Simple sequence repeats (SSRs) are popular and important molecular markers that exist widely in plants. Here, we conducted a comprehensive identification and comparative analysis of SSRs in 14 tree species. A total of 16, 298 SSRs were identified from 429, 449 genes, and primers were successfully designed for 99.44% of the identified SSRs. Our analysis indicated that tri-nucleotide SSRs were the most abundant, with an average of ~834 per species. Functional enrichment analysis by combining SSR-containing genes in all species, revealed 50 significantly enriched terms, with most belonging to transcription factor families associated with plant development and abiotic stresses such as Myeloblastosis_DNA-bind_4 (Myb_DNA-bind_4), APETALA2 (AP2), and Fantastic Four meristem regulator (FAF). Further functional enrichment analysis showed that 48 terms related to abiotic stress regulation and floral development were significantly enriched in ten species, whereas no significantly enriched terms were found in four species. Interestingly, the largest number of enriched terms was detected in Citrus sinensis (L.) Osbeck, accounting for 54.17% of all significantly enriched functional terms. Finally, we analyzed AP2 and trihelix gene families (Myb_DNA-bind_4) due to their significant enrichment in SSR-containing genes. The results indicated that whole-genome duplication (WGD) and whole genome triplication (WGT) might have played major roles in the expansion of the AP2 gene family but only slightly affected the expansion of the trihelix gene family during evolution. In conclusion, the identification and comprehensive characterization of SSR markers will greatly facilitate future comparative genomics and functional genomics studies.

Evaluation of a Custom SNP Panel for Identifying and Rectifying of Misjudged Paternity in Deficiency Cases

Parentage testing is routinely performed by genotyping short tandem repeat (STR) through capillary electrophoresis in the present. However, ambiguous or even misjudged paternity based on STRs happens from time to time in cases where only one putative parent is available. We analyzed STR data of 7,818,969 unrelated pairs and 75 close-relative pairs and found that although the probability of a random false match between non-relatives was 4.22 × 10^-6, the incidence of false or ambiguous paternity results between children and first-degree relatives of their true parent was as high as 18.67%. These results highlight the risk of false inclusion of a relative or even non-relatives in parentage testing with STRs. We then validated all ambiguous STR results by targeted sequencing with a custom panel containing 4,830 individual identification single nucleotide polymorphisms (IISNP), found that the ratio of mismatch loci to total SNPs was 1.78-6.95% in close relatives compared with 10.93-13.49% in unrelated pairs. Last, we reported three real cases with undetermined paternity by STRs and rectified them by dissecting with our IISNP panel. These results suggested that high-density IISNP panel can be used to identify and rectify misjudged cases effectively.

Moving beyond disease to function: Physiological roles for polyglutamine-rich sequences in cell decisions

Glutamine-rich tracts, also known as polyQ domains, have received a great deal of attention for their role in multiple neurodegenerative diseases, including Huntington's disease (HD), spinocerebellar ataxia (SCA), and others [22], [27]. Expansions in the normal polyQ tracts are thus commonly linked to disease, but polyQ domains themselves play multiple important functional roles in cells that are being increasingly appreciated. The biochemical nature of these domains allows them to adopt a number of different structures and form large assemblies that enable environmental responsiveness, localized signaling, and cellular memory. In many cases, these involve the formation of condensates that have varied material states. In this review, we highlight known and emerging functional roles for polyQ tracts in normal cell physiology.

Large scale in silico characterization of repeat expansion variation in human genomes

Significant progress has been made in elucidating single nucleotide polymorphism diversity in the human population. However, the majority of the variation space in the genome is structural and remains partially elusive. One form of structural variation is tandem repeats (TRs). Expansion of TRs are responsible for over 40 diseases, but we hypothesize these represent only a fraction of the pathogenic repeat expansions that exist. Here we characterize long or expanded TR variation in 1,115 human genomes as well as a replication cohort of 2,504 genomes, identified using ExpansionHunter Denovo. We found that individual genomes typically harbor several rare, large TRs, generally in non-coding regions of the genome. We noticed that these large TRs are enriched in their proximity to Alu elements. The vast majority of these large TRs seem to be expansions of smaller TRs that are already present in the reference genome. We are providing this TR profile as a resource for comparison to undiagnosed rare disease genomes in order to detect novel disease-causing repeat expansions.

Rhein laden pH-responsive polymeric nanoparticles for treatment of osteoarthritis

Osteoarthritis (OA) is a condition associated with severe inflammation, cartilage destruction and degeneration of joints. Rhein (Rh) is an effective anti-inflammatory drug with proven efficacy in in-vitro and in-vivo models. pH sensitive Rh and NH₄HCO₃ laden poly (lactic-co-glycolic acid (PLGA) nanoparticles (NPs) (Rh-PLGA-NPs@NH₄) are developed for an effective treatment of OA. The Rh-PLGA-NPs@NH₄ are prepared along with Rh-PLGA-NPs as a control by double emulsion method. Rh-PLGA-NPs@NH₄ was characterized for their size, shape, morphology and encapsulation efficiency (EE). The effect of pH on release of Rh from Rh-PLGA-NPs@NH₄ was studied at different pH. Further, the cytotoxicity effect of Rh-PLGA-NPs@NH₄ on THP-1 cells were evaluated. Anti-inflammatory efficacy was evaluated on LPS stimulated THP-1 cells and the release of pro-inflammatory cytokines was evaluated and compared with control. The size of Rh-PLGA-NPs@NH₄ and Rh-PLGA-NPs was found to be 190.7 ± 1.2 nm and 134.6 ± 2.4 nm respectively with poly dispersity (PDI) 0.14 and 0.15. The zeta potential of Rh-PLGA-NPs@NH₄ was found to be -22 ± 1.12 mV. Rh-PLGA-NPs@NH₄ were uniform, smooth and spherical shape as confirmed using electron microscopy analysis. Rh-PLGA-NPs@NH₄ release the Rh more effectively in the low pH of synovial fluid environment (SFE). Rh-PLGA-NPs@NH₄ also significantly affect inflammatory cytokines TNF-α and IL-1β and reduced their release in LPS stimulated THP-1 cells. Reactive oxygen species (ROS), a mediator responsible for the cartilage collapse was also found to be reduced. Results proposes that Rh-PLGA-NPs could provide therapeutic solution to those patients who suffer from chronic joint ailments by reducing the progression of OA.

Relatively semi-conservative replication and a folded slippage model for short tandem repeats

Background

The ubiquitous presence of short tandem repeats (STRs) in virtually all genomes implicates their functional relevance, while a widely-accepted definition of STR is yet to be established. Previous studies majorly focus on relatively longer STRs, while shorter repeats were generally excluded. Herein, we have adopted a more generous criteria to define shorter repeats, which has led to the definition of a much larger number of STRs that lack prior analysis. Using this definition, we analyzed the short repeats in 55 randomly selected segments in 55 randomly selected genomic sequences from a fairly wide range of species covering animals, plants, fungi, protozoa, bacteria, archaea and viruses.

Results

Our analysis reveals a high percentage of short repeats in all 55 randomly selected segments, indicating that the universal presence of high-content short repeats could be a common characteristic of genomes across all biological kingdoms. Therefore, it is reasonable to assume a mechanism for continuous production of repeats that can make the replicating process relatively semi-conservative. We have proposed a folded replication slippage model that considers the geometric space of nucleotides and hydrogen bond stability to explain the mechanism more explicitly, with improving the existing straight-line slippage model. The folded slippage model can explain the expansion and contraction of mono- to hexa- nucleotide repeats with proper folding angles. Analysis of external forces in the folding template strands also suggests that expansion exists more commonly than contraction in the short tandem repeats.

Conclusion

The folded replication slippage model provides a reasonable explanation for the continuous occurrences of simple sequence repeats in genomes. This model also contributes to the explanation of STR-to-genome evolution and is an alternative model that complements semi-conservative replication.

Intragenic repeat expansion in the cell wall protein gene HPF1 controls yeast chronological aging

Aging varies among individuals due to both genetics and environment, but the underlying molecular mechanisms remain largely unknown. Using a highly recombined Saccharomyces cerevisiae population, we found 30 distinct quantitative trait loci (QTLs) that control chronological life span (CLS) in calorie-rich and calorie-restricted environments and under rapamycin exposure. Calorie restriction and rapamycin extended life span in virtually all genotypes but through different genetic variants. We tracked the two major QTLs to the cell wall glycoprotein genes FLO11 and HPF1 We found that massive expansion of intragenic tandem repeats within the N-terminal domain of HPF1 was sufficient to cause pronounced life span shortening. Life span impairment by HPF1 was buffered by rapamycin but not by calorie restriction. The HPF1 repeat expansion shifted yeast cells from a sedentary to a buoyant state, thereby increasing their exposure to surrounding oxygen. The higher oxygenation altered methionine, lipid, and purine metabolism, and inhibited quiescence, which explains the life span shortening. We conclude that fast-evolving intragenic repeat expansions can fundamentally change the relationship between cells and their environment with profound effects on cellular lifestyle and longevity.

Polymorphism of simple sequence repeats may quantitatively regulate gene transcription

The degree of polymorphism, i.e., DNA sequence divergence, of short AT-rich tandemly arranged simple sequence repeats at or near promoters and 5'- untranslated regions of mRNA may quantitatively regulate transcription of tissue-specific genes. Less polymorphic repeats allow greater gene expression. Preferential binding of hypophosphorylated H1 histone to these repeats may diminish binding of transcription factors. Preferential binding of hypophosphorylated high mobility group chromatin proteins would increase this binding. Shorter simple sequence repeats have undergone fewer point mutations than longer repeats, hence they are less polymorphic and more conserved. The role of transcribed simple sequence repeats in frog embryo germ layer determination is considered.

Genotype-phenotype correlation and description of two novel mutations in Iranian patients with glycogen storage disease 1b (GSD1b)

BACKGROUND

Glycogen storage disease (GSD) is a rare inborn error of the synthesis or degradation of glycogen metabolism. GSD1, the most common type of GSD, is categorized into GSD1a and GSD1b which caused by the deficiency of glucose-6-phosphatase (G6PC) and glucose-6-phosphate transporter (SLC37A4), respectively. The high rates of consanguineous marriages in Iran provide a desirable context to facilitate finding the homozygous pathogenic mutations. This study designates to evaluate the clinical and genetic characteristics of patients with GSD1b to assess the possible genotype-phenotype correlation.

RESULTS

Autozygosity mapping was performed on nineteen GSD suspected families to suggest the causative loci. The mapping was done using two panels of short tandem repeat (STR) markers linked to the corresponding genes. The patients with autozygous haplotype block for the markers flanking the genes were selected for direct sequencing. Six patients showed autozygosity in the candidate markers for SLC37A4. Three causative variants were detected. The recurrent mutation of c.1042_1043delCT (p.Leu348Valfs*53) and a novel missense mutation of c.365G > A (p.G122E) in the homozygous state were identified in the SLC37A4. In silico analysis was performed to predict the pathogenicity of the variants. A novel whole SLC37A4 gene deletion using long-range PCR and sequencing was confirmed as well. Severe and moderate neutropenia was observed in patients with frameshift and missense variants, respectively. The sibling with the whole gene deletion has shown both severe neutropenia and leukopenia.

CONCLUSIONS

The results showed that the hematological findings may have an appropriate correlation with the genotype findings. However, for a definite genotype-phenotype correlation, specifically for the clinical and biochemical phenotype, further studies with larger sample sizes are needed.

Mutational spectrum of autosomal recessive limb-girdle muscular dystrophies in a cohort of 112 Iranian patients and reporting of a possible founder effect

Background

Limb-girdle muscular dystrophies are a group of genetically heterogeneous diseases that are inherited in both autosomal dominant (LGMDD) and autosomal recessive forms (LGMDR), the latter is more common especially in populations with high consanguineous marriages like Iran. In the present study, we aimed to investigate the genetic basis of patients who are suspicious of being affected by LGMDR.

DNA samples of 60 families suspected of LGMD were extracted from their whole blood. Four short tandem repeat (STR) markers for each candidate genes related to LGMD R1 (calpain3 related)- R6 (δ-sarcoglycan-related) were selected, and all these 24 STRs were applied in two sets of multiplex PCR. After autozygosity mapping, Sanger sequencing and variant analysis were done. Predicting identified variants’ effect was performed using in-silico tools, and they were interpreted according to the American College of Medical Genomics and Genetics (ACMG) guideline. MLPA was used for those patients who had large deletions.

Fresh muscle specimens were taken from subjects and were evaluated using the conventional panel of histochemical stains.

Results

forty out of sixty families showed homozygote haplotypes in CAPN3, DYSF, SGCA, and SGCB genes. The exons and intron-exon boundaries of the relevant genes were sequenced and totally 38 mutations including CAPN3 (n = 15), DYSF (n = 9), SGCB (n = 11), and SGCA (n = 3) were identified. Five out of them were novel. The most prevalent form of LGMDs in our study was calpainopathy followed by sarcoglycanopathy in which beta-sarcoglycanopathy was the most common form amongst them. Exon 2 deletion in the SGCB gene was the most frequent mutation in this study.

We also reported evidence of a possible founder effect in families with mutations in DYSF and SGCB genes. We also detected a large consanguineous family suffered from calpainopathy who showed allelic heterogeneity.

Conclusions

This study can expand our knowledge about the genetic spectrum of LGMD in Iran, and also suggest the probable founder effects in some Iranian subpopulations which confirming it with more sample size can facilitate our genetic diagnosis and genetic counseling.

Development and comprehensive evaluation of a noninvasive prenatal paternity testing method through a scaled trial

PURPOSE

To eliminate the miscarriage risks caused by traditional invasive sampling methods, we develop a noninvasive prenatal paternity testing (NIPPT) method and evaluate its efficiency, reliability and sensitivity based on a scaled trial.

METHODS

We use maternal cell-free DNA and massive parallel sequencing to obtain NIPPT genotypes for parents and fetuses based on quality-controlled genome-wide single nucleotide polymorphisms (SNPs). In a preliminary testing, data from 14 pregnant women and 7 negative controls are used for setting threshold of fetal genotyping in reference to postpartum children. After that, those from 349 cases with pregnancies of 6-35 gestational weeks (GW) and 9 negative controls from non-pregnant women who have fertility experience previously are in-depth evaluated.

RESULTS

In all cases, the biological fathers have been successfully identified from unrelated with a combined paternity index (CPI) of 3.58 × 10¹⁸ - 1.46 × 10¹⁶⁵ for the cases versus 1.52 × 10^-22 - 2.30 × 10^-839 for the controls. For negative controls, fetal SNPs originating from previous pregnancies could not be detected. Our NIPPT results completely aligned with the invasive prenatal test results using PCR-CE STR methods.

CONCLUSION

NIPPT can be applied to determine paternity accurately from 6 weeks after conception until birth and may serve as an alternative prenatal paternity test advantageous to the currently-used methods.

De novo emergence and potential function of human-specific tandem repeats in brain-related loci

Tandem repeats (TRs) are widespread in the genomes of all living organisms. In eukaryotes, they are found in both coding and noncoding regions and have potential roles in the regulation of cellular processes such as transcription, translation and in the modification of protein structure. Recent studies have highlighted TRs as a key regulator of gene expression and a potential contributor to human evolution. Thus, TRs are emerging as an important source of variation that can result in differential gene expression at intra- and inter-species levels. In this study, we performed a genome-wide survey to identify TRs that have emerged in the human lineage. We further examined these loci to explore their potential functional significance for human evolution. We identified 152 human-specific TR (HSTR) loci containing a repeat unit of more than ten bases, with most of them showing a repeat count of two. Gene set enrichment analysis showed that HSTR-associated genes were associated with biological functions in brain development and synapse function. In addition, we compared gene expression of human HSTR loci with orthologues from non-human primates (NHP) in seven different tissues. Strikingly, the expression level of HSTR-associated genes in brain tissues was significantly higher in human than in NHP. These results suggest the possibility that de novo emergence of TRs could have resulted in altered gene expression in humans within a short-time frame and contributed to the rapid evolution of human brain function.

Characterization and haplotype study of 6 novel STR markers related to the KCNQ1 gene in heterogeneous cardiovascular disorders in the Iranian population

The KCNQ1 gene has a significant role in long QT syndrome, Jervell and Lange-Nielsen syndrome, familial atrial fibrillation, and short QT syndrome. Analyzing such heterogeneous disorders, six novel short tandem repeat (STR) markers around the KCNQ1 gene were found and evaluated in a healthy population, and other statistical traits of the markers were detected.Materials and methods: Using Tandem Repeats Finder (TRF) and Sequence-Based Estimation of Repeat Variability (SERV) software, STR markers were detected with valid tetra- and pentanucleotide repeats. The markers were investigated for a total of 60 unrelated Iranian healthy individuals and analyzed using GenAlEx 6.502 and Cervus 3.0.7.Results: A total of 77 haplotypes was detected, of which 25 haplotypes were unique and the others occurred at least two times. The number of haplotypes per locus ranged from 7 to 18 with the highest frequency of 69.2%, and the observed heterozygosity was calculated as 0.589. The power of discrimination ranged from 0.70 to 0.96. Five of the markers meet Hardy–Weinberg equilibrium.Conclusion: A novel panel of STR markers was described with high allele heterozygosity and segregation in every locus, which may lead to faster and more credible recognition of the disease-inducing KCNQ1 gene and allow it to be used for human identity testing and prenatal diagnosis.

Microsatellite Stability in STR Analysis Aspergillus fumigatus Depends on Number of Repeat Units

More than a decade ago a short tandem repeat-based typing method was developed for the fungus Aspergillus fumigatus. This STRAf assay is based on the analysis of nine short tandem repeat markers. Interpretation of this STRAf assay is complicated when there are only one or two differences in tandem repeat markers between isolates, as the stability of these markers is unknown. To determine the stability of these nine markers, a STRAf assay was performed on 73–100 successive generations of five clonally expanded A. fumigatus isolates. In a total of 473 generations we found five times an increase of one tandem repeat unit. Three changes were found in the trinucleotide repeat marker STRAf 3A, while the other two were found in the trinucleotide repeat marker STRAf 3C. The di- or tetranucleotide repeats were not altered. The altered STRAf markers 3A and 3C demonstrated the highest number of repeat units (≥50) as compared to the other markers (≤26). Altogether, we demonstrated that 7 of 9 STRAf markers remain stable for 473 generations and that the frequency of alterations in tandem repeats is positively correlated with the number of repeats. The potential low level instability of STRAf markers 3A and 3C should be taken into account when interpreting STRAf data during an outbreak.

Targeting high transcriptional control activity of long mononucleotide A-T repeats in cancer by Argonaute 1

Epigenetic regulatory changes alter the gene regulation function of DNA repeat elements in cancer and consequently promote malignant phenotypes. Some short tandem repeat sequences, distributed throughout the human genome, can play a role as cis-regulatory elements of the genes. Distributions of tandem long (≥10) and short (<10) A-T repeats in the genome are different depending on gene functions. Long repeats are more commonly found in housekeeping genes and may regulate genes in harmonious fashion. Mononucleotide A-repeats around transcription start sites interact with Argonaute proteins (AGOs) to regulate gene expression. miRNA-bound AGO alterations in cancer have been reported; consequently, these changes would affect genes containing mononucleotide A- and T-repeats. Here, we showed an unprecedented hallmark of gene regulation in cancer. We evaluated the gene expression profiles reported in the Gene Expression Omnibus and found a high density of 13-27 A-T repeats in the up-regulated genes in malignancies derived from the bladder, cervix, head and neck, ovary, vulva, breast, colon, liver, lung, prostate, kidney, thyroid, adrenal gland, bone, blood cells, muscle and brain. Transfection of cell-penetrating protein tag AGO1 containing poly uracils (CPP-AGO1-polyUs) to the lung cancer cell lines altered gene regulation depending on the presence of long A-T repeats. CPP-AGO1-polyUs limited cell proliferation and the ability of a cancer cell to grow into a colony in lung cancer cell lines. In conclusion, long A-T repeats up-regulated many genes in cancer that can be targeted by AGO1 to change the expression of many genes and limited cancer growth.

Mapping the landscape of tandem repeat variability by targeted long read single molecule sequencing in familial X-linked intellectual disability

BACKGROUND

The etiology of more than half of all patients with X-linked intellectual disability remains elusive, despite array-based comparative genomic hybridization, whole exome or genome sequencing. Since short read massive parallel sequencing approaches do not allow the detection of larger tandem repeat expansions, we hypothesized that such expansions could be a hidden cause of X-linked intellectual disability.

METHODS

We selectively captured over 1800 tandem repeats on the X chromosome and characterized them by long read single molecule sequencing in 3 families with idiopathic X-linked intellectual disability.

RESULTS

In male DNA samples, full tandem repeat length sequences were obtained for 88-93% of the targets and up to 99.6% of the repeats with a moderate guanine-cytosine content. Read length and analysis pipeline allow to detect cases of > 900 bp tandem repeat expansion. In one family, one repeat expansion co-occurs with down-regulation of the neighboring MIR222 gene. This gene has previously been implicated in intellectual disability and is apparently linked to FMR1 and NEFH overexpression associated with neurological disorders.

CONCLUSIONS

This study demonstrates the power of single molecule sequencing to measure tandem repeat lengths and detect expansions, and suggests that tandem repeat mutations may be a hidden cause of X-linked intellectual disability.

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.

Molecular genetics of a cohort of 635 cases of phenylketonuria in a consanguineous population

Phenylketonuria (PKU) is an inborn error of amino acid metabolism caused by mutations in the phenylalanine hydroxylase (PAH) gene, characterized by intellectual deficit and neuropsychiatric complications in untreated patients with estimated frequency of about one in 10,000 to 15,000 live births. PAH deficiency can be detected by neonatal screening in nearly all cases with hyperphenylalaninemia on a heel prick blood spot. Molecular testing of the PAH gene can then be performed in affected family members. Herein, we report molecular study of 635 patients genetically diagnosed with PKU from all ethnicities in Iran. The disease-causing mutations were found in 611 (96.22%) of cases. To the best of our knowledge, this is the most comprehensive molecular genetics study of PKU in Iran, identifying 100 distinct mutations in the PAH gene, including 15 previously unreported mutations. Interestingly, we found unique cases of PKU with uniparental disomy, germline mosaicism, and coinheritance with another Mendelian single-gene disorder that provides new insights for improving the genetic counseling, prenatal diagnosis (PND), and/or pre-implantation genetic diagnosis (PGD) for the inborn error of metabolism group of disorders.

Development of 17 novel microsatellite markers for Lycoris aurea and L. radiata (Amaryllidaceae) using next-generation sequencing

PREMISE OF THE STUDY

Lycoris is an ornamental and medicinal plant. We developed microsatellite markers for L. aurea and L. radiata simultaneously by using a hybrid between these two species.

METHODS AND RESULTS

Ion Torrent next-generation sequencing produced 1,784,504 reads. Testing 64 primer sets allowed for the development of 17 novel microsatellite markers: 16 for L. aurea, 10 for L. radiata, and nine common markers. Lycoris aurea had one to 12 alleles per locus and observed and expected heterozygosity levels of 0-0.923 and 0.038-0.809, respectively. Lycoris radiata had three to 12 alleles per locus and observed and expected heterozygosity levels of 0-0.909 and 0.127-0.797, respectively. Ten markers were cross-amplified for L. sprengeri.

CONCLUSIONS

Hybrid sequencing can facilitate the cost-effective development of molecular markers for parental species. The markers developed here are useful for studying Lycoris population structure.

Massive variation of short tandem repeats with functional consequences across strains of Arabidopsis thaliana

Short tandem repeat (STR) mutations may comprise more than half of the mutations in eukaryotic coding DNA, yet STR variation is rarely examined as a contributor to complex traits. We assessed this contribution across a collection of 96 strains of Arabidopsis thaliana, genotyping 2046 STR loci each, using highly parallel STR sequencing with molecular inversion probes. We found that 95% of examined STRs are polymorphic, with a median of six alleles per STR across these strains. STR expansions (large copy number increases) are found in most strains, several of which have evident functional effects. These include three of six intronic STR expansions we found to be associated with intron retention. Coding STRs were depleted of variation relative to noncoding STRs, and we detected a total of 56 coding STRs (11%) showing low variation consistent with the action of purifying selection. In contrast, some STRs show hypervariable patterns consistent with diversifying selection. Finally, we detected 133 novel STR-phenotype associations under stringent criteria, most of which could not be detected with SNPs alone, and validated some with follow-up experiments. Our results support the conclusion that STRs constitute a large, unascertained reservoir of functionally relevant genomic variation.

Biological Roles of Protein-Coding Tandem Repeats in the Yeast Candida Albicans

Tandem repeat (TR) DNA mutates faster than other DNA by insertion and deletion of repeats. Large parts of eukaryotic proteomes are encoded by ORFs containing protein-coding TRs (TR-ORFs, pcTRs) with largely unknown biological consequences. We explored these in the yeast Candida albicans, an opportunistic human pathogen. We found that almost half of C. albicans’ proteins are encoded by TR-ORFs. pcTR frequency differed only moderately between different gene (GO) categories. Bioinformatic predictions of genome-wide mutation rates and clade-specific differences in pcTR allele frequencies indicated that pcTRs (i) significantly increase the genome-wide mutation rate; (ii) significantly impact on fitness and (iii) allow the evolution of selectively advantageous clade-specific protein variants. Synonymous mutations reduced the repetitiveness of many amino acid repeat-encoding pcTRs. A survey, in 58 strains, revealed that in some pcTR regions in which repetitiveness was not significantly diminished by synonymous mutations the habitat predicted which alleles were present, suggesting roles of pcTR mutation in short-term adaptation and pathogenesis. In C. albicans pcTR mutation apparently is an important mechanism for mutational advance and possibly also rapid adaptation, with synonymous mutations providing a mechanism for adjusting mutation rates of individual pcTRs. Analyses of Arabidopsis and human pcTRs showed that the latter also occurs in other eukaryotes.

Genomic architecture of haddock (Melanogrammus aeglefinus) shows expansions of innate immune genes and short tandem repeats

Background

Increased availability of genome assemblies for non-model organisms has resulted in invaluable biological and genomic insight into numerous vertebrates, including teleosts. Sequencing of the Atlantic cod (Gadus morhua) genome and the genomes of many of its relatives (Gadiformes) demonstrated a shared loss of the major histocompatibility complex (MHC) II genes 100 million years ago. An improved version of the Atlantic cod genome assembly shows an extreme density of tandem repeats compared to other vertebrate genome assemblies. Highly contiguous assemblies are therefore needed to further investigate the unusual immune system of the Gadiformes, and whether the high density of tandem repeats found in Atlantic cod is a shared trait in this group.

Results

Here, we have sequenced and assembled the genome of haddock (Melanogrammus aeglefinus) – a relative of Atlantic cod – using a combination of PacBio and Illumina reads. Comparative analyses reveal that the haddock genome contains an even higher density of tandem repeats outside and within protein coding sequences than Atlantic cod. Further, both species show an elevated number of tandem repeats in genes mainly involved in signal transduction compared to other teleosts. A characterization of the immune gene repertoire demonstrates a substantial expansion of MCHI in Atlantic cod compared to haddock. In contrast, the Toll-like receptors show a similar pattern of gene losses and expansions. For the NOD-like receptors (NLRs), another gene family associated with the innate immune system, we find a large expansion common to all teleosts, with possible lineage-specific expansions in zebrafish, stickleback and the codfishes.

Conclusions

The generation of a highly contiguous genome assembly of haddock revealed that the high density of short tandem repeats as well as expanded immune gene families is not unique to Atlantic cod – but possibly a feature common to all, or most, codfishes. A shared expansion of NLR genes in teleosts suggests that the NLRs have a more substantial role in the innate immunity of teleosts than other vertebrates. Moreover, we find that high copy number genes combined with variable genome assembly qualities may impede complete characterization of these genes, i.e. the number of NLRs in different teleost species might be underestimates.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4616-y) contains supplementary material, which is available to authorized users.

In silico characterization of tandem repeats in Trichophyton rubrum and related dermatophytes provides new insights into their role in pathogenesis

Trichophyton rubrum is the most common etiological agent of dermatophytoses worldwide, which is able to degrade keratinized tissues. The sequencing of the genome of different dermatophyte species has provided a large amount of data, including tandem repeats that may play a role in genetic variability and in the pathogenesis of these fungi. Tandem repeats are adjacent DNA sequences of 2–200 nucleotides in length, which exert regulatory and adaptive functions. These repetitive DNA sequences are found in different classes of fungal proteins, especially those involved in cell adhesion, a determinant factor for the establishment of fungal infection. The objective of this study was to develop a Dermatophyte Tandem Repeat Database (DTRDB) for the storage and identification of tandem repeats in T. rubrum and six other dermatophyte species. The current version of the database contains 35 577 tandem repeats detected in 16 173 coding sequences. The repeats can be searched using entry parameters such as repeat unit length (nt—nucleotide), repeat number, variability score, and repeat sequence motif. These data were used to study the relative frequency and distribution of repeats in the sequences, as well as their possible functions in dermatophytes. A search of the database revealed that these repeats occur in 22–33% of genes transcribed in dermatophytes where they could be involved in the success of adaptation to the host tissue and establishment of infection. The repeats were detected in transcripts that are mainly related to three biological processes: regulation, adhesion, and metabolism. The database developed enables users to identify and analyse tandem repeat regions in target genes related to pathogenicity and fungal–host interactions in dermatophytes and may contribute to the discovery of new targets for the development of antifungal agents.

Database URL:http://comp.mch.ifsuldeminas.edu.br/dtrdb/

Systematic Profiling of Short Tandem Repeats in the Cattle Genome

Short tandem repeats (STRs), or microsatellites, are genetic variants with repetitive 2–6 base pair motifs in many mammalian genomes. Using high-throughput sequencing and experimental validations, we systematically profiled STRs in five Holsteins. We identified a total of 60,106 microsatellites and generated the first high-resolution STR map, representing a substantial pool of polymorphism in dairy cattle. We observed significant STRs overlap with functional genes and quantitative trait loci (QTL). We performed evolutionary and population genetic analyses using over 20,000 common dinucleotide STRs. Besides corroborating the well-established positive correlation between allele size and variance in allele size, these analyses also identified dozens of outlier STRs based on two anomalous relationships that counter expected characteristics of neutral evolution. And one STR locus overlaps with a significant region of a summary statistic designed to detect STR-related selection. Additionally, our results showed that only 57.1% of STRs located within SNP-based linkage disequilibrium (LD) blocks whereas the other 42.9% were out of blocks. Therefore, a substantial number of STRs are not tagged by SNPs in the cattle genome, likely due to STR's distinct mutation mechanism and elevated polymorphism. This study provides the foundation for future STR-based studies of cattle genome evolution and selection.

Single nucleotide polymorphisms and microsatellites in the canine glutathione S-transferase pi 1 (GSTP1) gene promoter

BACKGROUND

Genetic polymorphisms within the glutathione S-transferase P1 (GSTP1) gene affect the elimination of toxic xenobiotics by the GSTP1 enzyme. In dogs, exposure to environmental chemicals that may be GSTP1 substrates is associated with cancer. The objectives of this study were to investigate the genetic variability in the GSTP1 promoter in a diverse population of 278 purebred dogs, compare the incidence of any variants found between breeds, and predict their effects on gene expression. To provide information on ancestral alleles, a number of wolves, coyotes, and foxes were also sequenced.

RESULTS

Fifteen single nucleotide polymorphisms (SNPs) and two microsatellites were discovered. Three of these loci were only polymorphic in dogs while three other SNPs were unique to wolves and coyotes. The major allele at c.-46 is T in dogs but is C in the wild canids. The c.-185 delT variant was unique to dogs. The microsatellite located in the 5' untranslated region (5'UTR) was a highly polymorphic GCC tandem repeat, consisting of simple and compound alleles that varied in size from 10 to 22-repeat units. The most common alleles consisted of 11, 16, and 17-repeats. The 11-repeat allele was found in 10% of dogs but not in the other canids. Unequal recombination and replication slippage between similar and distinct alleles may be the mechanism for the multiple microsatellites observed. Twenty-eight haplotypes were constructed in the dog, and an additional 8 were observed in wolves and coyotes. While the most common haplotype acrossbreeds was the wild-type *1A(17), other prevalent haplotypes included *3A(11) in Greyhounds, *6A(16) in Labrador Retrievers, *9A(16) in Golden Retrievers, and *8A(19) in Standard Poodles. Boxers and Siberian Huskies exhibited minimal haplotypic diversity. Compared to the simple 16*1 allele, the compound 16*2 allele (found in 12% of dogs) may interfere with transcription factor binding and/or the stability of the GSTP1 transcript.

CONCLUSIONS

Dogs and other canids exhibit extensive variation in the GSTP1 promoter. Genetic polymorphisms within distinct haplotypes prevalent in certain breeds can affect GSTP1 expression and carcinogen detoxification, and thus may be useful as genetic markers for cancer in dogs.

Long noncoding RNAs in the model species Brachypodium distachyon

Eukaryotic genomes are pervasively transcribed and only a small portion of the transcribed sequences belongs to protein coding genes. High-throughput sequencing technology contributed to consolidate this perspective, allowing the identification of numerous noncoding RNAs with key roles in biological processes. Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nt with limited phylogenetic conservation, expressed at low levels and characterized by tissue/organ specific expression profiles. Although a large set of lncRNAs has been identified, the functional roles of lncRNAs are only beginning to be recognized and the molecular mechanism of lncRNA-mediated gene regulation remains largely unexplored, particularly in plants where their annotation and characterization are still incomplete. Using public and proprietary poly-(A)+ RNA-seq data as well as a collection of full length ESTs from several organs, developmental stages and stress conditions in three Brachypodium distachyon inbred lines, we describe the identification and the main features of thousands lncRNAs. Here we provide a genome-wide characterization of lncRNAs, highlighting their intraspecies conservation and describing their expression patterns among several organs/tissues and stress conditions. This work represents a fundamental resource to deepen our knowledge on long noncoding RNAs in C3 cereals, allowing the Brachypodium community to exploit these results in future research programs.

Variable number of tandem repeats of 9 Plasmodium vivax genes among Southeast Asian isolates

The variable number of tandem repeats (VNTRs) provides valuable information about both the functional and evolutionary aspects of genetic diversity. Comparative analysis of 3 Plasmodium falciparum genomes has shown that more than 9% of its open reading frames (ORFs) harbor VNTRs. Although microsatellites and VNTR genes of P. vivax were reported, the VNTR polymorphism of genes has not been examined widely. In this study, 230 P. vivax genes were analyzed for VNTRs by SERV, and 33 kinds of TR deletions or insertions from 29 P. vivax genes (12.6%) were found. Of these, 9 VNTR fragments from 8 P. vivax genes were used for PCR amplification and sequence analysis to examine the genetic diversity among 134 isolates from four Southeast Asian countries (China, Republic of Korea, Thailand, and Myanmar) with different malaria endemicity. We confirmed the existence of extensive polymorphism of VNTR fragments in field isolates. This detection provides several suitable markers for analysis of the molecular epidemiology of P. vivax field isolates.