Search

Genome-wide association studies (GWAS) have identified hundreds of SNPs responsible for variation in human quantitative traits. However, genome-wide-significant associations often fail to replicate across independent cohorts, in apparent inconsistency with their apparent strong effects in discovery cohorts. This limited success of replication raises pervasive questions about the utility of the GWAS field. We identify all 332 studies of quantitative traits from the NHGRI-EBI GWAS Database with attempted replication. We find that the majority of studies provide insufficient data to evaluate replication rates. The remaining papers replicate significantly worse than expected (p < 10−14), even when adjusting for regression-to-the-mean of effect size between discovery- and replication-cohorts termed the Winner’s Curse (p < 10−16). We show this is due in part to misreporting replication cohort-size as a maximum number, rather than per-locus one. In 39 studies accurately reporting per-locus cohort-size for attempted replication of 707 loci in samples with similar ancestry, replication rate matched expectation (predicted 458, observed 457, p = 0.94). In contrast, ancestry differences between replication and discovery (13 studies, 385 loci) cause the most highly-powered decile of loci to replicate worse than expected, due to difference in linkage disequilibrium. [ABSTRACT FROM AUTHOR]

In recent years, genome and exome sequencing studies have implicated a plethora of new disease genes with rare causal variants. Here, I review 150 exome sequencing studies that claim to have discovered that a disease can be caused by different rare variants in the same gene, and I determine whether their methods followed the current best-practice guidelines in the interpretation of their data. Specifically, I assess whether studies appropriately assess controls for rare variants throughout the entire gene or implicated region as opposed to only investigating the specific rare variants identified in the cases, and I assess whether studies present sufficient co-segregation data for statistically significant linkage. I find that the proportion of studies performing gene-based analyses has increased with time, but that even in 2015 fewer than 40% of the reviewed studies used this method, and only 10% presented statistically significant co-segregation data. Furthermore, I find that the genes reported in these papers are explaining a decreasing proportion of cases as the field moves past most of the low-hanging fruit, with 50% of the genes from studies in 2014 and 2015 having variants in fewer than 5% of cases. As more studies focus on genes explaining relatively few cases, the importance of performing appropriate gene-based analyses is increasing. It is becoming increasingly important for journal editors and reviewers to require stringent gene-based evidence to avoid an avalanche of misleading disease gene discovery papers. [ABSTRACT FROM AUTHOR]

Phosphohexomutase superfamily enzymes catalyze the reversible intramolecular transfer of a phosphoryl moiety on hexose sugars. Bacillus subtilis phosphoglucomutase PgcA catalyzes the reversible interconversion of glucose 6-phosphate (Glc-6-P) and glucose 1-phosphate (Glc-1-P), a precursor of UDP-glucose (UDP-Glc). B. subtilis phosphoglucosamine mutase (GlmM) is a member of the same enzyme superfamily that converts glucosamine 6-phosphate (GlcN-6-P) to glucosamine 1-phosphate (GlcN-1-P), a precursor of the amino sugar moiety of peptidoglycan. Here, we present evidence that B. subtilis PgcA possesses activity as a phosphoglucosamine mutase that contributes to peptidoglycan biosynthesis. This activity was made genetically apparent by the synthetic lethality of pgcA with glmR, a positive regulator of amino sugar biosynthesis, which can be specifically suppressed by overproduction of GlmM. A gain-of-function mutation in a substrate binding loop (PgcA G47S) increases this secondary activity and suppresses a glmR mutant. Our results demonstrate that bacterial phosphoglucomutases may possess secondary phosphoglucosamine mutase activity, and that this dual activity may provide some level of functional redundancy for the essential peptidoglycan biosynthesis pathway. [ABSTRACT FROM AUTHOR]

Antimicrobial resistance (AMR) is currently one of the most important challenges to the treatment of bacterial infections. A critical issue to combat AMR is to restrict its spread. In several instances, bacterial plasmids are involved in the global spread of AMR. Plasmids belonging to the incompatibility group (Inc)HI are widespread in Enterobacteriaceae and most of them express multiple antibiotic resistance determinants. They play a relevant role in the recent spread of colistin resistance. We present in this report novel findings regarding IncHI plasmid conjugation. Conjugative transfer in liquid medium of an IncHI plasmid requires expression of a plasmid-encoded, large-molecular-mass protein that contains an Ig-like domain. The protein, termed RSP, is encoded by a gene (ORF R0009) that maps in the Tra2 region of the IncHI1 R27 plasmid. The RSP protein is exported outside the cell by using the plasmid-encoded type IV secretion system that is also used for its transmission to new cells. Expression of the protein reduces cell motility and enables plasmid conjugation. Flagella are one of the cellular targets of the RSP protein. The RSP protein is required for a high rate of plasmid transfer in both flagellated and nonflagellated Salmonella cells. This effect suggests that RSP interacts with other cellular structures as well as with flagella. These unidentified interactions must facilitate mating pair formation and, hence, facilitate IncHI plasmid conjugation. Due to its location on the outer surfaces of the bacterial cell, targeting the RSP protein could be a means of controlling IncHI plasmid conjugation in natural environments or of combatting infections caused by AMR enterobacteria that harbor IncHI plasmids. [ABSTRACT FROM AUTHOR]

The Y chromosome harbors nine multi-copy ampliconic gene families expressed exclusively in testis. The gene copies within each family are >99% identical to each other, which poses a major challenge in evaluating their copy number. Recent studies demonstrated high variation in Y ampliconic gene copy number among humans. However, how this variation affects expression levels in human testis remains understudied. Here we developed a novel computational tool Ampliconic Copy Number Estimator (AmpliCoNE) that utilizes read sequencing depth information to estimate Y ampliconic gene copy number per family. We applied this tool to whole-genome sequencing data of 149 men with matched testis expression data whose samples are part of the Genotype-Tissue Expression (GTEx) project. We found that the Y ampliconic gene families with low copy number in humans were deleted or pseudogenized in non-human great apes, suggesting relaxation of functional constraints. Among the Y ampliconic gene families, higher copy number leads to higher expression. Within the Y ampliconic gene families, copy number does not influence gene expression, rather a high tolerance for variation in gene expression was observed in testis of presumably healthy men. No differences in gene expression levels were found among major Y haplogroups. Age positively correlated with expression levels of the HSFY and PRY gene families in the African subhaplogroup E1b, but not in the European subhaplogroups R1b and I1. We also found that expression of five out of six Y ampliconic gene families is coordinated with that of their non-Y (i.e. X or autosomal) homologs. Indeed, five ampliconic gene families had consistently lower expression levels when compared to their non-Y homologs suggesting dosage regulation, while the HSFY family had higher expression levels than its X homolog and thus lacked dosage regulation. [ABSTRACT FROM AUTHOR]

There has been much effort to prioritize genomic variants with respect to their impact on “function”. However, function is often not precisely defined: sometimes it is the disease association of a variant; on other occasions, it reflects a molecular effect on transcription or epigenetics. Here, we coupled multiple genomic predictors to build GRAM, a GeneRAlized Model, to predict a well-defined experimental target: the expression-modulating effect of a non-coding variant on its associated gene, in a transferable, cell-specific manner. Firstly, we performed feature engineering: using LASSO, a regularized linear model, we found transcription factor (TF) binding most predictive, especially for TFs that are hubs in the regulatory network; in contrast, evolutionary conservation, a popular feature in many other variant-impact predictors, has almost no contribution. Moreover, TF binding inferred from in vitro SELEX is as effective as that from in vivo ChIP-Seq. Second, we implemented GRAM integrating only SELEX features and expression profiles; thus, the program combines a universal regulatory score with an easily obtainable modifier reflecting the particular cell type. We benchmarked GRAM on large-scale MPRA datasets, achieving AUROC scores of 0.72 in GM12878 and 0.66 in a multi-cell line dataset. We then evaluated the performance of GRAM on targeted regions using luciferase assays in the MCF7 and K562 cell lines. We noted that changing the insertion position of the construct relative to the reporter gene gave very different results, highlighting the importance of carefully defining the exact prediction target of the model. Finally, we illustrated the utility of GRAM in fine-mapping causal variants and developed a practical software pipeline to carry this out. In particular, we demonstrated in specific examples how the pipeline could pinpoint variants that directly modulate gene expression within a larger linkage-disequilibrium block associated with a phenotype of interest (e.g., for an eQTL). [ABSTRACT FROM AUTHOR]

