Every microsatellite is different: Intrinsic DNA features dictate mutagenesis of common microsatellites present in the human genome

Beyond PD(L)-1 Blockade in Microsatellite-Instable Cancers: Current Landscape of Immune Co-Inhibitory Receptor Targeting

High microsatellite instability (MSI-H) derives from genomic hypermutability due to deficient mismatch repair function. Colorectal (CRC) and endometrial cancers (EC) are the tumor types that more often present MSI-H. Anti-PD(L)-1 antibodies have been demonstrated to be agnostically effective in patients with MSI-H cancer, but 50-60% of them do not respond to single-agent treatment, highlighting the necessity of expanding their treatment opportunities. Ipilimumab (anti-CTLA4) is the only immune checkpoint inhibitor (ICI) non-targeting PD(L)-1 that has been approved so far by the FDA for MSI-H cancer, namely, CRC in combination with nivolumab. Anti-TIM3 antibody LY3321367 showed interesting clinical activity in combination with anti-PDL-1 antibody in patients with MSI-H cancer not previously treated with anti-PD(L)-1. In contrast, no clinical evidence is available for anti-LAG3, anti-TIGIT, anti-BTLA, anti-ICOS and anti-IDO1 antibodies in MSI-H cancers, but clinical trials are ongoing. Other immunotherapeutic strategies under study for MSI-H cancers include vaccines, systemic immunomodulators, STING agonists, PKM2 activators, T-cell immunotherapy, LAIR-1 immunosuppression reversal, IL5 superagonists, oncolytic viruses and IL12 partial agonists. In conclusion, several combination therapies of ICIs and novel strategies are emerging and may revolutionize the treatment paradigm of MSI-H patients in the future. A huge effort will be necessary to find reliable immune biomarkers to personalize therapeutical decisions.

Genome-wide identification and characterization of microsatellite markers in Bactrian Camel

Simple sequence repeats (SSRs) have been widely used for parentage testing, marker-assisted selection, and evolutionary studies. The insufficient availability of SSR markers in Bactrian camels partially accounts for the lack of systematic breeding. Therefore, we aimed to establish a comprehensive SSR dataset for the Bactrian camel. Our approach involved genome searching to locate every SSR in the genome, SSR-enriched sequencing to acquire polymorphism information, and literature research to collect published data. The resulting dataset contains 213,711 SSRs and details their characteristics, including genome coordinates, motifs, lengths, annotations, PCR primers, and polymorphism information. The dataset reveals a biased distribution of SSRs in the Bactrian camel genome, reflecting the mutation mechanism and complex evolution of SSRs. In practice, we successfully demonstrated the utility of the dataset through parentage testing using 15 randomly selected SSRs. This comprehensive dataset can facilitate systematic breeding and enable QTL mapping and GWAS of the Bactrian camel.

A hypervariable intron of the STAYGREEN locus provides excellent discrimination among Pisum fulvum accessions and reveals evidence for a relatively recent hybridization event with Pisum sativum

An analysis of 82 non-synonymous Pisum fulvum accessions for sequence variation in a fragment of the STAYGREEN (SGR) locus revealed 57 alleles, most of which differed in indel structure. Eight additional P. fulvum accessions, each supposedly synonymous with a different accession of the initial group, were also analyzed. In every case the paired synonymous accessions possessed the same SGR sequence but varied slightly for a 6-trait morphological phenotype, indicating that SGR sequence is a much more reliable indicator of accession identity than is a morphological characterization. SGR sequence analysis confirmed our previous finding that P. fulvum accessions separate into two allele groups. This division was not supported by results of previous studies that were based on sequences distributed across the entire genome, suggesting that the division may have been produced by selection at a nearby locus and that the SGR phylogeny may not be good indicator of overall relationships within the species. One P. fulvum accession, PI 595941 (=JI1796), displayed an SGR sequence outside the variation typical of the species. Instead, its allele resembled alleles limited to a set of Pisum sativum landraces from the Middle East, suggesting hybridization between ancestors of PI 595941 and some primitive form of domesticated P. sativum. With one exception from the extreme northwest corner of Israel, P. fulvum accessions collected north of latitude 35.5° N were fixed for alleles from group A. These northern accessions also displayed greatly reduced SGR sequence diversity compared to group A accessions collected from other regions, suggesting that the northern populations may represent recent extensions of the range of the species. Group B accessions were distributed from Lake Tiberias south and were generally sympatric with the southern group A accessions. Although group B accessions occupied a smaller area than group A, the SGR sequence diversity in this group (28 alleles in 33 accessions) exceeded that for group A.

Sequence variants affecting the genome-wide rate of germline microsatellite mutations

Microsatellites are polymorphic tracts of short tandem repeats with one to six base-pair (bp) motifs and are some of the most polymorphic variants in the genome. Using 6084 Icelandic parent-offspring trios we estimate 63.7 (95% CI: 61.9-65.4) microsatellite de novo mutations (mDNMs) per offspring per generation, excluding one bp repeats motifs (homopolymers) the estimate is 48.2 mDNMs (95% CI: 46.7-49.6). Paternal mDNMs occur at longer repeats than maternal ones, which are in turn larger with a mean size of 3.4 bp vs 3.1 bp for paternal ones. mDNMs increase by 0.97 (95% CI: 0.90-1.04) and 0.31 (95% CI: 0.25-0.37) per year of father's and mother's age at conception, respectively. Here, we find two independent coding variants that associate with the number of mDNMs transmitted to offspring; The minor allele of a missense variant (allele frequency (AF) = 1.9%) in MSH2, a mismatch repair gene, increases transmitted mDNMs from both parents (effect: 13.1 paternal and 7.8 maternal mDNMs). A synonymous variant (AF = 20.3%) in NEIL2, a DNA damage repair gene, increases paternally transmitted mDNMs (effect: 4.4 mDNMs). Thus, the microsatellite mutation rate in humans is in part under genetic control.

The sequence of the repetitive motif influences the frequency of multistep mutations in Short Tandem Repeats

Microsatellites, or Short Tandem Repeats (STRs), are subject to frequent length mutations that involve the loss or gain of an integer number of repeats. This work aimed to investigate the correlation between STRs' specific repetitive motif composition and mutational dynamics, specifically the occurrence of single- or multistep mutations. Allelic transmission data, comprising 323,818 allele transfers and 1,297 mutations, were gathered for 35 Y-chromosomal STRs with simple structure. Six structure groups were established: ATT, CTT, TCTA/GATA, GAAA/CTTT, CTTTT, and AGAGAT, according to the repetitive motif present in the DNA leading strand of the markers. Results show that the occurrence of multistep mutations varies significantly among groups of markers defined by the repetitive motif. The group of markers with the highest frequency of multistep mutations was the one with repetitive motif CTTTT (25% of the detected mutations) and the lowest frequency corresponding to the group with repetitive motifs TCTA/GATA (0.93%). Statistically significant differences (α = 0.05) were found between groups with repetitive motifs with different lengths, as is the case of TCTA/GATA and ATT (p = 0.0168), CTT (p < 0.0001) and CTTTT (p < 0.0001), as well as between GAAA/CTTT and CTTTT (p = 0.0102). The same occurred between the two tetrameric groups GAAA/CTTT and TCTA/GATA (p < 0.0001) - the first showing 5.7 times more multistep mutations than the second. When considering the number of repeats of the mutated paternal alleles, statistically significant differences were found for alleles with 10 or 12 repeats, between GATA and ATT structure groups. These results, which demonstrate the heterogeneity of mutational dynamics across repeat motifs, have implications in the fields of population genetics, epidemiology, or phylogeography, and whenever STR mutation models are used in evolutionary studies in general.

