Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Restriction site heteroplasmy in the mitochondrial DNA of Brycon opalinus (Cuvier, 1819) (Characiformes, Characidae, Bryconiae)

Homoplasmy is a feature usually found in the mtDNA of higher animal taxa. On the other hand, the presence of two classes of mtDNA in the same cell or organism is rare and may appear in length or site variation. Data from mtDNA RFLP analysis of Brycon opalinus populations (Cuvier, 1819; Characiformes, Characidae, Bryconinae) revealed site heteroplasmy from endonuclease NheI digestion. Southern blotting hybridization was used to survey a total of 257 specimens with 24 restriction enzymes. Three different restriction fragment patterns of mtDNA were obtained from NheI digestion. Two individuals from hatchery broodstock were found to have two of them. NheI digests of heteroplasmic individuals yielded two fragments of approximately 1180 and 1260 bp. Despite the low frequency of this type of heteroplasmy in the whole B. opalinus population, the presence of site heteroplasmy in this species supports the evidence of this phenomenon in lower vertebrate groups.

Caracterização dos padrões genéticos de populações invasoras e naturalizadas de schizolobium parahyba (Caesalpinioideae – Fabaceae) por restriction-site associated dna-sequencing

Pós-graduação em Ciência Florestal - FCA

Genetic Structure of the Florida Key Tree Cactus, Pilosocereus robinii, using Restriction Site associated DNA (RAD) markers

Rare plant conservation efforts must utilize current genetic methods to ensure the evolutionary potential of populations is preserved. One such effort involves the Key Tree Cactus, Pilosocereus robinii, which is an endangered columnar cactus native to the Florida Keys. The populations have precipitously declined over the past decade because of habitat loss and increasing soil salinity from rising sea levels and storm surge. Next-generation DNA sequencing was used to assess the genetic structure of the populations. Twenty individuals representative of both wild and extirpated cacti were chosen for Restriction Site Associated DNA (RAD) analysis. Samples processed using the HindIII and NotIII restriction enzymes produced 82,382,440 high quality reads used for genetic mapping, from which 5,265 Single Nucleotide Polymorphisms (SNPs) were discovered. The analysis revealed that the Keys’ populations are closely related with little population differentiation. In addition, the populations display evidence of inbreeding and low genetic diversity.

Rapid Site-Directed Mutagenesis Using Two-PCR-Generated DNA Fragments Reproducing the Plasmid Template.

We describe a new rapid and efficient polymerase chain reaction (PCR)-based site-directed mutagenesis method. This procedure is effective with any plasmid and it employs four oligonucleotide primers. One primer contains the desired mutation, the second is oriented in the opposite direction (one of these two primers should be phosphorylated), and the third and fourth should be coding in complementary fashion for a unique restriction site to be introduced in a nonessential region. The method consists of two simultaneous PCR reactions; the PCR products are digested with the enzyme that recognizes the newly introduced unique restriction site and then ligased and used to transform competent bacteria. Additionally, the use of Dpn I facilitates the elimination of template DNA. The newly introduced restriction site is essential for ligation in the correct orientation of the two-PCR products and is further used for mutant screening. Resulting plasmids carry both the new restriction site and the desired mutation. Using this method, more than 20 mutants have already been generated (using two different kinds of templates); all these mutants were sequenced for the desired mutation and transfected into AtT-20 cells and the expressed mutant proteins encoded by the vector were assayed.

A FUNCTIONAL VARIANT OF HLA-DR1 DELINEATED BY T-LYMPHOCYTE CLONE REACTIVITY, PROTEIN CONFORMATION, AND GENOMIC RESTRICTION SITES

