Rafael Montiel

Polar and brown bear genomes reveal ancient admixture and demographic footprints of past climate change

Polar bears (PBs) are superbly adapted to the extreme Arctic environment and have become emblematic of the threat to biodiversity from global climate change. Their divergence from the lower-latitude brown bear provides a textbook example of rapid evolution of distinct phenotypes. However, limited mitochondrial and nuclear DNA evidence conflicts in the timing of PB origin as well as placement of the species within versus sister to the brown bear lineage. We gathered extensive genomic sequence data from contemporary polar, brown, and American black bear samples, in addition to a 130,000- to 110,000-y old PB, to examine this problem from a genome-wide perspective. Nuclear DNA markers reflect a species tree consistent with expectation, showing polar and brown bears to be sister species. However, for the enigmatic brown bears native to Alaskas Alexander Archipelago, we estimate that not only their mitochondrial genome, but also 5–10% of their nuclear genome, is most closely related to PBs, indicating ancient admixture between the two species. Explicit admixture analyses are consistent with ancient splits among PBs, brown bears and black bears that were later followed by occasional admixture. We also provide paleodemographic estimates that suggest bear evolution has tracked key climate events, and that PB in particular experienced a prolonged and dramatic decline in its effective population size during the last ca. 500,000 years. We demonstrate that brown bears and PBs have had sufficiently independent evolutionary histories over the last 4–5 million years to leave imprints in the PB nuclear genome that likely are associated with ecological adaptation to the Arctic environment.

Application and comparison of large-scale solution-based DNA capture-enrichment methods on ancient DNA

The development of second-generation sequencing technologies has greatly benefitted the field of ancient DNA (aDNA). Its application can be further exploited by the use of targeted capture-enrichment methods to overcome restrictions posed by low endogenous and contaminating DNA in ancient samples. We tested the performance of Agilents SureSelect and Mycroarrays MySelect in-solution capture systems on Illumina sequencing libraries built from ancient maize to identify key factors influencing aDNA capture experiments. High levels of clonality as well as the presence of multiple-copy sequences in the capture targets led to biases in the data regardless of the capture method. Neither method consistently outperformed the other in terms of average target enrichment, and no obvious difference was observed either when two tiling designs were compared. In addition to demonstrating the plausibility of capturing aDNA from ancient plant material, our results also enable us to provide useful recommendations for those planning targeted-sequencing on aDNA.

Genetic Structure and Origin of Peopling in The Azores Islands (Portugal): The View from mtDNA

The Azores islands (Portugal), uninhabited when discovered by Portuguese navigators in the fifteenth century, are located in the Atlantic Ocean 1500 km from the European mainland. The archipelago is formed by nine islands of volcanic origin that define three geographical groups: Eastern (S. Miguel and Sta. Maria), Central (Terceira, Faial, Pico, Graciosa and S. Jorge) and Western (Flores and Corvo). To improve the genetic characterisation of the Azorean population, and to clarify some aspects related to the history of settlement, a study of mtDNA was conducted in the population of the archipelago. The HVRI region was sequenced and specific RFLPs were screened in 146 samples obtained from unrelated individuals with Azorean ancestry (50 from the Eastern group, 60 from the Central group, and 37 from the Western group). Samples were classified into haplogroups based on the information obtained from both sequencing and RFLP analysis.

An apoptosis-inducing serine protease secreted by the entomopathogenic nematode Steinernema carpocapsae.

Steinernema carpocapsae is an insect parasitic nematode able to parasitise and kill the host within 48 h. Secreted products (ESP) of the parasitic stage of a virulent strain contain higher amounts of proteolytic activity than a low virulence strain, suggesting proteases are involved in virulence. From the ESP we purified a protein (Sc-SP-3) with a M(r) of 30 kDa and a pI of 7 that cleaved the synthetic substrate N-succinyl-Ala-Ala-Pro-Phe-pNA and was inhibited by phenylmethanesulfonyl fluoride, benzamidine and chymostatin, thus indicating that it belongs to the chymotrypsin-like serine protease family. Sc-SP-3 has a V(max) of 0.3 mM min(-1)ml(-1) and K(m) of 6.6 x 10(-4)M, with maximum activity at pH 8 and 40 degrees C. The full-length cDNA was obtained using degenerate oligonucleotides for serine proteases. This open reading frame encodes a preproprotein containing a putative signal peptide composed of 16 amino acid residues, a prodomain of 40 residues and a mature protease domain of 261 residues, including the catalytic triad His/Asp/Ser characteristic of trypsin-like serine proteases. The N-terminal sequence and the peptide masses fingerprint obtained by MALDI-TOF-MS for the purified protein matched the cDNA. Gene expression analysis by quantitative real-time-PCR showed that this gene is expressed only during the parasitic stage and that pre-invasive nematodes inside the mid-gut expressed higher amounts of Sc-SP-3 than those that already enter the haemocoel. Sc-SP-3 caused histolysis in the insect mid-gut. In vitro assays demonstrated that Sc-SP-3 digested extracellular proteins and induced apoptosis in Sf9 insect cells, thus suggesting Sc-SP-3 is a multifunctional chymotrypsin-like protease involved in pathogenesis.