Allele age has long been a focus of population genetic research, primarily because it can be an important clue to the fitness effects of an allele. By virtue of their effects on fitness, alleles under directional selection are expected to be younger than neutral alleles of the same frequency. We developed a new coalescent-based estimator of a close proxy for allele age, the time when a copy of an allele first shares common ancestry with other chromosomes in a sample not carrying that allele. The estimator performs well, including for the very rarest of alleles that occur just once in a sample, with a bias that is typically negative. The estimator is mostly insensitive to population demography and to factors that can arise in population genomic pipelines, including the statistical phasing of chromosomes. Applications to 1000 Genomes Data and UK10K genome data confirm predictions that singleton alleles that alter proteins are significantly younger than those that do not, with a greater difference in the larger UK10K dataset, as expected. The 1000 Genomes populations varied markedly in their distributions for singleton allele ages, suggesting that these distributions can be used to inform models of demographic history, including recent events that are only revealed by their impacts on the ages of very rare alleles. [ABSTRACT FROM AUTHOR]

One of the most frequently mutated proteins in human B-lineage leukemia is the transcription factor PAX5. These mutations often result in partial rather than complete loss of function of the transcription factor. While the functional dose of PAX5 has a clear connection to human malignancy, there is limited evidence for that heterozygote loss of PAX5 have a dramatic effect on the development and function of B-cell progenitors. One possible explanation comes from the finding that PAX5 mutated B-ALL often display complex karyotypes and additional mutations. Thus, PAX5 might be one component of a larger transcription factor network targeted in B-ALL. To investigate the functional network associated with PAX5 we used BioID technology to isolate proteins associated with this transcription factor in the living cell. This identified 239 proteins out of which several could be found mutated in human B-ALL. Most prominently we identified the commonly mutated IKZF1 and RUNX1, involved in the formation of ETV6-AML1 fusion protein, among the interaction partners. ChIP- as well as PLAC-seq analysis supported the idea that these factors share a multitude of target genes in human B-ALL cells. Gene expression analysis of mouse models and primary human leukemia suggested that reduced function of PAX5 increased the ability of an oncogenic form of IKZF1 or ETV6-AML to modulate gene expression. Our data reveals that PAX5 belong to a regulatory network frequently targeted by multiple mutations in B-ALL shedding light on the molecular interplay in leukemia cells. [ABSTRACT FROM AUTHOR]

The SWI/SNF-family chromatin remodeling protein ATRX is a tumor suppressor in sarcomas, gliomas and other malignancies. Its loss of function facilitates the alternative lengthening of telomeres (ALT) pathway in tumor cells, while it also affects Polycomb repressive complex 2 (PRC2) silencing of its target genes. To further define the role of inactivating ATRX mutations in carcinogenesis, we knocked out atrx in our previously reported p53/nf1-deficient zebrafish line that develops malignant peripheral nerve sheath tumors and gliomas. Complete inactivation of atrx using CRISPR/Cas9 was lethal in developing fish and resulted in an alpha-thalassemia-like phenotype including reduced alpha-globin expression. In p53/nf1-deficient zebrafish neither peripheral nerve sheath tumors nor gliomas showed accelerated onset in atrx+/- fish, but these fish developed various tumors that were not observed in their atrx+/+ siblings, including epithelioid sarcoma, angiosarcoma, undifferentiated pleomorphic sarcoma and rare types of carcinoma. These cancer types are included in the AACR Genie database of human tumors associated with mutant ATRX, indicating that our zebrafish model reliably mimics a role for ATRX-loss in the early pathogenesis of these human cancer types. RNA-seq of p53/nf1- and p53/nf1/atrx-deficient tumors revealed that down-regulation of telomerase accompanied ALT-mediated lengthening of the telomeres in atrx-mutant samples. Moreover, inactivating mutations in atrx disturbed PRC2-target gene silencing, indicating a connection between ATRX loss and PRC2 dysfunction in cancer development. [ABSTRACT FROM AUTHOR]

At the molecular level, the evolution of new traits can be broadly divided between changes in gene expression and changes in protein-coding sequence. For proteins, the evolution of novel functions is generally thought to proceed through sequential point mutations or recombination of whole functional units. In Saccharomyces, the uptake of the sugar maltotriose into the cell is the primary limiting factor in its utilization, but maltotriose transporters are relatively rare, except in brewing strains. No known wild strains of Saccharomyces eubayanus, the cold-tolerant parent of hybrid lager-brewing yeasts (Saccharomyces cerevisiae x S. eubayanus), are able to consume maltotriose, which limits their ability to fully ferment malt extract. In one strain of S. eubayanus, we found a gene closely related to a known maltotriose transporter and were able to confer maltotriose consumption by overexpressing this gene or by passaging the strain on maltose. Even so, most wild strains of S. eubayanus lack native maltotriose transporters. To determine how this rare trait could evolve in naive genetic backgrounds, we performed an adaptive evolution experiment for maltotriose consumption, which yielded a single strain of S. eubayanus able to grow on maltotriose. We mapped the causative locus to a gene encoding a novel chimeric transporter that was formed by an ectopic recombination event between two genes encoding transporters that are unable to import maltotriose. In contrast to classic models of the evolution of novel protein functions, the recombination breakpoints occurred within a single functional domain. Thus, the ability of the new protein to carry maltotriose was likely acquired through epistatic interactions between independently evolved substitutions. By acquiring multiple mutations at once, the transporter rapidly gained a novel function, while bypassing potentially deleterious intermediate steps. This study provides an illuminating example of how recombination between paralogs can establish novel interactions among substitutions to create adaptive functions. [ABSTRACT FROM AUTHOR]

Saccharomyces eubayanus is the non-S. cerevisiae parent of the lager-brewing hybrid S. pastorianus. In contrast to most S. cerevisiae and Frohberg-type S. pastorianus strains, S. eubayanus cannot utilize the α-tri-glucoside maltotriose, a major carbohydrate in brewer’s wort. In Saccharomyces yeasts, utilization of maltotriose is encoded by the subtelomeric MAL gene family, and requires transporters for maltotriose uptake. While S. eubayanus strain CBS 12357T harbors four SeMALT genes which enable uptake of the α-di-glucoside maltose, it lacks maltotriose transporter genes. In S. cerevisiae, sequence identity indicates that maltotriose and maltose transporters likely evolved from a shared ancestral gene. To study the evolvability of maltotriose utilization in S. eubayanus CBS 12357T, maltotriose-assimilating mutants obtained after UV mutagenesis were subjected to laboratory evolution in carbon-limited chemostat cultures on maltotriose-enriched wort. An evolved strain showed improved maltose and maltotriose fermentation in 7 L fermenter experiments on industrial wort. Whole-genome sequencing revealed a novel mosaic SeMALT413 gene, resulting from repeated gene introgressions by non-reciprocal translocation of at least three SeMALT genes. The predicted tertiary structure of SeMalT413 was comparable to the original SeMalT transporters, but overexpression of SeMALT413 sufficed to enable growth on maltotriose, indicating gene neofunctionalization had occurred. The mosaic structure of SeMALT413 resembles the structure of S. pastorianus maltotriose-transporter gene SpMTY1, which has high sequences identity to alternatingly S. cerevisiae MALx1, S. paradoxus MALx1 and S. eubayanus SeMALT3. Evolution of the maltotriose transporter landscape in hybrid S. pastorianus lager-brewing strains is therefore likely to have involved mechanisms similar to those observed in the present study. [ABSTRACT FROM AUTHOR]

The author discusses the study which describes the RAS superfamily conditions RALA involving GTPase function, Details about the study which emphasized genotype-phenotype correlations in RAS-related genes and RALA could involve additional mechanisms and pathways are offered. The author believes that the scope of developmental effects should be examined with germline RALA and ARF1 genes.