This research characterized a serologically indistinguishable form of HLA-DR1 that: (1) cannot stimulate some DR1-restricted or specific T-lymphocyte clones; (2) displays an unusual electrophoretic pattern on two dimensional gels; and (3) is marked by a polymorphic restriction site of the alpha gene. Inefficient stimulation of some DR1-restricted clones was a property of DR1$\sp{+}$ cells that shared HLA-B14 on the same haplotype and/or were carriers of 21-hydroxylase (21-OH) deficiency. Nonclassical 21-OH deficiency frequently demonstrates genetic linkage with HLA-B14;DR1 haplotypes and associates with duplications of C4B and one 21-OH gene. Cells having both stimulatory (DR1$\sb{\rm n}$) and nonstimulatory (DR1$\sb{\rm x}$) parental haplotypes did not mediate proliferation of these clones. However, heterozygous DR1$\sb{\rm x}$, 2 and DR1$\sb{\rm x}$, 7 cells were efficient stimulators of DR2 and DR7 specific clones, respectively, suggesting that a trans acting factor may modify DR1 alleles or products to yield a dominant DR1$\sb{\rm x}$ phenotype. Incompetent stimulator populations did not secrete an intercellular soluble or contact dependent suppressor factor nor did they express interleukin-2 receptors competing for T-cell growth factors. Two dimensional gel analysis of anti-DR immunoprecipitates revealed, in addition to normal DR$\alpha$ and DR$\beta$ chains, a 50kD species from DR1$\sb{\rm x}$ but not from the majority of DR1$\sb{\rm n}$ or non-DR1 cells. The 50kD structure was stable under reducing conditions in SDS and urea, had antigenic homology with DR, and dissociated after boiling into 34kD and 28kD peptide chains apparently identical with DR$\alpha$ and DR$\beta$ as shown by limited digest peptide maps. N-linked glycosylation and sialation of DRgp50 appeared to be unchanged from normal DR$\alpha$ and DR$\beta$. Bg1II digestion and $DR\alpha$ probing of DR1$\sb{\rm x}$ genomic DNA revealed a 4.5kb fragment while DR1$\sb{\rm n}$ DNA yielded 3.8 and 0.76kb fragments; all restriction sites mapped to the 3$\sp\prime$ untranslated region of $DR\alpha$. Collectively, these data suggest that DRgp50 represents a novel combinatorial association between constitutive chains of DR that may interfere with or compete for normal T cell receptor recognition of DR1 as both an alloantigen and restricting element. Furthermore, extensive chromosomal abnormalities previously mapped to the class III region of B14;DR1 haplotypes may extend into the adjacent class II region with consequent intrusion on immune function. ^

Sequence-related structural DNA anomalies affect restriction enzyme activity and DNA polymerase I binding

DNA serves as a target molecule for several types of enzymes and may assume a wide variety of structural motifs depending upon the local sequence. The BssHII restriction site (GC)3 resides in a 9bp region of alternating pyrimidine and purine residues within the &phis;X174 genome. Such sequences are known to demonstrate non-canonical helical behavior under the appropriate conditions. The kinetics of BssHII cleavage was investigated in supercoiled and linear plasmid DNA, and in a 323bp DNA fragment obtained via amplification of &phis;X174. The rate of enzyme cleavage was enhanced in the supercoiled form and in the presence of 50μM cobalt hexamine. Similarly, cobalt hexamine was also found to enhance TaqI activity directly adjacent to the (GC)3 region. ^ Initial DNA polymerase I binding studies (including a gel mobility shift assay and a protection assay) indicated a notable interaction between DNA polymerase I and the BssHII site. An in-depth study revealed that equilibrium binding of DNA polymerase I to the T7 RNA polymerase promoter was comparable to that of the (GC)3 site, however the strongest interaction was observed with a cruciform containing region. Increasing the ionic strength of the solution environment, including the addition of DNA polymerase I reaction buffer significantly decreased the equilibrium dissociation constant values. ^ It is suggested that the region within or around the BssHII site experiences a conformational change generating a novel structure under the influence of supercoiled tension or 50μM cobalt hexamine. It is proposed that this transition may enhance enzyme activity and binding by providing an initial enzyme-docking site—the rate-limiting step in restriction enzyme kinetics. The high binding potential of DNA polymerase I for each of the motifs described, is hypothesized to be due to recognition of the structural DNA anomalies by the 3′–5′ exonuclease domain. ^

Molecular typing of dengue virus type 2 in Brazil