Systems biology of the genomes' microsatellite signature of Orthopoxvirus including the Monkeypox virus

This study is an attempt to extract and analyse the microsatellites or simple sequence repeats (SSRs) from the genomes of eight species of the genus Orthopoxvirus. The average size of genomes included in the study was 205 kb while the GC% was 33% for all but one. A total of 10,584 SSRs and 854 cSSRs were observed. POX2 with the largest genome of 224.499 kb had maximum of 1493 SSRs and 121 cSSRs (compound SSR) while POX7 with the smallest genome of 185.578 kb had minimum incident SSRs and cSSRs at 1181 and 96, respectively. There was significant correlation between genome size and SSR incidence. Di-nucleotide repeats were the most prevalent (57.47%) followed by mono- at 33% and tri- at 8.6%. Mono-nucleotide SSRs were predominantly T (51%) and A (48.4%). A majority of 80.32% SSRs were in the coding region. The three most similar genomes as per heat map POX1, POX7 and POX5 (93% similarity) are adjacent to one another in the phylogenetic tree. Ankyrin/Ankyrin like protein and Kelch protein which are associated with host determination and divergence have the highest SSR density in almost all studied viruses. Thus, SSRs are involved in genome evolution and host determination of viruses.

ETV6 dependency in Ewing sarcoma by antagonism of EWS-FLI1-mediated enhancer activation

The EWS-FLI1 fusion oncoprotein deregulates transcription to initiate the paediatric cancer Ewing sarcoma. Here we used a domain-focused CRISPR screen to implicate the transcriptional repressor ETV6 as a unique dependency in this tumour. Using biochemical assays and epigenomics, we show that ETV6 competes with EWS-FLI1 for binding to select DNA elements enriched for short GGAA repeat sequences. Upon inactivating ETV6, EWS-FLI1 overtakes and hyper-activates these cis-elements to promote mesenchymal differentiation, with SOX11 being a key downstream target. We show that squelching of ETV6 with a dominant-interfering peptide phenocopies these effects and suppresses Ewing sarcoma growth in vivo. These findings reveal targeting of ETV6 as a strategy for neutralizing the EWS-FLI1 oncoprotein by reprogramming of genomic occupancy.

De novo variants and recombination at 4q35: Hints for preimplantation genetic testing in facioscapulohumeral muscular dystrophy

Facioscapulohumeral muscular dystrophy (FSHD) has been associated with the deletion of an integral number of 3.3 kb units of the polymorphic D4Z4 repeat array at 4q35. The prenatal identification of this defect can be carried out on chorionic villi or amniocytes, whereas preimplantation genetic testing for monogenic disorders (PGT-M) requires molecular markers linked to the D4Z4 allele of reduced size. In this context the reliability of this association is crucial. To test the informativeness of the nearby polymorphic markers we investigated recombination at 4q35 using the polymorphic markers D4S1523, D4S163 and D4S139 positioned at 0.55, 0.5 and 0.21 Mb proximal to the D4Z4 array respectively. We determined the probability of recombination events to occur in the D4Z4-D4S1523 interval considering 86 subjects belonging to 12 FSHD families and found a recombination frequency of 14% between D4Z4 and D4S1523. Our study also revealed the occurrence of de novo variants and germline mosaicism. These findings highlight the recombinogenic nature of the 4q subtelomere and indicate that caution should be taken when interpreting PGT-M results. It is advisable that a woman who underwent a PGT-M cycle undertakes a prenatal DNA analysis to confirm the size of the D4Z4 alleles carried by the fetus.

The whole chloroplast genome of Mesua ferrea: Insight into the dynamic pattern of evolution and its comparison with species from recently diverged families

Mesua ferrea is an important source of timber, oil and herbal medicines. In the present investigation, we assembled the whole chloroplast genome of M. ferrea of size 161.4 kb. The genome contained 86 protein-coding genes, 38 tRNAs, 8 rRNA genes and exhibited a characteristic quadripartite structural orientation, with two inverted repeats (27,614 bp) separated by an LSC (88,746 bp) region and an SSC (27,614 bp) (17,470 bp). Interestingly, no gene loss was identified in the M. ferrea genome, contrary to what has been observed in other Clusioid species. We compared the chloroplast genome of M. ferrea with the chloroplast genome of Bonnetia and Garcinia belonging to Bonnetiaceae and Clusiaceae families. Overall, the compared genomes possess a similar synteny of gene order except for a small inversion in Garcinia species. M. ferrea has the largest chloroplast genome size in Clusioid clade owing to the lengthening of the LSC, IR, and non-coding regions. Substantial differences were observed in population of simple sequence repeats (SSRs) and RNA editing sites among the studied genomes. A comparative assessment of chloroplast genomes revealed five highly divergence regions: rpl32, trnS-GCU_trnG-UCC, petN-psbM, psbZ_trnG-GCC and ccsA_ndhD among the analyzed sequences. Phylogenetic analyses and sequence homology search indicate that M. ferrea is closely related to the Garcinia species.

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.

Clinical prospects of WRN inhibition as a treatment for MSI tumours

The discovery of synthetic lethal interactions with genetic deficiencies in cancers has highlighted several candidate targets for drug development, with variable clinical success. Recent work has unveiled a promising synthetic lethal interaction between inactivation/inhibition of the WRN DNA helicase and tumours with microsatellite instability, a phenotype that arises from DNA mismatch repair deficiency. While these and further studies have highlighted the therapeutic potential of WRN inhibitors, compounds with properties suitable for clinical exploitation remain to be described. Furthermore, the complexities of MSI development and its relationship to cancer evolution pose challenges for clinical prospects. Here, we discuss possible paths of MSI tumour development, the viability of WRN inhibition as a strategy in different scenarios, and the necessary conditions to create a roadmap towards successful implementation of WRN inhibitors in the clinic.

Mismatch repair is a double-edged sword in the battle against microsatellite instability

Roughly 3% of the human genome consists of microsatellites or tracts of short tandem repeats (STRs). These STRs are often unstable, undergoing high-frequency expansions (increases) or contractions (decreases) in the number of repeat units. Some microsatellite instability (MSI) is seen at multiple STRs within a single cell and is associated with certain types of cancer. A second form of MSI is characterised by expansion of a single gene-specific STR and such expansions are responsible for a group of 40+ human genetic disorders known as the repeat expansion diseases (REDs). While the mismatch repair (MMR) pathway prevents genome-wide MSI, emerging evidence suggests that some MMR factors are directly involved in generating expansions in the REDs. Thus, MMR suppresses some forms of expansion while some MMR factors promote expansion in other contexts. This review will cover what is known about the paradoxical effect of MMR on microsatellite expansion in mammalian cells.