Complex chromosomal rearrangements (CCRs) are rearrangements involving more than two chromosomes or more than two breakpoints. Whole genome sequencing (WGS) allows for outstanding high resolution characterization on the nucleotide level in unique sequences of such rearrangements, but problems remain for mapping breakpoints in repetitive regions of the genome, which are known to be prone to rearrangements. Hence, multiple complementary WGS experiments are sometimes needed to solve the structures of CCRs. We have studied three individuals with CCRs: Case 1 and Case 2 presented with de novo karyotypically balanced, complex interchromosomal rearrangements (46,XX,t(2;8;15)(q35;q24.1;q22) and 46,XY,t(1;10;5)(q32;p12;q31)), and Case 3 presented with a de novo, extremely complex intrachromosomal rearrangement on chromosome 1. Molecular cytogenetic investigation revealed cryptic deletions in the breakpoints of chromosome 2 and 8 in Case 1, and on chromosome 10 in Case 2, explaining their clinical symptoms. In Case 3, 26 breakpoints were identified using WGS, disrupting five known disease genes. All rearrangements were subsequently analyzed using optical maps, linked-read WGS, and short-read WGS. In conclusion, we present a case series of three unique de novo CCRs where we by combining the results from the different technologies fully solved the structure of each rearrangement. The power in combining short-read WGS with long-molecule sequencing or optical mapping in these unique de novo CCRs in a clinical setting is demonstrated. [ABSTRACT FROM AUTHOR]

More than any other genome components, Transposable Elements (TEs) have the capacity to move across species barriers through Horizontal Transfer (HT), with substantial evolutionary consequences. Previous large-scale surveys, based on full-genomes comparisons, have revealed the transposition mode as an important predictor of HT rates variation across TE superfamilies. However, host biology could represent another major explanatory factor, one that needs to be investigated through extensive taxonomic sampling. Here we test this hypothesis using a field collection of 460 Arthropod species from Tahiti and surrounding islands. Through targeted massive parallel sequencing, we uncover patterns of HT in three widely-distributed TE superfamilies with contrasted modes of transposition. In line with earlier findings, the DNA transposons under study (TC1-Mariner) were found to transfer horizontally at the highest frequency, closely followed by the LTR superfamily (Copia), in contrast with the non-LTR superfamily (Jockey), that mostly diversifies through vertical inheritance and persists longer within genomes. Strikingly, across all superfamilies, we observe a marked excess of HTs in Lepidoptera, an insect order that also commonly hosts baculoviruses, known for their ability to transport host TEs. These results turn the spotlight on baculoviruses as major potential vectors of TEs in Arthropods, and further emphasize the importance of non-vertical TE inheritance in genome evolution. [ABSTRACT FROM AUTHOR]

The BRCA Challenge is a long-term data-sharing project initiated within the Global Alliance for Genomics and Health (GA4GH) to aggregate BRCA1 and BRCA2 data to support highly collaborative research activities. Its goal is to generate an informed and current understanding of the impact of genetic variation on cancer risk across the iconic cancer predisposition genes, BRCA1 and BRCA2. Initially, reported variants in BRCA1 and BRCA2 available from public databases were integrated into a single, newly created site, . The purpose of the BRCA Exchange is to provide the community with a reliable and easily accessible record of variants interpreted for a high-penetrance phenotype. More than 20,000 variants have been aggregated, three times the number found in the next-largest public database at the project’s outset, of which approximately 6,150 have expert classifications. The data set is based on shared information from existing clinical databases—Breast Cancer Information Core (BIC), ClinVar, and the Leiden Open Variation Database (LOVD)—as well as population databases, all linked to a single point of access. The BRCA Challenge has brought together the existing international Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium expert panel, along with expert clinicians, diagnosticians, researchers, and database providers, all with a common goal of advancing our understanding of BRCA1 and BRCA2 variation. Ongoing work includes direct contact with national centers with access to BRCA1 and BRCA2 diagnostic data to encourage data sharing, development of methods suitable for extraction of genetic variation at the level of individual laboratory reports, and engagement with participant communities to enable a more comprehensive understanding of the clinical significance of genetic variation in BRCA1 and BRCA2. [ABSTRACT FROM AUTHOR]

Sequencing of whole cancer genomes has revealed an abundance of recurrent mutations in gene-regulatory promoter regions, in particular in melanoma where strong mutation hotspots are observed adjacent to ETS-family transcription factor (TF) binding sites. While sometimes interpreted as functional driver events, these mutations are commonly believed to be due to locally inhibited DNA repair. Here, we first show that low-dose UV light induces mutations preferably at a known ETS promoter hotspot in cultured cells even in the absence of global or transcription-coupled nucleotide excision repair (NER). Further, by genome-wide mapping of cyclobutane pyrimidine dimers (CPDs) shortly after UV exposure and thus before DNA repair, we find that ETS-related mutation hotspots exhibit strong increases in CPD formation efficacy in a manner consistent with tumor mutation data at the single-base level. Analysis of a large whole genome cohort illustrates the widespread contribution of this effect to recurrent mutations in melanoma. While inhibited NER underlies a general increase in somatic mutation burden in regulatory elements including ETS sites, our data supports that elevated DNA damage formation at specific genomic bases is at the core of the prominent promoter mutation hotspots seen in skin cancers, thus explaining a key phenomenon in whole-genome cancer analyses. [ABSTRACT FROM AUTHOR]

Mitochondrial DNA (mtDNA) is useful to assist with identification of the source of a biological sample, or to confirm matrilineal relatedness. Although the autosomal genome is much larger, mtDNA has an advantage for forensic applications of multiple copy number per cell, allowing better recovery of sequence information from degraded samples. In addition, biological samples such as fingernails, old bones, teeth and hair have mtDNA but little or no autosomal DNA. The relatively low mutation rate of the mitochondrial genome (mitogenome) means that there can be large sets of matrilineal-related individuals sharing a common mitogenome. Here we present the simulation software that we use to describe the distribution of the number of mitogenomes in a population that match a given mitogenome, and investigate its dependence on population size and growth rate, and on a database count of the mitogenome. Further, we report on the distribution of the number of meioses separating pairs of individuals with matching mitogenome. Our results have important implications for assessing the weight of mtDNA profile evidence in forensic science, but mtDNA analysis has many non-human applications, for example in tracking the source of ivory. Our methods and software can also be used for simulations to help validate models of population history in human or non-human populations. [ABSTRACT FROM AUTHOR]

Morphogen gradients provide concentration-dependent positional information along polarity axes. Although the dynamics of the establishment of these gradients is well described, precision and noise in the downstream activation processes remain elusive. A simple paradigm to address these questions is the Bicoid morphogen gradient that elicits a rapid step-like transcriptional response in young fruit fly embryos. Focusing on the expression of the major Bicoid target, hunchback (hb), at the onset of zygotic transcription, we used the MS2-MCP approach which combines fluorescent labeling of nascent mRNA with live imaging at high spatial and temporal resolution. Removing 36 putative Zelda binding sites unexpectedly present in the original MS2 reporter, we show that the 750 bp of the hb promoter are sufficient to recapitulate endogenous expression at the onset of zygotic transcription. After each mitosis, in the anterior, expression is turned on to rapidly reach a plateau with all nuclei expressing the reporter. Consistent with a Bicoid dose-dependent activation process, the time period required to reach the plateau increases with the distance to the anterior pole. Despite the challenge imposed by frequent mitoses and high nuclei-to-nuclei variability in transcription kinetics, it only takes 3 minutes at each interphase for the MS2 reporter loci to distinguish subtle differences in Bicoid concentration and establish a steadily positioned and steep (Hill coefficient ~ 7) expression boundary. Modeling based on the cooperativity between the 6 known Bicoid binding sites in the hb promoter region, assuming rate limiting concentrations of the Bicoid transcription factor at the boundary, is able to capture the observed dynamics of pattern establishment but not the steepness of the boundary. This suggests that a simple model based only on the cooperative binding of Bicoid is not sufficient to describe the spatiotemporal dynamics of early hb expression. [ABSTRACT FROM AUTHOR]