Strain typing is a critical tool for molecular epidemiological analysis and can provide important information about the spread of dengue viruses. Here, we performed a molecular characterization of DEN-2 viruses isolated in Brazil during 1990-2000 from geographically and temporally distinct areas in order to investigate the genetic distribution of this serotype circulating in the country. Restriction site-specific polymerase chain reaction (RSS)-PCR presented the same pattern for all 52 Brazilian samples, showing the circulation of just one DEN-2 variant. Phylogenetic analysis using progressive pairwise alignments from 240-nucleotide sequences of the E/NS1 junction in 15 isolates showed that they belong to genotype III (Jamaica genotype).

Genomics of the divergence continuum in an African plant biodiversity hotspot, I: drivers of population divergence in Restio capensis (Restionaceae).

Understanding the drivers of population divergence, speciation and species persistence is of great interest to molecular ecology, especially for species-rich radiations inhabiting the world's biodiversity hotspots. The toolbox of population genomics holds great promise for addressing these key issues, especially if genomic data are analysed within a spatially and ecologically explicit context. We have studied the earliest stages of the divergence continuum in the Restionaceae, a species-rich and ecologically important plant family of the Cape Floristic Region (CFR) of South Africa, using the widespread CFR endemic Restio capensis (L.) H.P. Linder & C.R. Hardy as an example. We studied diverging populations of this morphotaxon for plastid DNA sequences and >14 400 nuclear DNA polymorphisms from Restriction site Associated DNA (RAD) sequencing and analysed the results jointly with spatial, climatic and phytogeographic data, using a Bayesian generalized linear mixed modelling (GLMM) approach. The results indicate that population divergence across the extreme environmental mosaic of the CFR is mostly driven by isolation by environment (IBE) rather than isolation by distance (IBD) for both neutral and non-neutral markers, consistent with genome hitchhiking or coupling effects during early stages of divergence. Mixed modelling of plastid DNA and single divergent outlier loci from a Bayesian genome scan confirmed the predominant role of climate and pointed to additional drivers of divergence, such as drift and ecological agents of selection captured by phytogeographic zones. Our study demonstrates the usefulness of population genomics for disentangling the effects of IBD and IBE along the divergence continuum often found in species radiations across heterogeneous ecological landscapes.

'Next generation' biogeography: towards understanding the drivers of species diversification and persistence

The drivers of species diversification and persistence are of great interest to current biogeography, especially in those global biodiversity hotspots' harbouring most of Earth's animal and plant life. Classical multispecies biogeographical work has yielded fascinating insights into broad-scale patterns of diversification, and DNA-based intraspecific phylogeographical studies have started to complement this picture at much finer temporal and spatial scales. The advent of novel next-generation sequencing (NGS) technologies provides the opportunity to greatly scale up the numbers of individuals, populations and species sampled, potentially merging intraspecific and interspecific approaches to biogeographical inference. Here, we outline these prospects and issues by using the example of an undisputed hotspot, the Cape of southern Africa. We outline the current state of knowledge on the biogeography of species diversification within the Cape, review the literature for phylogeographical evidence of its likely drivers and mechanisms, and suggest possible ways forward based on NGS approaches. We demonstrate the potential of these methods and current bioinformatic issues with the help of restriction-site-associated DNA (RAD) sequencing data for three highly divergent species of the Restionaceae, an important plant radiation in the Cape. A thorough understanding of the mechanisms that facilitate species diversification and persistence in spatially structured, species-rich environments will require the adoption of novel genomic and bioinformatic tools in biogeographical studies.

Temporal distribution of dengue virus serotypes in Colombian endemic area and dengue incidence: re-introduction of dengue-3 associated to mild febrile illness and primary infection