Estimations of Mutation Rates Depend on Population Allele Frequency Distribution: The Case of Autosomal Microsatellites

Microsatellites (or short-tandem repeats (STRs)) are widely used in anthropology and evolutionary studies. Their extensive polymorphism and rapid evolution make them the ideal genetic marker for dating events, such as the age of a gene or a population. This usage requires the estimation of mutation rates, which are usually estimated by counting the observed Mendelian incompatibilities in one-generation familial configurations (typically parent(s)–child duos or trios). Underestimations are inevitable when using this approach, due to the occurrence of mutational events that do not lead to incompatibilities with the parental genotypes (‘hidden’ or ‘covert’ mutations). It is known that the likelihood that one mutation event leads to a Mendelian incompatibility depends on the mode of genetic transmission considered, the type of familial configuration (duos or trios) considered, and the genotype(s) of the progenitor(s). In this work, we show how the magnitude of the underestimation of autosomal microsatellite mutation rates varies with the populations’ allele frequency distribution spectrum. The Mendelian incompatibilities approach (MIA) was applied to simulated parent(s)/offspring duos and trios in different populational scenarios. The results showed that the magnitude and type of biases depend on the population allele frequency distribution, whatever the type of familial data considered, and are greater when duos, instead of trios, are used to obtain the estimates. The implications for molecular anthropology are discussed and a simple framework is presented to correct the naïf estimates, along with an informatics tool for the correction of incompatibility rates obtained through the MIA.

CLCC1 c. 75C>A Mutation in Pakistani Derived Retinitis Pigmentosa Families Likely Originated With a Single Founder Mutation 2,000-5,000 Years Ago

Background: A CLCC1 c. 75C > A (p.D25E) mutation has been associated with autosomal recessive pigmentosa in patients in and from Pakistan. CLCC1 is ubiquitously expressed, and knockout models of this gene in zebrafish and mice are lethal in the embryonic period, suggesting that possible retinitis pigmentosa mutations in this gene might be limited to those leaving partial activity. In agreement with this hypothesis, the mutation is the only CLCC1 mutation associated with retinitis pigmentosa to date, and all identified patients with this mutation share a common SNP haplotype surrounding the mutation, suggesting a common founder.

Methods: SNPs were genotyped by a combination of WGS and Sanger sequencing. The original founder haplotype, and recombination pathways were delineated by examination to minimize recombination events. Mutation age was estimated by four methods including an explicit solution, an iterative approach, a Bayesian approach and an approach based solely on ancestral segment lengths using high density SNP data.

Results: All members of each of the nine families studied shared a single autozygous SNP haplotype for the CLCC1 region ranging from approximately 1–3.5 Mb in size. The haplotypes shared by the families could be derived from a single putative ancestral haplotype with at most two recombination events. Based on the haplotype and Gamma analysis, the estimated age of the founding mutation varied from 79 to 196 generations, or approximately 2,000–5,000 years, depending on the markers used in the estimate. The DMLE (Bayesian) estimates ranged from 2,160 generations assuming a population growth rate of 0–309 generations assuming a population growth rate of 2% with broad 95% confidence intervals.

Conclusion: These results provide insight into the origin of the CLCC1 mutation in the Pakistan population. This mutation is estimated to have occurred 2000–5,000 years ago and has been transmitted to affected families of Pakistani origin in geographically dispersed locations around the world. This is the only mutation in CLCC1 identified to date, suggesting that the CLCC1 gene is under a high degree of constraint, probably imposed by functional requirements for this gene during embryonic development.

Improved Models of Coalescence Ages of Y-DNA Haplogroups

Databases of commercial DNA-testing companies now contain more customers with sequenced DNA than any completed academic study, leading to growing interest from academic and forensic entities. An important result for both these entities and the test takers themselves is how closely two individuals are related in time, as calculated through one or more molecular clocks. For Y-DNA, existing interpretations of these clocks are insufficiently accurate to usefully measure relatedness in historic times. In this article, I update the methods used to calculate coalescence ages (times to most-recent common ancestor, or TMRCAs) using a new, probabilistic statistical model that includes Y-SNP, Y-STR and ancilliary historical data, and provide examples of its use.

Next generation sequencing of Y-STRs in father-son pairs and comparison with traditional capillary electrophoresis

To evaluate the promising advantages of massively parallel sequencing (MPS) in our casework, we analysed a total of 33 Y-chromosomal short tandem repeats (Y-STRs) with traditional capillary electrophoresis (CE) and 25 Y-STRs using the newer MPS technology. We studied the outcome of both technologies in 64 father-son pairs using stock and custom-designed kits. Current MPS technology confirmed the 13 mutational events observed with CE and improved our understanding of the complex nature of STR mutations. By detecting isometric sequence variants between unrelated males, we show that sequencing Y-STRs using MPS can boost discrimination power.

Gynecological Cancers Caused by Deficient Mismatch Repair and Microsatellite Instability

Mutations in mismatch repair genes leading to mismatch repair (MMR) deficiency (dMMR) and microsatellite instability (MSI) have been implicated in multiple types of gynecologic malignancies. Endometrial carcinoma represents the largest group, with approximately 30% of these cancers caused by dMMR/MSI. Thus, testing for dMMR is now routine for endometrial cancer. Somatic mutations leading to dMMR account for approximately 90% of these cancers. However, in 5-10% of cases, MMR protein deficiency is due to a germline mutation in the mismatch repair genes MLH1, MSH2, MSH6, PMS2, or EPCAM. These germline mutations, known as Lynch syndrome, are associated with an increased risk of both endometrial and ovarian cancer, in addition to colorectal, gastric, urinary tract, and brain malignancies. So far, gynecological cancers with dMMR/MSI are not well characterized and markers for detection of MSI in gynecological cancers are not well defined. In addition, currently advanced endometrial cancers have a poor prognosis and are treated without regard to MSI status. Elucidation of the mechanism causing dMMR/MSI gynecological cancers would aid in diagnosis and therapeutic intervention. Recently, a new immunotherapy was approved for the treatment of solid tumors with MSI that have recurred or progressed after failing traditional treatment strategies. In this review, we summarize the MMR defects and MSI observed in gynecological cancers, their prognostic value, and advances in therapeutic strategies to treat these cancers.

The Biological Significance of Multi-copy Regions and Their Impact on Variant Discovery

Identification of genetic variants via high-throughput sequencing (HTS) technologies has been essential for both fundamental and clinical studies. However, to what extent the genome sequence composition affects variant calling remains unclear. In this study, we identified 63,897 multi-copy sequences (MCSs) with a minimum length of 300 bp, each of which occurs at least twice in the human genome. The 151,749 genomic loci (multi-copy regions, or MCRs) harboring these MCSs account for 1.98% of the genome and are distributed unevenly across chromosomes. MCRs containing the same MCS tend to be located on the same chromosome. Gene Ontology (GO) analyses revealed that 3800 genes whose UTRs or exons overlap with MCRs are enriched for Golgi-related cellular component terms and various enzymatic activities in the GO biological function category. MCRs are also enriched for loci that are sensitive to neocarzinostatin-induced double-strand breaks. Moreover, genetic variants discovered by genome-wide association studies and recorded in dbSNP are significantly underrepresented in MCRs. Using simulated HTS datasets, we show that false variant discovery rates are significantly higher in MCRs than in other genomic regions. These results suggest that extra caution must be taken when identifying genetic variants in the MCRs via HTS technologies.

