Development of novel multiplex microsatellite polymerase chain reactions to enable high-throughput population genetic studies of Schistosoma haematobium

© 2015 Webster et al. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. The attached file is the published version of the article

DNA sequence diversity and the efficiency of natural selection in animal mitochondrial DNA

Selection is expected to be more efficient in species that are more diverse because both the efficiency of natural selection and DNA sequence diversity are expected to depend upon the effective population size. We explore this relationship across a data set of 751 mammal species for which we have mitochondrial polymorphism data. We introduce a method by which we can examine the relationship between our measure of the efficiency of natural selection, the nonsynonymous relative to the synonymous nucleotide site diversity (πN/πS), and synonymous nucleotide diversity (πS), avoiding the statistical non-independence between the two quantities. We show that these two variables are strongly negatively and linearly correlated on a log scale. The slope is such that as πS doubles, πN/πS is reduced by 34%. We show that the slope of this relationship differs between the two phylogenetic groups for which we have the most data, rodents and bats, and that it also differs between species with high and low body mass, and between those with high and low mass-specific metabolic rate

60,000 years of interactions between Central and Eastern Africa documented by major African mitochondrial haplogroup L2

Mitochondrial DNA (mtDNA) haplogroup L2 originated in Western Africa but is nowadays spread across the entire continent. L2 movements were previously postulated to be related to the Bantu expansion, but L2 expansions eastwards probably occurred much earlier. By reconstructing the phylogeny of L2 (44 new complete sequences) we provide insights on the complex net of within-African migrations in the last 60 thousand years (ka). Results show that lineages in Southern Africa cluster with Western/Central African lineages at a recent time scale, whereas, eastern lineages seem to be substantially more ancient. Three moments of expansion from a Central African source are associated to L2: (1) one migration at 70-50 ka into Eastern or Southern Africa, (2) postglacial movements (15-10 ka) into Eastern Africa; and (3) the southward Bantu Expansion in the last 5 ka. The complementary population and L0a phylogeography analyses indicate no strong evidence of mtDNA gene flow between eastern and southern populations during the later movement, suggesting low admixture between Eastern African populations and the Bantu migrants. This implies that, at least in the early stages, the Bantu expansion was mainly a demic diffusion with little incorporation of local populations.This research received support from the European project “A European Initial Training Network on the History, Archaeology, and New Genetics of the Trans-Atlantic Slave Trade (EUROTAST)” (EU project: 290344). PSo is supported by FCT (the Portuguese Foundation for Science and Technology), European Social Fund, Programa Operacional Potencial Humano and the FCT Investigator Programme (IF/01641/2013). IPATIMUP integrates the i3S Research Unit, which is partially supported by FCT. This work is funded by FEDER funds through the Operational Programme for Competitiveness FactorsCOMPETE and National Funds through FCT, under the project “PEst-C/SAU/LA0003/2013”. FCT/MEC supports CBMA through Portuguese funds (PIDDAC) - PEst-OE/BIA/UI4050/2014. NORTE-07-0162FEDER-00018 (Contributos para o reforço da capacidade do IPATIMUP enquanto actor do sistema regional de inovação) and NORTE-07-0162-FEDER-000067 (Reforço e consolidação da capacidade infraestrutural do IPATIMUP para o sistema regional de inovação), both supported by Programa Operacional Regional do Norte (ON.2 – O Novo Norte), through FEDER funds under the Quadro de Referência Estratégico Nacional (QREN)

Conservation of a pH-sensitive structure in the C-terminal region of spider silk extends across the entire silk gene family

Spiders produce multiple silks with different physical properties that allow them to occupy a diverse range of ecological niches, including the underwater environment. Despite this functional diversity, past molecular analyses show a high degree of amino acid sequence similarity between C-terminal regions of silk genes that appear to be independent of the physical properties of the resulting silks; instead, this domain is crucial to the formation of silk fibres. Here we present an analysis of the C-terminal domain of all known types of spider silk and include silk sequences from the spider Argyroneta aquatica, which spins the majority of its silk underwater. Our work indicates that spiders have retained a highly conserved mechanism of silk assembly, despite the extraordinary diversification of species, silk types and applications of silk over 350 million years. Sequence analysis of the silk C-terminal domain across the entire gene family shows the conservation of two uncommon amino acids that are implicated in the formation of a salt bridge, a functional bond essential to protein assembly. This conservation extends to the novel sequences isolated from A. aquatica. This finding is relevant to research regarding the artificial synthesis of spider silk, suggesting that synthesis of all silk types will be possible using a single process

