EvolMarkers: a database for mining exon and intron markers for evolution, ecology and conservation studies

Recent innovations in next-generation sequencing have lowered the cost of genome projects. Nevertheless, sequencing entire genomes for all representatives in a study remains expensive and unnecessary for most studies in ecology, evolution and conservation. It is still more cost-effective and efficient to target and sequence single-copy nuclear gene markers for such studies. Many tools have been developed for identifying nuclear markers, but most of these have focused on particular taxonomic groups. We have built a searchable database, EvolMarkers, for developing single-copy coding sequence (CDS) and exon-primed-intron-crossing (EPIC) markers that is designed to work across a broad range of phylogenetic divergences. The database is made up of single-copy CDS derived from BLAST searches of a variety of metazoan genomes. Users can search the database for different types of markers (CDS or EPIC) that are common to different sets of input species with different divergence characteristics. EvolMarkers can be applied to any taxonomic group for which genome data are available for two or more species. We included 82 genomes in the first version of EvolMarkers and have found the methods to be effective across Placozoa, Cnidaria, Arthropod, Nematoda, Annelida, Mollusca, Echinodermata, Hemichordata, Chordata and plants. We demonstrate the effectiveness of searching for CDS markers within annelids and show how to find potentially useful intronic markers within the lizard Anolis.     

Developing EPIC markers for chalcidoid Hymenoptera from EST and genomic data

Increasing numbers of phylogeographic studies make comparative inferences about the histories of co-distributed species. Although the aims of such studies are best achieved by jointly analysing sequences from multiple loci in a model-based framework, such data currently exist for few nonmodel systems. We used existing genomic data and expressed sequence tags (ESTs) for Hymenoptera and other insects to design intron-crossing primers for 40 loci, mainly ribosomal proteins, for chalcidoid parasitoids. Amplification success was scored on a range of taxa associated with two natural communities; oak galls and figs. Taxa were chosen at increasing distance from Nasonia, which was used for primer design, (i) within Pteromalids, (ii) within Chalcidoidea (Eupelmidae, Eulophidae, Eurytomidae, Ormyridae, Torymidae) and (iii) for a selection of distantly related gall and fig wasps (Cynipidae, Agaonidae). To assess the utility of these loci for phylogeographic and population genetic studies, we compared genetic diversity between Western Palaearctic refugia for two species. Our results show that it is feasible to design a large number of exon-primed-intron-crossing (EPIC) loci that may be informative about phylogeographic history within species but amplify across a large taxonomic range.     

SHLD2/FAM35A co‐operates with REV7 to coordinate DNA double‐strand break repair pathway choice

DNA double‐strand breaks (DSBs) can be repaired by two major pathways: non‐homologous end‐joining (NHEJ) and homologous recombination (HR). DNA repair pathway choice is governed by the opposing activities of 53BP1, in complex with its effectors RIF1 and REV7, and BRCA1. However, it remains unknown how the 53BP1/RIF1/REV7 complex stimulates NHEJ and restricts HR to the S/G2 phases of the cell cycle. Using a mass spectrometry (MS)‐based approach, we identify 11 high‐confidence REV7 interactors and elucidate the role of SHLD2 (previously annotated as FAM35A and RINN2) as an effector of REV7 in the NHEJ pathway. FAM35A depletion impairs NHEJ‐mediated DNA repair and compromises antibody diversification by class switch recombination (CSR) in B cells. FAM35A accumulates at DSBs in a 53BP1‐, RIF1‐, and REV7‐dependent manner and antagonizes HR by limiting DNA end resection. In fact, FAM35A is part of a larger complex composed of REV7 and SHLD1 (previously annotated as C20orf196 and RINN3), which promotes NHEJ and limits HR. Together, these results establish SHLD2 as a novel effector of REV7 in controlling the decision‐making process during DSB repair.     

Controlling the time evolution of mAb N‐linked glycosylation ‐ Part II: Model‐based predictions

N‐linked glycosylation is known to be a crucial factor for the therapeutic efficacy and safety of monoclonal antibodies (mAbs) and many other glycoproteins. The nontemplate process of glycosylation is influenced by external factors which have to be tightly controlled during the manufacturing process. In order to describe and predict mAb N‐linked glycosylation patterns in a CHO‐S cell fed‐batch process, an existing dynamic mathematical model has been refined and coupled to an unstructured metabolic model. High‐throughput cell culture experiments carried out in miniaturized bioreactors in combination with intracellular measurements of nucleotide sugars were used to tune the parameter configuration of the coupled models as a function of extracellular pH, manganese and galactose addition. The proposed modeling framework is able to predict the time evolution of N‐linked glycosylation patterns during a fed‐batch process as a function of time as well as the manipulated variables. A constant and varying mAb N‐linked glycosylation pattern throughout the culture were chosen to demonstrate the predictive capability of the modeling framework, which is able to quantify the interconnected influence of media components and cell culture conditions. Such a model‐based evaluation of feeding regimes using high‐throughput tools and mathematical models gives rise to a more rational way to control and design cell culture processes with defined glycosylation patterns. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:1135–1148, 2016     