Exhaled Breath Condensate: A Non-Invasive Source for Tracking of Genetic and Epigenetic Alterations in Lung Diseases

Lung diseases have been recognized as an extensive cause of morbidity and mortality in the worldwide. The high degree of clinical heterogeneity and nonspecific initial symptoms of lung diseases contribute to a delayed diagnosis. So, the molecular and genomic profiling play a pivotal role in promoting the pulmonary diseases. Exhaled breath condensate (EBC) as a novel and potential method for sampling the respiratory epithelial lining fluid is to assess the inflammatory and oxidative stress biomarkers, drugs and genetic alterations in the pathophysiologic processes of lung diseases. The recent studies on the analysis of EBC from both a genetic and epigenetic point of view were searched from database and reviewed. This review provides an overview of the current findings in the tracking of genomic and epigenetic alterations which are potentially effective in better management of cancer detection. In addition, respiratory microbiota DNA using EBC samples in association with pulmonary disease especially lung cancer were investigated. Various studies have concluded that EBC has a great potential for analysis of nuclear and mitochondrial DNA alterations as well as epigenetic modifications and identification of respiratory microbiome. Next-generation sequencing (NGS) based genomic profiling of EBC samples is recommended as a promising approach to establish personalized based prevention, diagnosis, treatment and post-treatment follow-ups for patients with lung diseases especially lung cancer.

Intraspecific Variation in Microsatellite Mutation Profiles in Daphnia magna

Microsatellite loci (tandem repeats of short nucleotide motifs) are highly abundant in eukaryotic genomes and often used as genetic markers because they can exhibit variation both within and between populations. Although widely recognized for their mutability and utility, the mutation rates of microsatellites have only been empirically estimated in a few species, and have rarely been compared across genotypes and populations within a species. Here, we investigate the dynamics of microsatellite mutation over long- and short-time periods by quantifying the starting abundance and mutation rates for microsatellites for six different genotypes of Daphnia magna, an aquatic microcrustacean, collected from three populations (Finland, Germany, and Israel). Using whole-genome sequences of these six starting genotypes, descendent mutation accumulation (MA) lines, and large population controls (non-MA lines), we find each genotype exhibits a distinctive initial microsatellite profile which clusters according to the population-of-origin. During the period of MA, we observe motif-specific, highly variable, and rapid microsatellite mutation rates across genotypes of D. magna, the average of which is order of magnitude greater than the recently reported rate observed in a single genotype of the congener, Daphnia pulex. In our experiment, genotypes with more microsatellites starting out exhibit greater losses and those with fewer microsatellites starting out exhibit greater gains—a context-dependent mutation bias that has not been reported previously. We discuss how genotype-specific mutation rates and spectra, in conjunction with evolutionary forces, can shape both the differential accumulation of repeat content in the genome and the evolution of mutation rates.

Defects in the GINS complex increase the instability of repetitive sequences via a recombination-dependent mechanism

Faithful replication and repair of DNA lesions ensure genome maintenance. During replication in eukaryotic cells, DNA is unwound by the CMG helicase complex, which is composed of three major components: the Cdc45 protein, Mcm2-7, and the GINS complex. The CMG in complex with DNA polymerase epsilon (CMG-E) participates in the establishment and progression of the replisome. Impaired functioning of the CMG-E was shown to induce genomic instability and promote the development of various diseases. Therefore, CMG-E components play important roles as caretakers of the genome. In Saccharomyces cerevisiae, the GINS complex is composed of the Psf1, Psf2, Psf3, and Sld5 essential subunits. The Psf1-1 mutant form fails to interact with Psf3, resulting in impaired replisome assembly and chromosome replication. Here, we show increased instability of repeat tracts (mononucleotide, dinucleotide, trinucleotide and longer) in yeast psf1-1 mutants. To identify the mechanisms underlying this effect, we analyzed repeated sequence instability using derivatives of psf1-1 strains lacking genes involved in translesion synthesis, recombination, or mismatch repair. Among these derivatives, deletion of RAD52, RAD51, MMS2, POL32, or PIF1 significantly decreased DNA repeat instability. These results, together with the observed increased amounts of single-stranded DNA regions and Rfa1 foci suggest that recombinational mechanisms make important contributions to repeat tract instability in psf1-1 cells. We propose that defective functioning of the CMG-E complex in psf1-1 cells impairs the progression of DNA replication what increases the contribution of repair mechanisms such as template switch and break-induced replication. These processes require sequence homology search which in case of a repeated DNA tract may result in misalignment leading to its expansion or contraction.

Processes that ensure genome stability are crucial for all organisms to avoid mutations and decrease the risk of diseases. The coordinated activity of mechanisms underlying the maintenance of high-fidelity DNA duplication and repair is critical to deal with the malfunction of replication forks or DNA damage. Repeated sequences in DNA are particularly prone to instability; these sequences undergo expansions or contractions, leading in humans to various neurological, neurodegenerative, and neuromuscular disorders. A mutant form of one of the noncatalytic subunits of active DNA helicase complex impairs DNA replication. Here, we show that this form also significantly increases the instability of mononucleotide, dinucleotide, trinucleotide and longer repeat tracts. Our results suggest that in cells that harbor a mutated variant of the helicase complex, continuation of DNA replication is facilitated by recombination processes, and this mechanism can be highly mutagenic during repair synthesis through repetitive regions, especially regions that form secondary structures. Our results indicate that proper functioning of the DNA helicase complex is crucial for maintenance of the stability of repeated DNA sequences, especially in the context of recently described disorders in which mutations or deregulation of the human homologs of genes encoding DNA helicase subunits were observed.

Population-Based Cytogenetic Banding Analysis and Phylogenetic Relationships of the Neotropical Fungus-Farming Ant Trachymyrmex holmgreni Wheeler, 1925

Trachymyrmex is one of the most species-rich genera within fungus-farming ants and presents intraspecific cytogenetic polymorphisms as well as possible cryptic species. This ant genus is currently paraphyletic. Therefore, to unravel systematic and taxonomic misunderstandings, it is necessary to incorporate new information. We aimed to cytogenetically and genetically examine Trachymyrmex holmgreni populations from southern and northern Brazil to identify intraspecific chromosomal variations that support incipient speciation and reveal the species' position in a molecular phylogeny. Our cytogenetic approach did not show population variation in the mapping of both 18S rDNA and the TTAGG(6) motif, presenting instead a pattern characteristic of correlated species. However, the clustered pattern of the microsatellite GA(15) showed significant differences among populations: a well-defined block in each homologue, distinctly irregular signs between homologues, and blocks in 2 pairs of homologues. Our phylogenetic reconstruction yielded unexpected results, grouping representatives of 3 former morphological groups into 1 clade, namely T. urichii, T. papulatus, and T. holmgreni. Previously, it was suggested that northern and southern populations of T. holmgreni may be undergoing incipient speciation, but we can only indicate that the southernmost population differs prominently from the others in its distribution pattern of the microsatellite GA(15). Our study also supports the uniformity of karyotypes and repetitive DNA from both telomeric sequences and ribosomal DNA in Trachymyrmex studied here. In addition, we clarify some phylogenetic uncertainties within the genus and suggest further relevant systematic changes. Finally, additional studies utilizing other probes and additional populations may allow the detection of hidden genetic variation.