Structural features of genomes, including the three-dimensional arrangement of DNA in the nucleus, are increasingly seen as key contributors to the regulation of gene expression. However, studies on how genome structure and nuclear organisation influence transcription have so far been limited to a handful of model species. This narrow focus limits our ability to draw general conclusions about the ways in which three-dimensional structures are encoded, and to integrate information from three-dimensional data to address a broader gamut of biological questions. Here, we generate a complete and gapless genome sequence for the filamentous fungus, Epichloë festucae. We use Hi-C data to examine the three-dimensional organisation of the genome, and RNA-seq data to investigate how Epichloë genome structure contributes to the suite of transcriptional changes needed to maintain symbiotic relationships with the grass host. Our results reveal a genome in which very repeat-rich blocks of DNA with discrete boundaries are interspersed by gene-rich sequences that are almost repeat-free. In contrast to other species reported to date, the three-dimensional structure of the genome is anchored by these repeat blocks, which act to isolate transcription in neighbouring gene-rich regions. Genes that are differentially expressed in planta are enriched near the boundaries of these repeat-rich blocks, suggesting that their three-dimensional orientation partly encodes and regulates the symbiotic relationship formed by this organism. [ABSTRACT FROM AUTHOR]

Short insertions, deletions (INDELs) and larger structural variants have been increasingly employed in genetic association studies, but few improvements over SNP-based association have been reported. In order to understand why this might be the case, we analysed two publicly available datasets and observed that 63% of INDELs called in A. thaliana and 64% in D. melanogaster populations are misrepresented as multiple alleles with different functional annotations, i.e. where the same underlying variant is represented by inconsistent alignments leading to different variant calls. To address this issue, we have developed the software Irisas to reclassify and re-annotate these variants, which we then used for single-locus tests of association. We also integrated them to predict the functional impact of SNPs, INDELs, and structural variants for burden testing. Using both approaches, we re-analysed the genetic architecture of complex traits in A. thaliana and D. melanogaster. Heritability analysis using SNPs alone explained on average 27% and 19% of phenotypic variance for A. thaliana and D. melanogaster respectively. Our method explained an additional 11% and 3%, respectively. We also identified novel trait loci that previous SNP-based association studies failed to map, and which contain established candidate genes. Our study shows the value of the association test with INDELs and integrating multiple types of variants in association studies in plants and animals. [ABSTRACT FROM AUTHOR]

Synonymous mutations continue to be filtered out from most large-scale cancer genome studies, but several lines of evidence suggest they can play driver roles in neoplastic disease. We investigated a case of an aggressive, apparently sporadic medullary thyroid carcinoma (MTC) harboring a somatic RET p.Cys634Arg mutation (a known MTC driver). A germ-line RET substitution (p.Cys630=) had also been found but was considered clinically irrelevant because of its synonymous nature. Next generation sequencing (NGS) of the tumor tissues revealed that the RET mutations were in cis. There was no evidence of gene amplification. Expression analysis found an increase of RET transcript in p.Cys630=;p.Cys634Arg patient compared with that found in 7 MTCs harboring p.Cys634 mutations. Minigene expression assays demonstrated that the presence of the synonymous RET mutation was sufficient to explain the increased RET mRNA level. In silico analyses and RNA immunoprecipitation experiments showed that the p.Cys630 = variant created new exonic splicing enhancer motifs that enhanced SRp55 recruitment to the mutant allele, leading to more efficient maturation of its pre-mRNA and an increased abundance of mature mRNA encoding a constitutively active RET receptor. These findings document a novel mechanism by which synonymous mutations can contribute to cancer progression. [ABSTRACT FROM AUTHOR]

Post-mitotic cell separation is one of the most prominent events in the life cycle of eukaryotic cells, but the molecular underpinning of this fundamental biological process is far from being concluded and fully characterized. We use budding yeast Saccharomyces cerevisiae as a model and demonstrate AMN1 as a major gene underlying post-mitotic cell separation in a natural yeast strain, YL1C. Specifically, we define a novel 11-residue domain by which Amn1 binds to Ace2. Moreover, we demonstrate that Amn1 induces proteolysis of Ace2 through the ubiquitin proteasome system and in turn, down-regulates Ace2’s downstream target genes involved in hydrolysis of the primary septum, thus leading to inhibition of cell separation and clumping of haploid yeast cells. Using ChIP assays and site-specific mutation experiments, we show that Ste12 and the a1-α12 heterodimer are two direct regulators of AMN1. Specifically, a1-α2, a diploid-specific heterodimer, prevents Ste12 from inactivating AMN1 through binding to its promoter. This demonstrates how the Amn1-governed cell separation is highly cell type dependent. Finally, we show that AMN1368D from YL1C is a dominant allele in most strains of S. cerevisiae and evolutionarily conserved in both genic structure and phenotypic effect in two closely related yeast species, K. lactis and C. glabrata. [ABSTRACT FROM AUTHOR]

Recursive splicing (RS) is an evolutionarily conserved process of removing long introns via multiple steps of splicing. It was first discovered in Drosophila and recently proven to occur also in humans. The detailed mechanism of recursive splicing is not well understood, in particular, whether it is kinetically coupled with transcription. To investigate the dynamic process that underlies recursive splicing, we systematically characterized 342 RS sites in three human cell types using published time-series data that monitored synchronized Pol II elongation and nascent RNA production with 4-thiouridine labeling. We found that half of the RS events occurred post-transcriptionally with long delays. For at least 18–47% RS introns, we detected RS junction reads only after detecting canonical splicing junction reads, supporting the notion that these introns were removed by both recursive splicing and canonical splicing. Furthermore, the choice of which splicing mechanism was used showed cell type specificity. Our results suggest that recursive splicing supplements, rather than replaces, canonical splicing for removing long introns. [ABSTRACT FROM AUTHOR]

The clinical spectrum of ciliopathies affecting motile cilia spans impaired mucociliary clearance in the respiratory system, laterality defects including heart malformations, infertility and hydrocephalus. Using linkage analysis and whole exome sequencing, we identified two recessive loss-of-function MNS1 mutations in five individuals from four consanguineous families: 1) a homozygous nonsense mutation p.Arg242* in four males with laterality defects and infertility and 2) a homozygous nonsense mutation p.Gln203* in one female with laterality defects and recurrent respiratory infections additionally carrying homozygous mutations in DNAH5. Consistent with the laterality defects observed in these individuals, we found Mns1 to be expressed in mouse embryonic ventral node. Immunofluorescence analysis further revealed that MNS1 localizes to the axonemes of respiratory cilia as well as sperm flagella in human. In-depth ultrastructural analyses confirmed a subtle outer dynein arm (ODA) defect in the axonemes of respiratory epithelial cells resembling findings reported in Mns1-deficient mice. Ultrastructural analyses in the female carrying combined mutations in MNS1 and DNAH5 indicated a role for MNS1 in the process of ODA docking (ODA-DC) in the distal respiratory axonemes. Furthermore, co-immunoprecipitation and yeast two hybrid analyses demonstrated that MNS1 dimerizes and interacts with the ODA docking complex component CCDC114. Overall, we demonstrate that MNS1 deficiency in humans causes laterality defects (situs inversus) and likely male infertility and that MNS1 plays a role in the ODA-DC assembly. [ABSTRACT FROM AUTHOR]

Cellular responses to changing environments frequently involve rapid reprogramming of the transcriptome. Regulated changes in mRNA degradation rates can accelerate reprogramming by clearing or stabilizing extant transcripts. Here, we measured mRNA stability using 4-thiouracil labeling in the budding yeast Saccharomyces cerevisiae during a nitrogen upshift and found that 78 mRNAs are subject to destabilization. These transcripts include Nitrogen Catabolite Repression (NCR) and carbon metabolism mRNAs, suggesting that mRNA destabilization is a mechanism for targeted reprogramming of the transcriptome. To explore the molecular basis of destabilization we implemented a SortSeq approach to screen the pooled deletion collection library for trans factors that mediate rapid GAP1 mRNA repression. We combined low-input multiplexed arcode sequencing with branched-DNA single-molecule mRNA ISH and luorescence-activated cell sorting (BFF) to identify the Lsm1-7p/Pat1p complex and general mRNA decay machinery as important for GAP1 mRNA clearance. We also find that the decapping modulators EDC3 and SCD6, translation factor eIF4G2, and the 5’ UTR of GAP1 are factors that mediate rapid repression of GAP1 mRNA, suggesting that translational control may impact the post-transcriptional fate of mRNAs in response to environmental changes. [ABSTRACT FROM AUTHOR]