We have investigated the temporal distribution of dengue (DEN) virus serotypes in the department (state) of Santander, Colombia, in relation to dengue incidence, infection pattern, and severity of disease. Viral isolation was attended on a total of 1452 acute serum samples collected each week from 1998 to 2004. The infection pattern was evaluated in 596 laboratory-positive dengue cases using an IgG ELISA, and PRNT test. The dengue incidence was documented by the local health authority. Predominance of DEN-1 in 1998 and DEN-3 re-introduction and predominance in 2001-2003 coincided with outbreaks. Predominance of DEN-2 in 2000-2001 coincided with more dengue hemorrhagic fever (DHF). DEN-4 was isolated in 2000-2001 and 2004 but was not predominant. There was an annual increase of primary dengue infections (from 13.7 to 81.4%) that correlated with frequency of DEN-3 (r = 0.83; P = 0.038). From the total number of primary dengue infections DEN-3 (81.3%) was the most frequent serotype. DHF was more frequent in DEN-2 infected patients than in DEN-3 infected patients: 27.5 vs 10.9% (P < 0.05). DEN-3 viruses belonged to subtype C (restriction site-specific-polymerase chain reaction) like viruses isolated in Sri-Lanka and other countries in the Americas. Our findings show the importance of continuous virological surveillance to identify the risk factors of dengue epidemics and severity.

Association of bGH and Pit-1 gene variants with milk production traits in dairy Gyr bulls

The objective of this study was to obtain genetic marker information in the Gyr breed by analyzing bGH and Pit-1 gene polymorphisms and to verify their association with milk production traits. One sample including 40 Gyr bulls were genotyped at two bGH gene restriction sites (bGH- AluI and bGH-MspI) and at one restriction site in the Pit-1 gene (Pit-1 HinfI). The bGH-MspI(-) allele was favorable for fat milk percentage. The heterozigous Pit-1 HinfI (+/-) bulls were superior for fat milk production, in relation to homozigous Pit-1 HinfI (+/+). The Pit-1 and bGH genes are strong candidates in the dairy cattle QTL search, and zebuine populations are promising samples for this purpose.

Hybridization Capture Using RAD Probes (hyRAD), a New Tool for Performing Genomic Analyses on Collection Specimens.

In the recent years, many protocols aimed at reproducibly sequencing reduced-genome subsets in non-model organisms have been published. Among them, RAD-sequencing is one of the most widely used. It relies on digesting DNA with specific restriction enzymes and performing size selection on the resulting fragments. Despite its acknowledged utility, this method is of limited use with degraded DNA samples, such as those isolated from museum specimens, as these samples are less likely to harbor fragments long enough to comprise two restriction sites making possible ligation of the adapter sequences (in the case of double-digest RAD) or performing size selection of the resulting fragments (in the case of single-digest RAD). Here, we address these limitations by presenting a novel method called hybridization RAD (hyRAD). In this approach, biotinylated RAD fragments, covering a random fraction of the genome, are used as baits for capturing homologous fragments from genomic shotgun sequencing libraries. This simple and cost-effective approach allows sequencing of orthologous loci even from highly degraded DNA samples, opening new avenues of research in the field of museum genomics. Not relying on the restriction site presence, it improves among-sample loci coverage. In a trial study, hyRAD allowed us to obtain a large set of orthologous loci from fresh and museum samples from a non-model butterfly species, with a high proportion of single nucleotide polymorphisms present in all eight analyzed specimens, including 58-year-old museum samples. The utility of the method was further validated using 49 museum and fresh samples of a Palearctic grasshopper species for which the spatial genetic structure was previously assessed using mtDNA amplicons. The application of the method is eventually discussed in a wider context. As it does not rely on the restriction site presence, it is therefore not sensitive to among-sample loci polymorphisms in the restriction sites that usually causes loci dropout. This should enable the application of hyRAD to analyses at broader evolutionary scales.

Sequence diversity in the coat protein gene of Lettuce big-vein associated virus and Mirafiori lettuce big-vein virus infecting lettuce in Brazil

Lettuce big vein associated virus (LBVaV) and Mirafiori lettuce big vein virus (MLBVV) have been found in mixed infection in Brazil causing the lettuce big vein disease. Analysis of part of the coat protein (CP) gene of Brazilian isolates of LBVaV collected from lettuce, showed at least 93% amino acid sequence identity with other LBVaV isolates. Genetic diversity among MLBVV CP sequences was higher when compared to LBVaV CP sequences, with amino acid sequence identity ranging between 91% to 100%. Brazilian isolates of MLBVV belong to subgroup A, with one RsaI restriction site on the coat protein gene. There is no indication for a possible geografical origin for the Brazilian isolates of LBVaV and MLBVV.