Loss of genetic diversity, recovery and allele surfing in a colonizing parasite, Geomydoecus aurei

Understanding the genetic consequences of changes in species distributions has wide-ranging implications for predicting future outcomes of climate change, for protecting threatened or endangered populations and for understanding the history that has led to current genetic patterns within species. Herein, we examine the genetic consequences of range expansion over a 25-year period in a parasite (Geomydoecus aurei) that is in the process of expanding its geographic range via invasion of a novel host. By sampling the genetics of 1,935 G. aurei lice taken from 64 host individuals collected over this time period using 12 microsatellite markers, we test hypotheses concerning linear spatial expansion, genetic recovery time and allele surfing. We find evidence of decreasing allelic richness (AR) with increasing distance from the source population, supporting a linear, stepping stone model of spatial expansion that emphasizes the effects of repeated bottleneck events during colonization. We provide evidence of post-bottleneck genetic recovery, with average AR of infrapopulations increasing about 30% over the 225-generation span of time observed directly in this study. Our estimates of recovery rate suggest, however, that recovery has plateaued and that this population may not reach genetic diversity levels of the source population without further immigration from the source population. Finally, we employ a grid-based sampling scheme in the region of ongoing population expansion and provide empirical evidence for the power of allele surfing to impart genetic structure on a population, even under conditions of selective neutrality and in a place that lacks strong barriers to gene flow.

Human copy number variants are enriched in regions of low mappability

Copy number variants (CNVs) are known to affect a large portion of the human genome and have been implicated in many diseases. Although whole-genome sequencing (WGS) can help identify CNVs, most analytical methods suffer from limited sensitivity and specificity, especially in regions of low mappability. To address this, we use PopSV, a CNV caller that relies on multiple samples to control for technical variation. We demonstrate that our calls are stable across different types of repeat-rich regions and validate the accuracy of our predictions using orthogonal approaches. Applying PopSV to 640 human genomes, we find that low-mappability regions are approximately 5 times more likely to harbor germline CNVs, in stark contrast to the nearly uniform distribution observed for somatic CNVs in 95 cancer genomes. In addition to known enrichments in segmental duplication and near centromeres and telomeres, we also report that CNVs are enriched in specific types of satellite and in some of the most recent families of transposable elements. Finally, using this comprehensive approach, we identify 3455 regions with recurrent CNVs that were missing from existing catalogs. In particular, we identify 347 genes with a novel exonic CNV in low-mappability regions, including 29 genes previously associated with disease.

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.

Integrated analysis sheds light on evolutionary trajectories of young transcription start sites in the human genome

Understanding the molecular mechanisms and evolution of the gene regulatory system remains a major challenge in biology. Transcription start sites (TSSs) are especially interesting because they are central to initiating gene expression. Previous studies revealed widespread transcription initiation and fast turnover of TSSs in mammalian genomes. Yet, how new TSSs originate and how they evolve over time remain poorly understood. To address these questions, we analyzed ∼200,000 human TSSs by integrating evolutionary (inter- and intra-species) and functional genomic data, particularly focusing on evolutionarily young TSSs that emerged in the primate lineage. TSSs were grouped according to their evolutionary age using sequence alignment information as a proxy. Comparisons of young and old TSSs revealed that (1) new TSSs emerge through a combination of intrinsic factors, like the sequence properties of transposable elements and tandem repeats, and extrinsic factors such as their proximity to existing regulatory modules; (2) new TSSs undergo rapid evolution that reduces the inherent instability of repeat sequences associated with a high propensity of TSS emergence; and (3) once established, the transcriptional competence of surviving TSSs is gradually enhanced, with evolutionary changes subject to temporal (fewer regulatory changes in younger TSSs) and spatial constraints (fewer regulatory changes in more isolated TSSs). These findings advance our understanding of how regulatory innovations arise in the genome throughout evolution and highlight the genomic robustness and evolvability in these processes.

Comprehensive Characterization of Simple Sequence Repeats in Eggplant (Solanum melongena L.) Genome and Construction of a Web Resource

We have characterized the simple sequence repeat (SSR) markers of the eggplant (Solanum melongena) using a recent high quality sequence of its whole genome. We found nearly 133,000 perfect SSRs, a density of 125.5 SSRs/Mbp, and also about 178,400 imperfect SSRs. Of the perfect SSRs, 15.6% were complex, with two stretches of repeats separated by an intervening block of <100 nt. Di- and trinucleotide SSRs accounted, respectively, for 43 and 37% of the total. The SSRs were classified according to their number of repeats and overall length, and were assigned to their linkage group. We found 2,449 of the perfect SSRs in 2,086 genes, with an overall density of 18.5 SSRs/Mbp across the gene space; 3,524 imperfect SSRs were present in 2,924 genes at a density of 26.7 SSRs/Mbp. Putative functions were assigned via ontology to genes containing at least one SSR. Using this data we developed an "Eggplant Microsatellite DataBase" (EgMiDB) which permits identification of SSR markers in terms of their location on the genome, type of repeat (perfect vs. imperfect), motif type, sequence, repeat number and genomic/gene context. It also suggests forward and reverse primers. We employed an in silico PCR analysis to validate these SSR markers, using as templates two CDS sets and three assembled transcriptomes obtained from diverse eggplant accessions.

Role of specialized DNA polymerases in the limitation of replicative stress and DNA damage transmission

Replication stress is a strong and early driving force for genomic instability and tumor development. Beside replicative DNA polymerases, an emerging group of specialized DNA polymerases is involved in the technical assistance of the replication machinery in order to prevent replicative stress and its deleterious consequences. During S-phase, altered progression of the replication fork by endogenous or exogenous impediments induces replicative stress, causing cells to reach mitosis with genomic regions not fully duplicated. Recently, specific mechanisms to resolve replication intermediates during mitosis with the aim of limiting DNA damage transmission to daughter cells have been identified. In this review, we detail the two major actions of specialized DNA polymerases that limit DNA damage transmission: the prevention of replicative stress by non-B DNA replication and the recovery of stalled replication forks.

Regulation and Modulation of Human DNA Polymerase δ Activity and Function