Identification of Spiroplasma insolitum symbionts in Anopheles gambiae

Background: Insect symbionts have the potential to block the transmission of vector-borne diseases by their hosts. The advancement of a symbiont-based transmission blocking strategy for malaria requires the identification and study of Anopheles symbionts. Methods: High throughput 16S amplicon sequencing was used to profile the bacteria associated with Anopheles gambiae sensu lato and identify potential symbionts. The polymerase chain reaction (PCR) with specific primers were subsequently used to monitor symbiont prevalence in field populations, as well as symbiont transmission patterns. Results: We report the discovery of the bacterial symbiont, Spiroplasma, in Anopheles gambiae in Kenya. We determine that geographically dispersed Anopheles gambiae populations in Kenya are infected with Spiroplasma at low prevalence levels. Molecular phylogenetics indicates that this Anopheles gambiae associated Spiroplasma is a member of the insolitum clade. We demonstrate that this symbiont is stably maternally transmitted across at least two generations and does not significantly affect the fecundity or egg to adult survival of its host. Conclusions: In diverse insect species, Spiroplasma has been found to render their host resistant to infection by pathogens. The identification of a maternally transmitted strain of Spiroplasma in Anopheles gambiae may therefore open new lines of investigation for the development of symbiont-based strategies for blocking malaria transmission

A new look at the LTR retrotransposon content of the chicken genome

BACKGROUND: LTR retrotransposons contribute approximately 10 % of the mammalian genome, but it has been previously reported that there is a deficit of these elements in the chicken relative to both mammals and other birds. A novel LTR retrotransposon classification pipeline, LocaTR, was developed and subsequently utilised to re-examine the chicken LTR retrotransposon annotation, and determine if the proposed chicken deficit is biologically accurate or simply a technical artefact. RESULTS: Using LocaTR 3.01 % of the chicken galGal4 genome assembly was annotated as LTR retrotransposon-derived elements (nearly double the previous annotation), including 1,073 that were structurally intact. Element distribution is significantly correlated with chromosome size and is non-random within each chromosome. Elements are significantly depleted within coding regions and enriched in gene sparse areas of the genome. Over 40 % of intact elements are found in clusters, unrelated by age or genera, generally in poorly recombining regions. The transcription of most LTR retrotransposons were suppressed or incomplete, but individual domain and full length retroviral transcripts were produced in some cases, although mostly with regularly interspersed stop codons in all reading frames. Furthermore, RNAseq data from 23 diverse tissues enabled greater characterisation of the co-opted endogenous retrovirus Ovex1. This gene was shown to be expressed ubiquitously but at variable levels across different tissues. LTR retrotransposon content was found to be very variable across the avian lineage and did not correlate with either genome size or phylogenetic position. However, the extent of previous, species-specific LTR retrotransposon annotation appears to be a confounding factor. CONCLUSIONS: Use of the novel LocaTR pipeline has nearly doubled the annotated LTR retrotransposon content of the chicken genome compared to previous estimates. Further analysis has described element distribution, clustering patterns and degree of expression in a variety of adult tissues, as well as in three embryonic stages. This study also enabled better characterisation of the co-opted gamma retroviral envelope gene Ovex1. Additionally, this work suggests that there is no deficit of LTR retrotransposons within the Galliformes relative to other birds, or to mammalian genomes when scaled for the three-fold difference in genome size

Three Thousand Years of Continuity in the Maternal Lineages of Ancient Sheep (Ovis aries) in Estonia

lthough sheep (Ovis aries) have been one of the most exploited domestic animals in Estonia since the Late Bronze Age, relatively little is known about their genetic history. Here, we explore temporal changes in Estonian sheep populations and their mitochondrial genetic diversity over the last 3000 years. We target a 558 base pair fragment of the mitochondrial hypervariable region in 115 ancient sheep from 71 sites in Estonia (c. 1200 BC – AD 1900s), 19 ancient samples from Latvia, Russia, Poland and Greece (6800 BC – AD 1700), as well as 44 samples of modern Kihnu native sheep breed. Our analyses revealed: (1) 49 mitochondrial haplotypes, associated with sheep haplogroups A and B; (2) high haplotype diversity in Estonian ancient sheep; (3) continuity in mtDNA haplotypes through time; (4) possible population expansion during the first centuries of the Middle Ages (associated with the establishment of the new power regime related to 13th century crusades); (5) significant difference in genetic diversity between ancient populations and modern native sheep, in agreement with the beginning of large-scale breeding in the 19th century and population decline in local sheep. Overall, our results suggest that in spite of the observed fluctuations in ancient sheep populations, and changes in the natural and historical conditions, the utilisation of local sheep has been constant in the territory of Estonia, displaying matrilineal continuity from the Middle Bronze Age through the Modern Period, and into modern native sheep