Cyanobacteria are a monophyletic phylogenetic group of global importance and have received considerable attention as potential host organisms for the renewable synthesis of chemical bulk products from atmospheric CO2. The cyanobacterial phylum exhibits enormous metabolic diversity with respect to morphology, lifestyle and habitat. As yet, however, research has mostly focused on few model strains and cyanobacterial diversity is insufficiently understood. In this respect, the increasing availability of fully sequenced bacterial genomes opens new and unprecedented opportunities to investigate the genetic inventory of organisms in the context of their pan-genome. Here, we seek understand cyanobacterial diversity using a comparative genome analysis of 77 fully sequenced and assembled cyanobacterial genomes. We use phylogenetic profiling to analyze the co-occurrence of clusters of likely ortholog genes (CLOGs) and reveal novel functional associations between CLOGs that are not captured by co-localization of genes. Going beyond pair-wise co-occurrences, we propose a network approach that allows us to identify modules of co-occurring CLOGs. The extracted modules exhibit a high degree of functional coherence and reveal known as well as previously unknown functional associations. We argue that the high functional coherence observed for the modules is a consequence of the similar-yet-diverse nature of cyanobacteria. Our approach highlights the importance of a multi-strain analysis to understand gene functions and environmental adaptations, with implications beyond the cyanobacterial phylum. The analysis is augmented with a simple toolbox that facilitates further analysis to investigate the co-occurrence neighborhood of specific CLOGs of interest. [ABSTRACT FROM AUTHOR]

House dust mites are common pests with an unusual evolutionary history, being descendants of a parasitic ancestor. Transition to parasitism is frequently accompanied by genome rearrangements, possibly to accommodate the genetic change needed to access new ecology. Transposable element (TE) activity is a source of genomic instability that can trigger large-scale genomic alterations. Eukaryotes have multiple transposon control mechanisms, one of which is RNA interference (RNAi). Investigation of the dust mite genome failed to identify a major RNAi pathway: the Piwi-associated RNA (piRNA) pathway, which has been replaced by a novel small-interfering RNAs (siRNAs)-like pathway. Co-opting of piRNA function by dust mite siRNAs is extensive, including establishment of TE control master loci that produce siRNAs. Interestingly, other members of the Acari have piRNAs indicating loss of this mechanism in dust mites is a recent event. Flux of RNAi-mediated control of TEs highlights the unusual arc of dust mite evolution. [ABSTRACT FROM AUTHOR]

Nephron progenitor number determines nephron endowment; a reduced nephron count is linked to the onset of kidney disease. Several transcriptional regulators including Six2, Wt1, Osr1, Sall1, Eya1, Pax2, and Hox11 paralogues are required for specification and/or maintenance of nephron progenitors. However, little is known about the regulatory intersection of these players. Here, we have mapped nephron progenitor-specific transcriptional networks of Six2, Hoxd11, Osr1, and Wt1. We identified 373 multi-factor associated ‘regulatory hotspots’ around genes closely associated with progenitor programs. To examine their functional significance, we deleted ‘hotspot’ enhancer elements for Six2 and Wnt4. Removal of the distal enhancer for Six2 leads to a ~40% reduction in Six2 expression. When combined with a Six2 null allele, progeny display a premature depletion of nephron progenitors. Loss of the Wnt4 enhancer led to a significant reduction of Wnt4 expression in renal vesicles and a mildly hypoplastic kidney, a phenotype also enhanced in combination with a Wnt4 null mutation. To explore the regulatory landscape that supports proper target gene expression, we performed CTCF ChIP-seq to identify insulator-boundary regions. One such putative boundary lies between the Six2 and Six3 loci. Evidence for the functional significance of this boundary was obtained by deep sequencing of the radiation-induced Brachyrrhine (Br) mutant allele. We identified an inversion of the Six2/Six3 locus around the CTCF-bound boundary, removing Six2 from its distal enhancer regulation, but placed next to Six3 enhancer elements which support ectopic Six2 expression in the lens where Six3 is normally expressed. Six3 is now predicted to fall under control of the Six2 distal enhancer. Consistent with this view, we observed ectopic Six3 in nephron progenitors. 4C-seq supports the model for Six2 distal enhancer interactions in wild-type and Br/+ mouse kidneys. Together, these data expand our view of the regulatory genome and regulatory landscape underpinning mammalian nephrogenesis. [ABSTRACT FROM AUTHOR]

Normal development requires the right splice variants to be made in the right tissues at the right time. The core splicing machinery is engaged in all splicing events, but which precise splice variant is made requires the choice between alternative splice sites—for this to occur, a set of splicing factors (SFs) must recognize and bind to short RNA motifs in the pre-mRNA. In C. elegans, there is known to be extensive variation in splicing patterns across development, but little is known about the targets of each SF or how multiple SFs combine to regulate splicing. Here we combine RNA-seq with in vitro binding assays to study how 4 different C. elegans SFs, ASD-1, FOX-1, MEC-8, and EXC-7, regulate splicing. The 4 SFs chosen all have well-characterised biology and well-studied loss-of-function genetic alleles, and all contain RRM domains. Intriguingly, while the SFs we examined have varied roles in C. elegans development, they show an unexpectedly high overlap in their targets. We also find that binding sites for these SFs occur on the same pre-mRNAs more frequently than expected suggesting extensive combinatorial control of splicing. We confirm that regulation of splicing by multiple SFs is often combinatorial and show that this is functionally significant. We also find that SFs appear to combine to affect splicing in two modes—they either bind in close proximity within the same intron or they appear to bind to separate regions of the intron in a conserved order. Finally, we find that the genes whose splicing are regulated by multiple SFs are highly enriched for genes involved in the cytoskeleton and in ion channels that are key for neurotransmission. Together, this shows that specific classes of genes have complex combinatorial regulation of splicing and that this combinatorial regulation is critical for normal development to occur. [ABSTRACT FROM AUTHOR]

The article focuses on the initiative "Up For A Challenge (U4C)—Stimulating Innovation in Breast Cancer Genetic Epidemiology" by the American organization National Cancer Institute (NCI). Topics discussed include understanding the role of disease risk through genome-wide association studies (GWAS) of the genetic variants of breast cancer, review by National Institutes of Health (NIH) of entries, and chart on U4C entries received.

While traditional forensic genetics has been oriented towards using human DNA in criminal investigation and civil court cases, it currently presents a much wider application range, including not only legal situations sensu stricto but also and, increasingly often, to preemptively avoid judicial processes. Despite some difficulties, current forensic genetics is progressively incorporating the analysis of nonhuman genetic material to a greater extent. The analysis of this material—including other animal species, plants, or microorganisms—is now broadly used, providing ancillary evidence in criminalistics in cases such as animal attacks, trafficking of species, bioterrorism and biocrimes, and identification of fraudulent food composition, among many others. Here, we explore how nonhuman forensic genetics is being revolutionized by the increasing variety of genetic markers, the establishment of faster, less error-burdened and cheaper sequencing technologies, and the emergence and improvement of models, methods, and bioinformatics facilities. [ABSTRACT FROM AUTHOR]

Box C/D snoRNAs are known to guide site-specific ribose methylation of ribosomal RNA. Here, we demonstrate a novel and unexpected role for box C/D snoRNAs in guiding 18S rRNA acetylation in yeast. Our results demonstrate, for the first time, that the acetylation of two cytosine residues in 18S rRNA catalyzed by Kre33 are guided by two orphan box C/D snoRNAs–snR4 and snR45 –not known to be involved in methylation in yeast. We identified Kre33 binding sites on these snoRNAs as well as on the 18S rRNA, and demonstrate that both snR4 and snR45 establish extended bipartite complementarity around the cytosines targeted for acetylation, similar to pseudouridylation pocket formation by the H/ACA snoRNPs. We show that base pairing between these snoRNAs and 18S rRNA require the putative helicase activity of Kre33, which is also needed to aid early pre-rRNA processing. Compared to yeast, the number of orphan box C/D snoRNAs in higher eukaryotes is much larger and we hypothesize that several of these may be involved in base-modifications. [ABSTRACT FROM AUTHOR]