This review focuses on the regulation and modulation of human DNA polymerase δ (Pol δ). The emphasis is on the mechanisms that regulate the activity and properties of Pol δ in DNA repair and replication. The areas covered are the degradation of the p12 subunit of Pol δ, which converts it from a heterotetramer (Pol δ4) to a heterotrimer (Pol δ3), in response to DNA damage and also during the cell cycle. The biochemical mechanisms that lead to degradation of p12 are reviewed, as well as the properties of Pol δ4 and Pol δ3 that provide insights into their functions in DNA replication and repair. The second focus of the review involves the functions of two Pol δ binding proteins, polymerase delta interaction protein 46 (PDIP46) and polymerase delta interaction protein 38 (PDIP38), both of which are multi-functional proteins. PDIP46 is a novel activator of Pol δ4, and the impact of this function is discussed in relation to its potential roles in DNA replication. Several new models for the roles of Pol δ3 and Pol δ4 in leading and lagging strand DNA synthesis that integrate a role for PDIP46 are presented. PDIP38 has multiple cellular localizations including the mitochondria, the spliceosomes and the nucleus. It has been implicated in a number of cellular functions, including the regulation of specialized DNA polymerases, mitosis, the DNA damage response, mouse double minute 2 homolog (Mdm2) alternative splicing and the regulation of the NADPH oxidase 4 (Nox4).

CRISPR-Cas9-targeted fragmentation and selective sequencing enable massively parallel microsatellite analysis

Microsatellites are multi-allelic and composed of short tandem repeats (STRs) with individual motifs composed of mononucleotides, dinucleotides or higher including hexamers. Next-generation sequencing approaches and other STR assays rely on a limited number of PCR amplicons, typically in the tens. Here, we demonstrate STR-Seq, a next-generation sequencing technology that analyses over 2,000 STRs in parallel, and provides the accurate genotyping of microsatellites. STR-Seq employs in vitro CRISPR-Cas9-targeted fragmentation to produce specific DNA molecules covering the complete microsatellite sequence. Amplification-free library preparation provides single molecule sequences without unique molecular barcodes. STR-selective primers enable massively parallel, targeted sequencing of large STR sets. Overall, STR-Seq has higher throughput, improved accuracy and provides a greater number of informative haplotypes compared with other microsatellite analysis approaches. With these new features, STR-Seq can identify a 0.1% minor genome fraction in a DNA mixture composed of different, unrelated samples.

Classification and characterization of microsatellite instability across 18 cancer types

Microsatellite instability (MSI), the spontaneous loss or gain of nucleotides from repetitive DNA tracts, is a diagnostic phenotype for gastrointestinal, endometrial, and colorectal tumors, yet the landscape of instability events across a wider variety of cancer types remains poorly understood. To explore MSI across malignancies, we examined 5,930 cancer exomes from 18 cancer types at more than 200,000 microsatellite loci and constructed a genomic classifier for MSI. We identified MSI-positive tumors in 14 of the 18 cancer types. We also identified loci that were more likely to be unstable in particular cancer types, resulting in specific instability signatures that involved cancer-associated genes, suggesting that instability patterns reflect selective pressures and can potentially identify novel cancer drivers. We also observed a correlation between survival outcomes and the overall burden of unstable microsatellites, suggesting that MSI may be a continuous, rather than discrete, phenotype that is informative across cancer types. These analyses offer insight into conserved and cancer-specific properties of MSI and reveal opportunities for improved methods of clinical MSI diagnosis and cancer gene discovery.

MIPSTR: a method for multiplex genotyping of germline and somatic STR variation across many individuals

Short tandem repeats (STRs) are highly mutable genetic elements that often reside in regulatory and coding DNA. The cumulative evidence of genetic studies on individual STRs suggests that STR variation profoundly affects phenotype and contributes to trait heritability. Despite recent advances in sequencing technology, STR variation has remained largely inaccessible across many individuals compared to single nucleotide variation or copy number variation. STR genotyping with short-read sequence data is confounded by (1) the difficulty of uniquely mapping short, low-complexity reads; and (2) the high rate of STR amplification stutter. Here, we present MIPSTR, a robust, scalable, and affordable method that addresses these challenges. MIPSTR uses targeted capture of STR loci by single-molecule Molecular Inversion Probes (smMIPs) and a unique mapping strategy. Targeted capture and our mapping strategy resolve the first challenge; the use of single molecule information resolves the second challenge. Unlike previous methods, MIPSTR is capable of distinguishing technical error due to amplification stutter from somatic STR mutations. In proof-of-principle experiments, we use MIPSTR to determine germline STR genotypes for 102 STR loci with high accuracy across diverse populations of the plant A. thaliana. We show that putatively functional STRs may be identified by deviation from predicted STR variation and by association with quantitative phenotypes. Using DNA mixing experiments and a mutant deficient in DNA repair, we demonstrate that MIPSTR can detect low-frequency somatic STR variants. MIPSTR is applicable to any organism with a high-quality reference genome and is scalable to genotyping many thousands of STR loci in thousands of individuals.

The overdue promise of short tandem repeat variation for heritability

Short tandem repeat (STR) variation has been proposed as a major explanatory factor in the heritability of complex traits in humans and model organisms. However, we still struggle to incorporate STR variation into genotype-phenotype maps. We review here the promise of STRs in contributing to complex trait heritability and highlight the challenges that STRs pose due to their repetitive nature. We argue that STR variants are more likely than single-nucleotide variants to have epistatic interactions, reiterate the need for targeted assays to genotype STRs accurately, and call for more appropriate statistical methods in detecting STR-phenotype associations. Lastly, we suggest that somatic STR variation within individuals may serve as a read-out of disease susceptibility, and is thus potentially a valuable covariate for future association studies.

The conserved PFT1 tandem repeat is crucial for proper flowering in Arabidopsis thaliana

It is widely appreciated that short tandem repeat (STR) variation underlies substantial phenotypic variation in organisms. Some propose that the high mutation rates of STRs in functional genomic regions facilitate evolutionary adaptation. Despite their high mutation rate, some STRs show little to no variation in populations. One such STR occurs in the Arabidopsis thaliana gene PFT1 (MED25), where it encodes an interrupted polyglutamine tract. Although the PFT1 STR is large (∼270 bp), and thus expected to be extremely variable, it shows only minuscule variation across A. thaliana strains. We hypothesized that the PFT1 STR is under selective constraint, due to previously undescribed roles in PFT1 function. We investigated this hypothesis using plants expressing transgenic PFT1 constructs with either an endogenous STR or synthetic STRs of varying length. Transgenic plants carrying the endogenous PFT1 STR generally performed best in complementing a pft1 null mutant across adult PFT1-dependent traits. In stark contrast, transgenic plants carrying a PFT1 transgene lacking the STR phenocopied a pft1 loss-of-function mutant for flowering time phenotypes and were generally hypomorphic for other traits, establishing the functional importance of this domain. Transgenic plants carrying various synthetic constructs occupied the phenotypic space between wild-type and pft1 loss-of-function mutants. By varying PFT1 STR length, we discovered that PFT1 can act as either an activator or repressor of flowering in a photoperiod-dependent manner. We conclude that the PFT1 STR is constrained to its approximate wild-type length by its various functional requirements. Our study implies that there is strong selection on STRs not only to generate allelic diversity, but also to maintain certain lengths pursuant to optimal molecular function.

Clinical and biochemical function of polymorphic NR0B1 GGAA-microsatellites in Ewing sarcoma: a report from the Children's Oncology Group

BACKGROUND