Authenticating ancient human mitochondrial DNA.

The use of ancient DNA techniques in human studies has been hampered by problems of contamination with modern human DNA. The main problem has been that the object of study belongs to the same species as the observer, and the complete elimination of the contamination risk is seemingly unlikely. Contamination has even been detected in the most specialized laboratories in this field. In these kinds of studies it is therefore very important to detect contamination and to distinguish contaminants from authentic results. Here, we report the use of a strategy to authenticate the identity of ancient mitochondrial DNA (mtDNA), based on the previously established relationship between D-loop sequence substitutions and haplogroup-specific restriction site changes. Forty-four individuals from a 16th-century necropolis were analyzed, from which 28 control region sequences were obtained. These sequences were preclassified into haplogroups, according to the observed motifs. Subsequently, the DNA extracts from which the sequences were obtained, along with independent extracts of subsets of the same individuals, were subjected to restriction fragment length polymorphism (RFLP) analysis to compare and corroborate the results. Using this approach, 24 sequences were authenticated, while two were discarded because of result mismatches. The final distribution of the haplogroups in the sample, and the differences in the sequences, are two additional criteria of authentication.

Serine Protease-mediated Host Invasion by the Parasitic Nematode Steinernema carpocapsae

Steinernema carpocapsae is an insect parasitic nematode used in biological control, which infects insects penetrating by mouth and anus and invading the hemocoelium through the midgut wall. Invasion has been described as a key factor in nematode virulence and suggested to be mediated by proteases. A serine protease cDNA from the parasitic stage was sequenced (sc-sp-1); the recombinant protein was produced in an Escherichia coli system, and a native protein was purified from the secreted products. Both proteins were confirmed by mass spectrometry to be encoded by the sc-sp-1 gene. Sc-SP-1 has a pI of 8.7, a molecular mass of 27.3 kDa, a catalytic efficiency of 22.2 × 104 s−1 m−1 against N-succinyl-Ala-Ala-Pro-Phe-pNA, and is inhibited by chymostatin (IC 0.07) and PMSF (IC 0.73). Sc-SP-1 belongs to the chymotrypsin family, based on sequence and biochemical analysis. Only the nematode parasitic stage expressed sc-sp-1. These nematodes in the midgut lumen, prepared to invade the insect hemocoelium, expressed higher levels than those already in the hemocoelium. Moreover, parasitic nematode sense insect peritrophic membrane and hemolymph more quickly than they do other tissues, which initiates sc-sp-1 expression. Ex vivo, Sc-SP-1 was able to bind to insect midgut epithelium and to cause cell detachment from basal lamina. In vitro, Sc-SP-1 formed holes in an artificial membrane model (Matrigel), whereas Sc-SP-1 treated with PMSF did not, very likely because it hydrolyzes matrix glycoproteins. These findings highlight the S. carpocapsae-invasive process that is a key step in the parasitism thus opening new perspectives for improving nematode virulence to use in biological control.

Determination of Human Caucasian Mitochondrial DNA Haplogroups by Means of a Hierarchical Approach

In this paper we propose a hierarchical approach that allows the screening of mitochondrial DNA (mtDNA) haplogroups in populations that have essentially West Eurasian mtDNA backgrounds but that could have some non-West Eurasian contributions. To develop and validate this scheme, we used data on 18 coding region polymorphisms (17 analyzed by RFLP analysis and 1 by sequencing) and sequences of hypervariable segment I (HVSI) of the mtDNA control region from the Azores Islands (Portugal) population. The proposed scheme allows the characterization of almost all West Eurasian and African major clusters by means of RFLPs. Furthermore, the scheme includes information on situations in which sequencing is pertinent to defining a particular haplogroup. The validity of the scheme is ensured by (1) using relatively stable polymorphic positions, (2) screening more than one position to define a specific haplogroup, and (3) typing confirmatory positions. Dubious samples can be resolved by sequencing. The robustness of this approach was assessed by sequencing all samples for HVSI, taking advantage of the previously established relationships between RFLPs and control region sequence polymorphisms. The use of this hierarchical approach avoids the screening of unnecessary control region polymorphisms and therefore results in a more rapid and cost-efficient screening than one in which all polymorphic positions are analyzed. Even if this approach leads to a lower level of phylogeographic resolution than the sequencing of all samples, it allows us to define population movements on a continental level and can be applied, unlike sequencing all samples, with a low cost in any laboratory. [End Page 431]