In many animal embryos, a specific gene expression pattern is established along the animal-vegetal axis soon after zygotic transcription begins. In the embryo of the ascidian Ciona intestinalis, soon after the division that separates animal and vegetal hemispheres into distinct blastomeres, maternal Gata.a and β-catenin activate specific genes in the animal and vegetal blastomeres, respectively. On the basis of these initial distinct gene expression patterns, gene regulatory networks promote animal cells to become ectodermal tissues and vegetal cells to become endomesodermal tissues and a part of the nerve cord. In the vegetal hemisphere, β-catenin directly activates Foxd, an essential transcription factor gene for specifying endomesodermal fates. In the present study, we found that Foxd also represses the expression of genes that are activated specifically in the animal hemisphere, including Dmrt1, Prdm1-r.a (Bz1), Prdm1-r.b (Bz2), and Otx. A reporter assay showed that Dmrt1 expression was directly repressed by Foxd, and a chromatin immunoprecipitation assay showed that Foxd was bound to the upstream regions of Dmrt1, Prdm1-r.a, Prdm1-r.b, and Otx. Thus, Foxd has a dual function of activating specific gene expression in the vegetal hemisphere and of repressing the expression of genes that are normally expressed in the animal hemisphere. This dual function stabilizes the initial patterning along the animal-vegetal axis by β-catenin and Gata.a. [ABSTRACT FROM AUTHOR]

Sequencing of whole tumor genomes holds the promise of revealing functional somatic regulatory mutations, such as those described in the TERT promoter. Recurrent promoter mutations have been identified in many additional genes and appear to be particularly common in melanoma, but convincing functional data such as influence on gene expression has been more elusive. Here, we show that frequently recurring promoter mutations in melanoma occur almost exclusively at cytosines flanked by a distinct sequence signature, TTCCG, with TERT as a notable exception. In active, but not inactive, promoters, mutation frequencies for cytosines at the 5’ end of this ETS-like motif were considerably higher than expected based on a UV trinucleotide mutational signature. Additional analyses solidify this pattern as an extended context-specific mutational signature that mediates an exceptional position-specific vulnerability to UV mutagenesis, arguing against positive selection. We further use ultra-sensitive amplicon sequencing to demonstrate that cell cultures exposed to UV light quickly develop subclonal mutations specifically in affected positions. Our findings have implications for the interpretation of somatic mutations in regulatory regions, and underscore the importance of genomic context and extended sequence patterns to accurately describe mutational signatures in cancer. [ABSTRACT FROM AUTHOR]

Preterm birth (PTB) affects ~12% of pregnancies in the US. Despite its high mortality and morbidity, the molecular etiology underlying PTB has been unclear. Numerous studies have been devoted to identifying genetic factors in maternal and fetal genomes, but so far few genomic loci have been associated with PTB. By analyzing whole-genome sequencing data from 816 trio families, for the first time, we observed the role of fetal de novo mutations in PTB. We observed a significant increase in de novo mutation burden in PTB fetal genomes. Our genomic analyses further revealed that affected genes by PTB de novo mutations were dosage sensitive, intolerant to genomic deletions, and their mouse orthologs were likely developmentally essential. These genes were significantly involved in early fetal brain development, which was further supported by our analysis of copy number variants identified from an independent PTB cohort. Our study indicates a new mechanism in PTB occurrence independently contributed from fetal genomes, and thus opens a new avenue for future PTB research. [ABSTRACT FROM AUTHOR]

Calcineurin is a highly conserved Ca2+/calmodulin-dependent serine/threonine-specific protein phosphatase that orchestrates cellular Ca2+ signaling responses. In Cryptococcus neoformans, calcineurin is activated by multiple stresses including high temperature, and is essential for stress adaptation and virulence. The transcription factor Crz1 is a major calcineurin effector in Saccharomyces cerevisiae and other fungi. Calcineurin dephosphorylates Crz1, thereby enabling Crz1 nuclear translocation and transcription of target genes. Here we show that loss of Crz1 confers phenotypes intermediate between wild-type and calcineurin mutants, and demonstrate that deletion of the calcineurin docking domain results in the inability of Crz1 to translocate into the nucleus under thermal stress. RNA-sequencing revealed 102 genes that are regulated in a calcineurin-Crz1-dependent manner at 37°C. The majority of genes were down-regulated in cna1Δ and crz1Δ mutants, indicating these genes are normally activated by the calcineurin-Crz1 pathway at high temperature. About 58% of calcineurin-Crz1 target genes have unknown functions, while genes with known or predicted functions are involved in cell wall remodeling, calcium transport, and pheromone production. We identified three calcineurin-dependent response element motifs within the promoter regions of calcineurin-Crz1 target genes, and show that Crz1 binding to target gene promoters is increased upon thermal stress in a calcineurin-dependent fashion. Additionally, we found a large set of genes independently regulated by calcineurin, and Crz1 regulates 59 genes independently of calcineurin. Given the intermediate crz1Δ mutant phenotype, and our recent evidence for a calcineurin regulatory network impacting mRNA in P-bodies and stress granules independently of Crz1, calcineurin likely acts on factors beyond Crz1 that govern mRNA expression/stability to operate a branched transcriptional/post-transcriptional stress response network necessary for fungal virulence. Taken together, our findings reveal the core calcineurin-Crz1 stress response cascade is maintained from ascomycetes to a pathogenic basidiomycete fungus, but its output in C. neoformans appears to be adapted to promote fungal virulence. [ABSTRACT FROM AUTHOR]

Structural Maintenance of Chromosomes (SMC) family proteins participate in multisubunit complexes that govern chromosome structure and dynamics. SMC-containing condensin complexes create chromosome topologies essential for mitosis/meiosis, gene expression, recombination, and repair. Many eukaryotes have two condensin complexes (I and II); C. elegans has three (I, II, and the X-chromosome specialized condensin IDC) and their regulation is poorly understood. Here we identify a novel SMC-like protein, SMCL-1, that binds to C. elegans condensin SMC subunits, and modulates condensin functions. Consistent with a possible role as a negative regulator, loss of SMCL-1 partially rescued the lethal and sterile phenotypes of a hypomorphic condensin mutant, while over-expression of SMCL-1 caused lethality, chromosome mis-segregation, and disruption of condensin IDC localization on X chromosomes. Unlike canonical SMC proteins, SMCL-1 lacks hinge and coil domains, and its ATPase domain lacks conserved amino acids required for ATP hydrolysis, leading to the speculation that it may inhibit condensin ATPase activity. SMCL-1 homologs are apparent only in the subset of Caenorhabditis species in which the condensin I and II subunit SMC-4 duplicated to create the condensin IDC- specific subunit DPY-27, suggesting that SMCL-1 helps this lineage cope with the regulatory challenges imposed by evolution of a third condensin complex. Our findings uncover a new regulator of condensins and highlight how the duplication and divergence of SMC complex components in various lineages has created new proteins with diverse functions in chromosome dynamics. [ABSTRACT FROM AUTHOR]

Cooperative transcription factor binding at cis-regulatory sites in the genome drives robust eukaryotic gene expression, and many such sites must be coordinated to produce coherent transcriptional programs. The transcriptional program leading to motile cilia formation requires members of the DNA-binding forkhead (Fox) and Rfx transcription factor families and these factors co-localize to cilia gene promoters, but it is not clear how many cilia genes are regulated by these two factors, whether these factors act directly or indirectly, or how these factors act with specificity in the context of a 3-dimensional genome. Here, we use genome-wide approaches to show that cilia genes reside at the boundaries of topological domains and that these areas have low enhancer density. We show that the transcription factors Foxj1 and Rfx2 binding occurs in the promoters of more cilia genes than other known cilia transcription factors and that while Rfx2 binds directly to promoters and enhancers equally, Foxj1 prefers direct binding to enhancers and is stabilized at promoters by Rfx2. Finally, we show that Rfx2 and Foxj1 lie at the anchor endpoints of chromatin loops, suggesting that target genes are activated when Foxj1 bound at distal sites is recruited via a loop created by Rfx2 binding at both sites. We speculate that the primary function of Rfx2 is to stabilize distal enhancers with proximal promoters by operating as a scaffolding factor, bringing key regulatory domains bound by Foxj1 into close physical proximity and enabling coordinated cilia gene expression. [ABSTRACT FROM AUTHOR]