Reorganization of inter‐chromosomal interactions in the 2q37‐deletion syndrome

Chromosomes occupy distinct interphase territories in the three‐dimensional nucleus. However, how these chromosome territories are arranged relative to one another is poorly understood. Here, we investigated the inter‐chromosomal interactions between chromosomes 2q, 12, and 17 in human mesenchymal stem cells (MSCs) and MSC‐derived cell types by DNA‐FISH. We compared our findings in normal karyotypes with a three‐generation family harboring a 2q37‐deletion syndrome, featuring a heterozygous partial deletion of histone deacetylase 4 (HDAC4) on chr2q37. In normal karyotypes, we detected stable, recurring arrangements and interactions between the three chromosomal territories with a tissue‐specific interaction bias at certain loci. These inter‐chromosomal interactions were confirmed by Hi‐C. Interestingly, the disease‐related HDAC4 deletion resulted in displaced inter‐chromosomal arrangements and altered interactions between the deletion‐affected chromosome 2 and chromosome 12 and/or 17 in 2q37‐deletion syndrome patients. Our findings provide evidence for a direct link between a structural chromosomal aberration and altered interphase architecture that results in a nuclear configuration, supporting a possible molecular pathogenesis.     

Selective base excision repair of DNA damage by the non‐base‐flipping DNA glycosylase AlkC

DNA glycosylases preserve genome integrity and define the specificity of the base excision repair pathway for discreet, detrimental modifications, and thus, the mechanisms by which glycosylases locate DNA damage are of particular interest. Bacterial AlkC and AlkD are specific for cationic alkylated nucleobases and have a distinctive HEAT‐like repeat (HLR) fold. AlkD uses a unique non‐base‐flipping mechanism that enables excision of bulky lesions more commonly associated with nucleotide excision repair. In contrast, AlkC has a much narrower specificity for small lesions, principally N3‐methyladenine (3mA). Here, we describe how AlkC selects for and excises 3mA using a non‐base‐flipping strategy distinct from that of AlkD. A crystal structure resembling a catalytic intermediate complex shows how AlkC uses unique HLR and immunoglobulin‐like domains to induce a sharp kink in the DNA, exposing the damaged nucleobase to active site residues that project into the DNA. This active site can accommodate and excise N3‐methylcytosine (3mC) and N1‐methyladenine (1mA), which are also repaired by AlkB‐catalyzed oxidative demethylation, providing a potential alternative mechanism for repair of these lesions in bacteria.     

Distinct roles of ATR and DNA‐PKcs in triggering DNA damage responses in ATM‐deficient cells

The cellular response to DNA double‐strand breaks involves direct activation of ataxia telangiectasia mutated (ATM) and indirect activation of ataxia telangiectasia and Rad3 related (ATR) in an ATM/Mre11/cell‐cycle‐dependent manner. Here, we report that the crucial checkpoint signalling proteins—p53, structural maintainance of chromosomes 1 (SMC1), p53 binding protein 1 (53BP1), checkpoint kinase (Chk)1 and Chk2—are phosphorylated rapidly by ATR in an ATM/Mre11/cell‐cycle‐independent manner, albeit at low levels. We observed the sequential recruitment of replication protein A (RPA) and ATR to the sites of DNA damage in ATM‐deficient cells, which provides a mechanistic basis for the observed phosphorylations. The recruitment of ATR and consequent phosphorylations do not require Mre11 but are dependent on Exo1. We show that these low levels of phosphorylation are biologically important, as ATM‐deficient cells enforce an early G2/M checkpoint that is ATR‐dependent. ATR is also essential for the late G2 accumulation that is peculiar to irradiated ATM‐deficient cells. Interestingly, phosphorylation of KRAB associated protein 1 (KAP‐1), a protein involved in chromatin remodelling, is mediated by DNA‐dependent protein kinase catalytic subunit (DNA‐PKcs) in a spatio‐temporal manner in addition to ATM. We posit that ATM substrates involved in cell‐cycle checkpoint signalling can be minimally phosphorylated independently by ATR, while a small subset of proteins involved in chromatin remodelling are phosphorylated by DNA‐PKcs in addition to ATM.     

An inexpensive,automation‐friendly protocol for recovering high‐quality DNA

Although commercial kits are available for automated DNA extraction, ‘artisanal’ protocols are not. In this study, we present a silica‐based method that is sensitive, inexpensive and compliant with automation. The effectiveness of this protocol has now been tested on more than 5000 animal specimens with highly positive results.     