The genetics involved in Ewing sarcoma susceptibility and prognosis are poorly understood. EWS/FLI and related EWS/ETS chimeras upregulate numerous gene targets via promoter-based GGAA-microsatellite response elements. These microsatellites are highly polymorphic in humans, and preliminary evidence suggests EWS/FLI-mediated gene expression is highly dependent on the number of GGAA motifs within the microsatellite.

OBJECTIVES

Here we sought to examine the polymorphic spectrum of a GGAA-microsatellite within the NR0B1 promoter (a critical EWS/FLI target) in primary Ewing sarcoma tumors, and characterize how this polymorphism influences gene expression and clinical outcomes.

RESULTS

A complex, bimodal pattern of EWS/FLI-mediated gene expression was observed across a wide range of GGAA motifs, with maximal expression observed in constructs containing 20-26 GGAA motifs. Relative to white European and African controls, the NR0B1 GGAA-microsatellite in tumor cells demonstrated a strong bias for haplotypes containing 21-25 GGAA motifs suggesting a relationship between microsatellite function and disease susceptibility. This selection bias was not a product of microsatellite instability in tumor samples, nor was there a correlation between NR0B1 GGAA-microsatellite polymorphisms and survival outcomes.

CONCLUSIONS

These data suggest that GGAA-microsatellite polymorphisms observed in human populations modulate EWS/FLI-mediated gene expression and may influence disease susceptibility in Ewing sarcoma.

Microsatellite interruptions stabilize primate genomes and exist as population-specific single nucleotide polymorphisms within individual human genomes

Interruptions of microsatellite sequences impact genome evolution and can alter disease manifestation. However, human polymorphism levels at interrupted microsatellites (iMSs) are not known at a genome-wide scale, and the pathways for gaining interruptions are poorly understood. Using the 1000 Genomes Phase-1 variant call set, we interrogated mono-, di-, tri-, and tetranucleotide repeats up to 10 units in length. We detected ∼26,000–40,000 iMSs within each of four human population groups (African, European, East Asian, and American). We identified population-specific iMSs within exonic regions, and discovered that known disease-associated iMSs contain alleles present at differing frequencies among the populations. By analyzing longer microsatellites in primate genomes, we demonstrate that single interruptions result in a genome-wide average two- to six-fold reduction in microsatellite mutability, as compared with perfect microsatellites. Centrally located interruptions lowered mutability dramatically, by two to three orders of magnitude. Using a biochemical approach, we tested directly whether the mutability of a specific iMS is lower because of decreased DNA polymerase strand slippage errors. Modeling the adenomatous polyposis coli tumor suppressor gene sequence, we observed that a single base substitution interruption reduced strand slippage error rates five- to 50-fold, relative to a perfect repeat, during synthesis by DNA polymerases α, β, or η. Computationally, we demonstrate that iMSs arise primarily by base substitution mutations within individual human genomes. Our biochemical survey of human DNA polymerase α, β, δ, κ, and η error rates within certain microsatellites suggests that interruptions are created most frequently by low fidelity polymerases. Our combined computational and biochemical results demonstrate that iMSs are abundant in human genomes and are sources of population-specific genetic variation that may affect genome stability. The genome-wide identification of iMSs in human populations presented here has important implications for current models describing the impact of microsatellite polymorphisms on gene expression.

Microsatellites are short tandem repeat DNA sequences located throughout the human genome that display a high degree of inter-individual variation. This characteristic makes microsatellites an attractive tool for population genetics and forensics research. Some microsatellites affect gene expression, and mutations within such microsatellites can cause disease. Interruption mutations disrupt the perfect repeated array and are frequently associated with altered disease risk, but they have not been thoroughly studied in human genomes. We identified interrupted mono-, di-, tri- and tetranucleotide MSs (iMS) within individual genomes from African, European, Asian and American population groups. We show that many iMSs, including some within disease-associated genes, are unique to a single population group. By measuring the conservation of microsatellites between human and chimpanzee genomes, we demonstrate that interruptions decrease the probability of microsatellite mutations throughout the genome. We demonstrate that iMSs arise in the human genome by single base changes within the DNA, and provide biochemical data suggesting that these stabilizing changes may be created by error-prone DNA polymerases. Our genome-wide study supports the model in which iMSs act to stabilize individual genomes, and suggests that population-specific differences in microsatellite architecture may be an avenue by which genetic ancestry impacts individual disease risk.

Differential recombination dynamics within the MHC of macaque species

A panel of 15 carefully selected microsatellites (short tandem repeats, STRs) has allowed us to study segregation and haplotype stability in various macaque species. The STRs span the major histocompatibility complex (MHC) region and map in more detail from the centromeric part of the Mhc-A to the DR region. Two large panels of Indian rhesus and Indonesian/Indochinese cynomolgus macaques have been subjected to pedigree analysis, allowing the definition of 161 and 36 different haplotypes and the physical mapping of 10 and 5 recombination sites, respectively. Although most recombination sites within the studied section of the Indian rhesus monkey MHC are situated between the Mhc-A and Mhc-B regions, the resulting recombination rate for this genomic segment is low and similar to that in humans. In contrast, in Indonesian/Indochinese macaques, two recombination sites, which appear to be absent in rhesus macaques, map between the class III and II regions. As a result, the mean recombination frequency of the core MHC, Mhc-A to class II, is higher in Indonesian/Indochinese cynomolgus than in Indian rhesus macaques, but as such is comparable to that in humans. The present communication demonstrates that the dynamics of recombination 'hot/cold spots' in the MHC, as well as their frequencies, may differ substantially between highly related macaque species.

Estimation for general birth-death processes

Birth-death processes (BDPs) are continuous-time Markov chains that track the number of "particles" in a system over time. While widely used in population biology, genetics and ecology, statistical inference of the instantaneous particle birth and death rates remains largely limited to restrictive linear BDPs in which per-particle birth and death rates are constant. Researchers often observe the number of particles at discrete times, necessitating data augmentation procedures such as expectation-maximization (EM) to find maximum likelihood estimates. For BDPs on finite state-spaces, there are powerful matrix methods for computing the conditional expectations needed for the E-step of the EM algorithm. For BDPs on infinite state-spaces, closed-form solutions for the E-step are available for some linear models, but most previous work has resorted to time-consuming simulation. Remarkably, we show that the E-step conditional expectations can be expressed as convolutions of computable transition probabilities for any general BDP with arbitrary rates. This important observation, along with a convenient continued fraction representation of the Laplace transforms of the transition probabilities, allows for novel and efficient computation of the conditional expectations for all BDPs, eliminating the need for truncation of the state-space or costly simulation. We use this insight to derive EM algorithms that yield maximum likelihood estimation for general BDPs characterized by various rate models, including generalized linear models. We show that our Laplace convolution technique outperforms competing methods when they are available and demonstrate a technique to accelerate EM algorithm convergence. We validate our approach using synthetic data and then apply our methods to cancer cell growth and estimation of mutation parameters in microsatellite evolution.

Allelic imbalance in 1p, 7q, 9p, 11p, 12q and 16q regions in non-small cell lung carcinoma and its clinical association: a pilot study