Jujube (Ziziphus jujuba Mill.) belongs to the Rhamnaceae family and is a popular fruit tree species with immense economic and nutritional value. Here, we report a draft genome of the dry jujube cultivar ‘Junzao’ and the genome resequencing of 31 geographically diverse accessions of cultivated and wild jujubes (Ziziphus jujuba var. spinosa). Comparative analysis revealed that the genome of ‘Dongzao’, a fresh jujube, was ~86.5 Mb larger than that of the ‘Junzao’, partially due to the recent insertions of transposable elements in the ‘Dongzao’ genome. We constructed eight proto-chromosomes of the common ancestor of Rhamnaceae and Rosaceae, two sister families in the order Rosales, and elucidated the evolutionary processes that have shaped the genome structures of modern jujubes. Population structure analysis revealed the complex genetic background of jujubes resulting from extensive hybridizations between jujube and its wild relatives. Notably, several key genes that control fruit organic acid metabolism and sugar content were identified in the selective sweep regions. We also identified S-locus genes controlling gametophytic self-incompatibility and investigated haplotype patterns of the S locus in the jujube genomes, which would provide a guideline for parent selection for jujube crossbreeding. This study provides valuable genomic resources for jujube improvement, and offers insights into jujube genome evolution and its population structure and domestication. [ABSTRACT FROM AUTHOR]

DNA methylation in plants is traditionally partitioned into CG, CHG and CHH contexts (with H any nucleotide but G). By investigating DNA methylation patterns in trinucleotide contexts in four angiosperm species, we show that such a representation hides spatial and functional partitioning of different methylation pathways and is incomplete. CG methylation (mCG) is largely context-independent whereas, at CHG motifs, there is under-representation of mCCG in pericentric regions of A. thaliana and tomato and throughout the chromosomes of maize and rice. In A. thaliana the biased representation of mCCG in heterochromatin is related to specificities of H3K9 methyltransferase SUVH family members. At CHH motifs there is an over-representation of different variant forms of mCHH that, similarly to mCCG hypomethylation, is partitioned into the pericentric regions of the two dicots but dispersed in the monocot chromosomes. The over-represented mCHH motifs in A. thaliana associate with specific types of transposon including both class I and II elements. At mCHH the contextual bias is due to the involvement of various chromomethyltransferases whereas the context-independent CHH methylation in A. thaliana and tomato is mediated by the RNA-directed DNA methylation process that is most active in the gene-rich euchromatin. This analysis therefore reveals that the sequence context of the methylome of plant genomes is informative about the mechanisms associated with maintenance of methylation and the overlying chromatin structure. [ABSTRACT FROM AUTHOR]

The maintenance of chromosome integrity is crucial for genetic stability. However, programmed chromosome fragmentations are known to occur in many organisms, and in the ciliate Tetrahymena the five germline chromosomes are fragmented into hundreds of minichromosomes during somatic nuclear differentiation. Here, we showed that there are different fates of these minichromosomes after chromosome breakage. Among the 326 somatic minichromosomes identified using genomic data, 50 are selectively eliminated from the mature somatic genome. Interestingly, many and probably most of these minichromosomes are eliminated during the growth period between 6 and 20 doublings right after conjugation. Genes with potential conjugation-specific functions are found in these minichromosomes. This study revealed a new mode of programmed DNA elimination in ciliates similar to those observed in parasitic nematodes, which could play a role in developmental gene regulation. [ABSTRACT FROM AUTHOR]

Assembly of kinetochore complexes, involving greater than one hundred proteins, is essential for chromosome segregation and genome stability. Neocentromeres, or new centromeres, occur when kinetochores assemble de novo, at DNA loci not previously associated with kinetochore proteins, and they restore chromosome segregation to chromosomes lacking a functional centromere. Neocentromeres have been observed in a number of diseases and may play an evolutionary role in adaptation or speciation. However, the consequences of neocentromere formation on chromosome missegregation rates, gene expression, and three-dimensional (3D) nuclear structure are not well understood. Here, we used Candida albicans, an organism with small, epigenetically-inherited centromeres, as a model system to study the functions of twenty different neocentromere loci along a single chromosome, chromosome 5. Comparison of neocentromere properties relative to native centromere functions revealed that all twenty neocentromeres mediated chromosome segregation, albeit to different degrees. Some neocentromeres also caused reduced levels of transcription from genes found within the neocentromere region. Furthermore, like native centromeres, neocentromeres clustered in 3D with active/functional centromeres, indicating that formation of a new centromere mediates the reorganization of 3D nuclear architecture. This demonstrates that centromere clustering depends on epigenetically defined function and not on the primary DNA sequence, and that neocentromere function is independent of its distance from the native centromere position. Together, the results show that a neocentromere can form at many loci along a chromosome and can support the assembly of a functional kinetochore that exhibits native centromere functions including chromosome segregation accuracy and centromere clustering within the nucleus. [ABSTRACT FROM AUTHOR]

The increasing application of gene panels for familial cancer susceptibility disorders will probably lead to an increased proposal of susceptibility gene candidates. Using ERCC2 DNA repair gene as an example, we show that proof of a possible role in cancer susceptibility requires a detailed dissection and characterization of the underlying mutations for genes with diverse cellular functions (in this case mainly DNA repair and basic cellular transcription). In case of ERCC2, panel sequencing of 1345 index cases from 587 German, 405 Lithuanian and 353 Czech families with breast and ovarian cancer (BC/OC) predisposition revealed 25 mutations (3 frameshift, 2 splice-affecting, 20 missense), all absent or very rare in the ExAC database. While 16 mutations were unique, 9 mutations showed up repeatedly with population-specific appearance. Ten out of eleven mutations that were tested exemplarily in cell-based functional assays exert diminished excision repair efficiency and/or decreased transcriptional activation capability. In order to provide evidence for BC/OC predisposition, we performed familial segregation analyses and screened ethnically matching controls. However, unlike the recently published RECQL example, none of our recurrent ERCC2 mutations showed convincing co-segregation with BC/OC or significant overrepresentation in the BC/OC cohort. Interestingly, we detected that some deleterious founder mutations had an unexpectedly high frequency of > 1% in the corresponding populations, suggesting that either homozygous carriers are not clinically recognized or homozygosity for these mutations is embryonically lethal. In conclusion, we provide a useful resource on the mutational landscape of ERCC2 mutations in hereditary BC/OC patients and, as our key finding, we demonstrate the complexity of correct interpretation for the discovery of “bonafide” breast cancer susceptibility genes. [ABSTRACT FROM AUTHOR]

The hydrolytic deamination of adenosine to inosine (A-to-I editing) in precursor mRNA induces variable gene products at the post-transcription level. How and to what extent A-to-I RNA editing diversifies transcriptome is not fully characterized in the evolution, and very little is known about the selective constraints that drive the evolution of RNA editing events. Here we present a study on A-to-I RNA editing, by generating a global profile of A-to-I editing for a phylogeny of seven Drosophila species, a model system spanning an evolutionary timeframe of approximately 45 million years. Of totally 9281 editing events identified, 5150 (55.5%) are located in the coding sequences (CDS) of 2734 genes. Phylogenetic analysis places these genes into 1,526 homologous families, about 5% of total gene families in the fly lineages. Based on conservation of the editing sites, the editing events in CDS are categorized into three distinct types, representing events on singleton genes (type I), and events not conserved (type II) or conserved (type III) within multi-gene families. While both type I and II events are subject to purifying selection, notably type III events are positively selected, and highly enriched in the components and functions of the nervous system. The tissue profiles are documented for three editing types, and their critical roles are further implicated by their shifting patterns during holometabolous development and in post-mating response. In conclusion, three A-to-I RNA editing types are found to have distinct evolutionary dynamics. It appears that nervous system functions are mainly tested to determine if an A-to-I editing is beneficial for an organism. The coding plasticity enabled by A-to-I editing creates a new class of binary variations, which is a superior alternative to maintain heterozygosity of expressed genes in a diploid mating system. [ABSTRACT FROM AUTHOR]