Method for preparing single-stranded DNA templates for Pyrosequencing using vector ligation and universal biotinylated primers

In Pyrosequencing, the addition of nucleotides to a primer-template hybrid is monitored by enzymatic conversion of chemical energy into detectable light. The technique yields both qualitative and quantitative sequence information because the chemical energy is released by a stoichiometric split off of pyrophosphates from incorporated deoxynucleotide triphosphates and a defined nucleotide dispensation order is given. Because Pyrosequencing works best if single-stranded DNA templates are used, template generation usually requires PCR with a target-specific biotinylated primer and a subsequent purification involving interaction of the biotin label with immobilized streptavidin. To circumvent the need for numerous and expensive template-specific biotinylated primers, we developed a method that uses the ligation of amplified DNA fragments into a plasmid vector, thereby facilitating subsequent PCR using a universal vector-specific biotinylated primer. This approach allows easy and straightforward isolation of single-stranded templates of any PCR product. As a proof of principle, we used the method for genotyping two single-nucleotide polymorphisms in the human genes CARD15 and A2M and for characterization of four multisite variations in the human DEFB104 gene.     

Sequence and evolution of rhesus monkey alphoid DNA

Summary Analysis of rhesus monkey alphoid DNA suggests that it arose by tandem duplication of an ancestral monomer unit followed by independent variation within two adjacent monomers (one becoming more divergent than the other) before their amplification as a dimer unit to produce tandem arrays. The rhesus monkey alphoid DNA is a tandemly repeated, 343-bp dimer; the consensus dimer is over 98% homologous to the alphoid dimers reported for baboon and bonnet monkey, 81% homologous to the African green monkey alpha monomer, and less than 70% homologous to the more divergent human alphoid DNAs. The consensus dimer consists of two wings (I and II, 172 and 171 bp, respectively) that are only 70% homologous to each other, but share seven regions of exact homology. These same regions are highly conserved among the consensus sequences of the other cercopithecid alphoid DNAs. The three alpha-protein binding sites reported for African green monkey alpha DNA by F. Strauss and A. Varshavsky (Cell 37: 889–901, 1984) occur in wings I and II, but with one site altered in wing I. Two cloned dimer segments are 98% homologous to the consensus, each containing 8 single-base-pair differences within the 343-bp segment. Surprisingly, 37% of these differences occur in regions that are evolutionarily conserved in the alphoid consensus sequences, including the alpha-protein binding sites. Sequence variation in this highly repetitive DNA family may produce unique nucleosomal architectures for different members of an alphoid array. These unique architectures may modulate the evolution of these repetitive DNAs and may produce unique centromeric characteristics in primate chromosomes.     

A genome‐wide association study using selective DNA pooling identifies candidate markers for fertility in Holstein cattle

The decline in the reproductive efficiency of dairy cows, especially those with high producing potential, has become a challenging problem. In this study, a selective DNA pooling approach was applied to a cow population whose oocytes were fertilized and cultured to obtain phenotypic records of fertilization rate and blastocyst rate. Using a stringent 5% genome‐wide significance level, 22 and five single nucleotide polymorphisms (SNPs) were found to be associated with fertilization rate and blastocyst rate, respectively. SNPs that showed significant association in selective DNA pooling were further evaluated by individual genotyping. Interestingly, the majority of the SNP associations were confirmed by individual genotyping, testifying to the effectiveness of selective DNA pooling using a high‐density SNP genotyping array. This study is the first application of the selective DNA pooling approach using the BovineSNP50 array in cattle.     

Structural insight into dynamic bypass of the major cisplatin‐DNA adduct by Y‐family polymerase Dpo4

Y‐family DNA polymerases bypass Pt‐GG, the cisplatin‐DNA double‐base lesion, contributing to the cisplatin resistance in tumour cells. To reveal the mechanism, we determined three structures of the Y‐family DNA polymerase, Dpo4, in complex with Pt‐GG DNA. The crystallographic snapshots show three stages of lesion bypass: the nucleotide insertions opposite the 3′G (first insertion) and 5′G (second insertion) of Pt‐GG, and the primer extension beyond the lesion site. We observed a dynamic process, in which the lesion was converted from an open and angular conformation at the first insertion to a depressed and nearly parallel conformation at the subsequent reaction stages to fit into the active site of Dpo4. The DNA translocation‐coupled conformational change may account for additional inhibition on the second insertion reaction. The structures illustrate that Pt‐GG disturbs the replicating base pair in the active site, which reduces the catalytic efficiency and fidelity. The in vivo relevance of Dpo4‐mediated Pt‐GG bypass was addressed by a dpo‐4 knockout strain of Sulfolobus solfataricus, which exhibits enhanced sensitivity to cisplatin and proteomic alterations consistent with genomic stress.