In lung cancer pathogenesis, genetic instability, i.e., loss of heterozygosity (LOH) and microsatellite instability (MSI) is a frequent molecular event, occurring at an early stage of cancerogenesis. The presence of LOH/MSI in non-small cell lung carcinoma (NSCLC) was found in many chromosomal regions, but exclusive of 3p their diagnostic value remains controversial. In this study we focused on other than 3p regions-1p31.2, 7q32.2, 9p21.3, 11p15.5, 12q23.2 and 16q22-the loci of many oncogenes and tumour suppressor genes. To analyze the potential role of LOH/MSI involved in NSCLC pathogenesis we allelotyped a panel of 13 microsatellite markers in a group of 56 cancer specimens. Our data demonstrate the presence of allelic loss for all (13) analyzed markers. Total LOH/MSI frequency in NSCLC was the highest for chromosomal region 11p15.5 (25.84 %), followed by 9p21.3 and 1p31.2 (19.87 and 16.67 % respectively). A statistically significant increase of total LOH/MSI frequency was detected for the 11p15.5 region (p = 0.0301; χ(2) test). The associations of total LOH/MSI frequency: 1) increase in 11p15.5 region (p = 0.047; χ(2) test) and 2) decrease in 7q32.2 region (p = 0.037; χ(2) test) have been statistically significant in AJCC III (American Joint Committee on Cancer Staging). In Fractional Allele Loss (FAL) index analysis, the correlation with cigarette addiction has been statistically significant. The increased amount of cigarettes smoked (pack years) in a lifetime correlates with increasing FAL (p = 0.024; Kruskal-Wallis test). These results demonstrate that LOH/MSI alternation in studied chromosomal regions is strongly influenced by tobacco smoking but do not seem to be pivotal NSCLC diagnostic marker with prognostic impact.

A unique epigenetic signature is associated with active DNA replication loci in human embryonic stem cells

The cellular epigenetic landscape changes as pluripotent stem cells differentiate to somatic cells or when differentiated cells transform to a cancerous state. These epigenetic changes are commonly correlated with differences in gene expression. Whether active DNA replication is also associated with distinct chromatin environments in these developmentally and phenotypically diverse cell types has not been known. Here, we used BrdU-seq to map active DNA replication loci in human embryonic stem cells (hESCs), normal primary fibroblasts and a cancer cell line, and correlated these maps to the epigenome. In all cell lines, the majority of BrdU peaks were enriched in euchromatin and at DNA repetitive elements, especially at microsatellite repeats, and coincided with previously determined replication origins. The most prominent BrdU peaks were shared between all cells but a sizable fraction of the peaks were specific to each cell type and associated with cell type-specific genes. Surprisingly, the BrdU peaks that were common to all cell lines were associated with H3K18ac, H3K56ac, and H4K20me1 histone marks only in hESCs but not in normal fibroblasts or cancer cells. Depletion of the histone acetyltransferases for H3K18 and H3K56 dramatically decreased the number and intensity of BrdU peaks in hESCs. Our data reveal a unique epigenetic signature that distinguishes active replication loci in hESCs from normal somatic or malignant cells.

Mutation rates, spectra, and genome-wide distribution of spontaneous mutations in mismatch repair deficient yeast

DNA mismatch repair is a highly conserved DNA repair pathway. In humans, germline mutations in hMSH2 or hMLH1, key components of mismatch repair, have been associated with Lynch syndrome, a leading cause of inherited cancer mortality. Current estimates of the mutation rate and the mutational spectra in mismatch repair defective cells are primarily limited to a small number of individual reporter loci. Here we use the yeast Saccharomyces cerevisiae to generate a genome-wide view of the rates, spectra, and distribution of mutation in the absence of mismatch repair. We performed mutation accumulation assays and next generation sequencing on 19 strains, including 16 msh2 missense variants implicated in Lynch cancer syndrome. The mutation rate for DNA mismatch repair null strains was approximately 1 mutation per genome per generation, 225-fold greater than the wild-type rate. The mutations were distributed randomly throughout the genome, independent of replication timing. The mutation spectra included insertions/deletions at homopolymeric runs (87.7%) and at larger microsatellites (5.9%), as well as transitions (4.5%) and transversions (1.9%). Additionally, repeat regions with proximal repeats are more likely to be mutated. A bias toward deletions at homopolymers and insertions at (AT)_n microsatellites suggests a different mechanism for mismatch generation at these sites. Interestingly, 5% of the single base pair substitutions might represent double-slippage events that occurred at the junction of immediately adjacent repeats, resulting in a shift in the repeat boundary. These data suggest a closer scrutiny of tumor suppressors with homopolymeric runs with proximal repeats as the potential drivers of oncogenesis in mismatch repair defective cells.

Chikungunya virus 3' untranslated region: adaptation to mosquitoes and a population bottleneck as major evolutionary forces

The 3' untranslated genome region (UTR) of arthropod-borne viruses is characterized by enriched direct repeats (DRs) and stem-loop structures. Despite many years of theoretical and experimental study, on-going positive selection on the 3'UTR had never been observed in 'real-time,' and the role of the arbovirus 3'UTR remains poorly understood. We observed a lineage-specific 3'UTR sequence pattern in all available Asian lineage of the mosquito-borne alphavirus, chikungunya virus (CHIKV) (1958-2009), including complicated mutation and duplication patterns of the long DRs. Given that a longer genome is usually associated with less efficient replication, we hypothesized that the fixation of these genetic changes in the Asian lineage 3'UTR was due to their beneficial effects on adaptation to vectors or hosts. Using reverse genetic methods, we examined the functional importance of each direct repeat. Our results suggest that adaptation to mosquitoes, rather than to mammalian hosts, is a major evolutionary force on the CHIKV 3'UTR. Surprisingly, the Asian 3'UTR appeared to be inferior to its predicted ancestral sequence for replication in both mammals and mosquitoes, suggesting that its fixation in Asia was not a result of directional selection. Rather, it may have resulted from a population bottleneck during its introduction from Africa to Asia. We propose that this introduction of a 3'UTR with deletions led to genetic drift and compensatory mutations associated with the loss of structural/functional constraints, followed by two independent beneficial duplications and fixation due to positive selection. Our results provide further evidence that the limited epidemic potential of the Asian CHIKV strains resulted from founder effects that reduced its fitness for efficient transmission by mosquitoes there.

The role of variable DNA tandem repeats in bacterial adaptation

DNA tandem repeats (TRs), also designated as satellite DNA, are inter- or intragenic nucleotide sequences that are repeated two or more times in a head-to-tail manner. Because TR tracts are prone to strand-slippage replication and recombination events that cause the TR copy number to increase or decrease, loci containing TRs are hypermutable. An increasing number of examples illustrate that bacteria can exploit this instability of TRs to reversibly shut down or modulate the function of specific genes, allowing them to adapt to changing environments on short evolutionary time scales without an increased overall mutation rate. In this review, we discuss the prevalence and distribution of inter- and intragenic TRs in bacteria and the mechanisms of their instability. In addition, we review evidence demonstrating a role of TR variations in bacterial adaptation strategies, ranging from immune evasion and tissue tropism to the modulation of environmental stress tolerance. Nevertheless, while bioinformatic analysis reveals that most bacterial genomes contain a few up to several dozens of intra- and intergenic TRs, only a small fraction of these have been functionally studied to date.