The genome is often described as the information repository of an organism. Whether millions or billions of letters of DNA, its transmission across generations confers the principal medium for inheritance of organismal traits. Several emerging areas of research demonstrate that this definition is an oversimplification. Here, we explore ways in which a deeper understanding of genomic diversity and cell physiology is challenging the concepts of physical permanence attached to the genome as well as its role as the sole information source for an organism. [ABSTRACT FROM AUTHOR]

Cataracts are a significant public health problem with no proven methods for prevention. Discovery of novel disease mechanisms to delineate new therapeutic targets is of importance in cataract prevention and therapy. Herein, we report that mutations in the RagA GTPase (RRAGA), a key regulator of the mechanistic rapamycin complex 1 (mTORC1), are associated with autosomal dominant cataracts. We performed whole exome sequencing in a family with autosomal dominant juvenile-onset cataracts, and identified a novel p.Leu60Arg mutation in RRAGA that co-segregated with the disease, after filtering against the dbSNP database, and at least 123,000 control chromosomes from public and in-house exome databases. In a follow-up direct screening of RRAGA in another 22 families and 142 unrelated patients with congenital or juvenile-onset cataracts, RRAGA was found to be mutated in two unrelated patients (p.Leu60Arg and c.-16G>A respectively). Functional studies in human lens epithelial cells revealed that the RRAGA mutations exerted deleterious effects on mTORC1 signaling, including increased relocation of RRAGA to the lysosomes, up-regulated mTORC1 phosphorylation, down-regulated autophagy, altered cell growth or compromised promoter activity. These data indicate that the RRAGA mutations, associated with autosomal dominant cataracts, play a role in the disease by acting through disruption of mTORC1 signaling. [ABSTRACT FROM AUTHOR]

Drechmeria coniospora is an obligate fungal pathogen that infects nematodes via the adhesion of specialized spores to the host cuticle. D. coniospora is frequently found associated with Caenorhabditis elegans in environmental samples. It is used in the study of the nematode’s response to fungal infection. Full understanding of this bi-partite interaction requires knowledge of the pathogen’s genome, analysis of its gene expression program and a capacity for genetic engineering. The acquisition of all three is reported here. A phylogenetic analysis placed D. coniospora close to the truffle parasite Tolypocladium ophioglossoides, and Hirsutella minnesotensis, another nematophagous fungus. Ascomycete nematopathogenicity is polyphyletic; D. coniospora represents a branch that has not been molecularly characterized. A detailed in silico functional analysis, comparing D. coniospora to 11 fungal species, revealed genes and gene families potentially involved in virulence and showed it to be a highly specialized pathogen. A targeted comparison with nematophagous fungi highlighted D. coniospora-specific genes and a core set of genes associated with nematode parasitism. A comparative gene expression analysis of samples from fungal spores and mycelia, and infected C. elegans, gave a molecular view of the different stages of the D. coniospora lifecycle. Transformation of D. coniospora allowed targeted gene knock-out and the production of fungus that expresses fluorescent reporter genes. It also permitted the initial characterisation of a potential fungal counter-defensive strategy, involving interference with a host antimicrobial mechanism. This high-quality annotated genome for D. coniospora gives insights into the evolution and virulence of nematode-destroying fungi. Coupled with genetic transformation, it opens the way for molecular dissection of D. coniospora physiology, and will allow both sides of the interaction between D. coniospora and C. elegans, as well as the evolutionary arms race that exists between pathogen and host, to be studied. [ABSTRACT FROM AUTHOR]

Pemphigus vulgaris (PV) is a life-threatening autoimmune mucocutaneous blistering disease caused by disruption of intercellular adhesion due to auto-antibodies directed against epithelial components. Treatment is limited to immunosuppressive agents, which are associated with serious adverse effects. The propensity to develop the disease is in part genetically determined. We therefore reasoned that the delineation of PV genetic basis may point to novel therapeutic strategies. Using a genome-wide association approach, we recently found that genetic variants in the vicinity of the ST18 gene confer a significant risk for the disease. Here, using targeted deep sequencing, we identified a PV-associated variant residing within the ST18 promoter region (p<0.0002; odds ratio = 2.03). This variant was found to drive increased gene transcription in a p53/p63-dependent manner, which may explain the fact that ST18 is up-regulated in the skin of PV patients. We then discovered that when overexpressed, ST18 stimulates PV serum-induced secretion of key inflammatory molecules and contributes to PV serum-induced disruption of keratinocyte cell-cell adhesion, two processes previously implicated in the pathogenesis of PV. Thus, the present findings indicate that ST18 may play a direct role in PV and consequently represents a potential target for the treatment of this disease. [ABSTRACT FROM AUTHOR]

The exclusive localization of the histone H3 variant CENP-A to centromeres is essential for accurate chromosome segregation. Ubiquitin-mediated proteolysis helps to ensure that CENP-A does not mislocalize to euchromatin, which can lead to genomic instability. Consistent with this, overexpression of the budding yeast CENP-ACse4 is lethal in cells lacking Psh1, the E3 ubiquitin ligase that targets CENP-ACse4 for degradation. To identify additional mechanisms that prevent CENP-ACse4 misincorporation and lethality, we analyzed the genome-wide mislocalization pattern of overexpressed CENP-ACse4 in the presence and absence of Psh1 by chromatin immunoprecipitation followed by high throughput sequencing. We found that ectopic CENP-ACse4 is enriched at promoters that contain histone H2A.ZHtz1 nucleosomes, but that H2A.ZHtz1 is not required for CENP-ACse4 mislocalization. Instead, the INO80 complex, which removes H2A.ZHtz1 from nucleosomes, promotes the ectopic deposition of CENP-ACse4. Transcriptional profiling revealed gene expression changes in the psh1Δ cells overexpressing CENP-ACse4. The down-regulated genes are enriched for CENP-ACse4 mislocalization to promoters, while the up-regulated genes correlate with those that are also transcriptionally up-regulated in an htz1Δ strain. Together, these data show that regulating centromeric nucleosome localization is not only critical for maintaining centromere function, but also for ensuring accurate promoter function and transcriptional regulation. [ABSTRACT FROM AUTHOR]

Vascular endothelial growth factor (VEGF) is an angiogenic and neurotrophic factor, secreted by endothelial cells, known to impact various physiological and disease processes from cancer to cardiovascular disease and to be pharmacologically modifiable. We sought to identify novel loci associated with circulating VEGF levels through a genome-wide association meta-analysis combining data from European-ancestry individuals and using a dense variant map from 1000 genomes imputation panel. Six discovery cohorts including 13,312 samples were analyzed, followed by in-silico and de-novo replication studies including an additional 2,800 individuals. A total of 10 genome-wide significant variants were identified at 7 loci. Four were novel loci (5q14.3, 10q21.3, 16q24.2 and 18q22.3) and the leading variants at these loci were rs114694170 (MEF2C, P = 6.79x10-13), rs74506613 (JMJD1C, P = 1.17x10-19), rs4782371 (ZFPM1, P = 1.59x10-9) and rs2639990 (ZADH2, P = 1.72x10-8), respectively. We also identified two new independent variants (rs34528081, VEGFA, P = 1.52x10-18; rs7043199, VLDLR-AS1, P = 5.12x10-14) at the 3 previously identified loci and strengthened the evidence for the four previously identified SNPs (rs6921438, LOC100132354, P = 7.39x10-1467; rs1740073, C6orf223, P = 2.34x10-17; rs6993770, ZFPM2, P = 2.44x10-60; rs2375981, KCNV2, P = 1.48x10-100). These variants collectively explained up to 52% of the VEGF phenotypic variance. We explored biological links between genes in the associated loci using Ingenuity Pathway Analysis that emphasized their roles in embryonic development and function. Gene set enrichment analysis identified the ERK5 pathway as enriched in genes containing VEGF associated variants. eQTL analysis showed, in three of the identified regions, variants acting as both cis and trans eQTLs for multiple genes. Most of these genes, as well as some of those in the associated loci, were involved in platelet biogenesis and functionality, suggesting the importance of this process in regulation of VEGF levels. This work also provided new insights into the involvement of genes implicated in various angiogenesis related pathologies in determining circulating VEGF levels. The understanding of the molecular mechanisms by which the identified genes affect circulating VEGF levels could be important in the development of novel VEGF-related therapies for such diseases. [ABSTRACT FROM AUTHOR]