Genome-wide microRNA screening in Nile tilapia reveals pervasive isomiRs’ transcription, sex-biased arm switching and increasing complexity of expression throughout development

MicroRNAs (miRNAs) are key regulators of gene expression in multicellular organisms. The elucidation of miRNA function and evolution depends on the identification and characterization of miRNA repertoire of strategic organisms, as the fast-evolving cichlid fishes. Using RNA-seq and comparative genomics we carried out an in-depth report of miRNAs in Nile tilapia (Oreochromis niloticus), an emergent model organism to investigate evo-devo mechanisms. Five hundred known miRNAs and almost one hundred putative novel vertebrate miRNAs have been identified, many of which seem to be teleost-specific, cichlid-specific or tilapia-specific. Abundant miRNA isoforms (isomiRs) were identified with modifications in both 5p and 3p miRNA transcripts. Changes in arm usage (arm switching) of nine miRNAs were detected in early development, adult stage and even between male and female samples. We found an increasing complexity of miRNA expression during ontogenetic development, revealing a remarkable synchronism between the rate of new miRNAs recruitment and morphological changes. Overall, our results enlarge vertebrate miRNA collection and reveal a notable differential ratio of miRNA arms and isoforms influenced by sex and developmental life stage, providing a better picture of the evolutionary and spatiotemporal dynamics of miRNAs

A Genome Wide Association Study of arabinoxylan content in 2-row spring barley grain

In barley endosperm arabinoxylan (AX) is the second most abundant cell wall polysaccharide and in wheat it is the most abundant polysaccharide in the starchy endosperm walls of the grain. AX is one of the main contributors to grain dietary fibre content providing several health benefits including cholesterol and glucose lowering effects, and antioxidant activities. Due to its complex structural features, AX might also affect the downstream applications of barley grain in malting and brewing. Using a high pressure liquid chromatography (HPLC) method we quantified AX amounts in mature grain in 128 spring 2-row barley accessions. Amounts ranged from ~ 5.2 μg/g to ~ 9 μg/g. We used this data for a Genome Wide Association Study (GWAS) that revealed three significant quantitative trait loci (QTL) associated with grain AX levels which passed a false discovery threshold (FDR) and are located on two of the seven barley chromosomes. Regions underlying the QTLs were scanned for genes likely to be involved in AX biosynthesis or turnover, and strong candidates, including glycosyltransferases from the GT43 and GT61 families and glycoside hydrolases from the GH10 family, were identified. Phylogenetic trees of selected gene families were built based on protein translations and were used to examine the relationship of the barley candidate genes to those in other species. Our data reaffirms the roles of existing genes thought to contribute to AX content, and identifies novel QTL (and candidate genes associated with them) potentially influencing the AX content of barley grain. One potential outcome of this work is the deployment of highly associated single nucleotide polymorphisms markers in breeding programs to guide the modification of AX abundance in barley grain

Balanced gene losses, duplications and intensive rearrangements led to an unusual regularly sized genome in Arbutus unedo chloroplasts

Completely sequenced plastomes provide a valuable source of information about the duplication, loss, and transfer events of chloroplast genes and phylogenetic data for resolving relationships among major groups of plants. Moreover, they can also be useful for exploiting chloroplast genetic engineering technology. Ericales account for approximately six per cent of eudicot diversity with 11,545 species from which only three complete plastome sequences are currently available. With the aim of increasing the number of ericalean complete plastome sequences, and to open new perspectives in understanding Mediterranean plant adaptations, a genomic study on the basis of the complete chloroplast genome sequencing of Arbutus unedo and an updated phylogenomic analysis of Asteridae was implemented. The chloroplast genome of A. unedo shows extensive rearrangements but a medium size (150,897 nt) in comparison to most of angiosperms. A number of remarkable distinct features characterize the plastome of A. unedo: five-fold dismissing of the SSC region in relation to most angiosperms; complete loss or pseudogenization of a number of essential genes; duplication of the ndhH-D operon and its location within the two IRs; presence of large tandem repeats located near highly re-arranged regions and pseudogenes. All these features outline the primary evolutionary split between Ericaceae and other ericalean families. The newly sequenced plastome of A. unedo with the available asterid sequences allowed the resolution of some uncertainties in previous phylogenies of Asteridae