Increased transcript diversity: novel splicing variants of Machado–Joseph Disease gene (ATXN3)

Machado–Joseph disease (MJD) is a late-onset neurodegenerative disorder that presents clinical heterogeneity not completely explained by its causative mutation. MJD is caused by an expansion of a CAG tract at exon 10 of the ATXN3 gene (14q32.1), which encodes for ataxin-3. The main goal of this study was to analyze the occurrence of alternative splicing at the ATXN3 gene, by sequencing a total of 415 cDNAs clones (from 20 MJD patients and 14 controls). Two novel exons are described for the ATXN3 gene. Fifty-six alternative splicing variants, generated by four types of splicing events, were observed. From those variants, 50 were not previously described, and 26 were only found in MJD patients samples. Most of the variants (85.7%) present frameshift, which leads to the appearance of premature stop codons. Thirty-seven of the observed variants constitute good targets to nonsense-mediated decay, the remaining are likely to be translated into at least 20 different isoforms. The presence of ataxin-3 domains was assessed, and consequences of domain disruption are discussed. The present study demonstrates high variability in the ATXN3 gene transcripts, providing a basis for further investigation on the contribution of alternative splicing to the MJD pathogenic process, as well as to the larger group of the polyglutamine disorders.

Transcripts analysis of the entomopathogenic nematode Steinernema carpocapsae induced in vitro with insect haemolymph

Steinernema carpocapsae is an insect parasitic nematode widely used in pest control programs. The efficacy of this nematode in controlling insects has been found to be related to the pathogenicity of the infective stage. In order to study the parasitic mechanisms exhibited by this parasite, a cDNA library of the induced S. carpocapsae parasitic phase was generated. A total of 2500 clones were sequenced and 2180 high-quality ESTs were obtained from this library. Cluster analysis generated a total of 1592 unique sequences including 1393 singletons. About 63% of the unique sequences had significant hits (e</=1e-05) to the non-redundant protein database. The remaining sequences most likely represent putative novel protein coding genes. Comparative analysis identified 377 homologs in C. elegans, 431 in C. briggsae and 75 in other nematodes. Classification of the predicted proteins revealed involvement in diverse cellular, metabolic and extracellular functions. One hundred and nineteen clusters were predicted to encode putative secreted proteins such as proteases, proteases inhibitors, lectins, saposin-like proteins, acetylcholinesterase, anti-oxidants, and heat-shock proteins, which can possibly have host interactions. This dataset provides a basis for genomic studies towards a better understanding of the events that occur in the parasitic process of this entomopathogenic nematode, including invasion of the insect haemocoelium, adaptations to insect innate immunity and stress responses, and production of virulence factors. The identification of key genes in the parasitic process provides useful tools for the improvement of S. carpocapsae as a biological agent.

Population Genetics of Wild-Type CAG Repeats in the Machado-Joseph Disease Gene in Portugal

Objective: To gain insights on the molecular mechanisms of mutation that led to the emergence of expanded alleles in the MJD gene, by studying the behavior of wild-type alleles and testing the association of its distribution with the representation of the disease. Methods: The number of CAG motifs in the MJD gene was determined in a representative sample of 1000 unrelated individuals. Associations between the repeat size and the epidemiological representation of MJD were tested. Results: The allelic profile of the total sample was in the normal range (13–41 repeats), with mode (CAG)23. No intermediate alleles were present. Allelic size distribution showed a negative skew. The correlation between the epidemiological representation of MJD in each district and the frequency of small, medium and large normal alleles was not significant. Further correlations performed grouping the districts also failed to produce significant results. Conclusions: The absence of association between the size of the repeats and the representation of MJD demonstrates that prevalence is not an indirect reflection of the frequency of large normal alleles. Globally the results obtained are in accordance with a model that postulates the occurrence of a few mutations on the basis of most of the MJD cases worldwide.