CHIIMP: An automated high‐throughput microsatellite genotyping platform reveals greater allelic diversity in wild chimpanzees

Short tandem repeats (STRs), also known as microsatellites, are commonly used to noninvasively genotype wild‐living endangered species, including African apes. Until recently, capillary electrophoresis has been the method of choice to determine the length of polymorphic STR loci. However, this technique is labor intensive, difficult to compare across platforms, and notoriously imprecise. Here we developed a MiSeq‐based approach and tested its performance using previously genotyped fecal samples from long‐term studied chimpanzees in Gombe National Park, Tanzania. Using data from eight microsatellite loci as a reference, we designed a bioinformatics platform that converts raw MiSeq reads into locus‐specific files and automatically calls alleles after filtering stutter sequences and other PCR artifacts. Applying this method to the entire Gombe population, we confirmed previously reported genotypes, but also identified 31 new alleles that had been missed due to sequence differences and size homoplasy. The new genotypes, which increased the allelic diversity and heterozygosity in Gombe by 61% and 8%, respectively, were validated by replicate amplification and pedigree analyses. This demonstrated inheritance and resolved one case of an ambiguous paternity. Using both singleplex and multiplex locus amplification, we also genotyped fecal samples from chimpanzees in the Greater Mahale Ecosystem in Tanzania, demonstrating the utility of the MiSeq‐based approach for genotyping nonhabituated populations and performing comparative analyses across field sites. The new automated high‐throughput analysis platform (available at https://github.com/ShawHahnLab/chiimp ) will allow biologists to more accurately and effectively determine wildlife population size and structure, and thus obtain information critical for conservation efforts.     

High‐throughput sequencing and graph‐based cluster analysis facilitate microsatellite development from a highly complex genome

Despite recent advances in high‐throughput sequencing, difficulties are often encountered when developing microsatellites for species with large and complex genomes. This probably reflects the close association in many species of microsatellites with cryptic repetitive elements. We therefore developed a novel approach for isolating polymorphic microsatellites from the club‐legged grasshopper (Gomphocerus sibiricus), an emerging quantitative genetic and behavioral model system. Whole genome shotgun Illumina MiSeq sequencing was used to generate over three million 300 bp paired‐end reads, of which 67.75% were grouped into 40,548 clusters within RepeatExplorer. Annotations of the top 468 clusters, which represent 60.5% of the reads, revealed homology to satellite DNA and a variety of transposable elements. Evaluating 96 primer pairs in eight wild‐caught individuals, we found that primers mined from singleton reads were six times more likely to amplify a single polymorphic microsatellite locus than primers mined from clusters. Our study provides experimental evidence in support of the notion that microsatellites associated with repetitive elements are less likely to successfully amplify. It also reveals how advances in high‐throughput sequencing and graph‐based repetitive DNA analysis can be leveraged to isolate polymorphic microsatellites from complex genomes.     

A proposal for standardization in forensic equine DNA typing: allele nomenclature for 17 equine‐specific STR loci

In this study, a proposal is presented for the allele nomenclature of 17 polymorphic STR loci (AHT4, AHT5, ASB2, ASB17, ASB23, CA425, HMS1, HMS2, HMS3, HMS6, HMS7, HTG4, HTG6, HTG7, HTG10, LEX3 and VHL20) for equine genotyping (Equus caballus). The nomenclature is based on sequence data of the polymorphic region of the STR loci as recommended by the DNA commission of the International Society for Forensic Genetics for human DNA typing. For each STR locus, several alleles were selected and animals homozygous for those alleles were subjected to sequence analysis. The alleles of the 17 STR loci consisted either of simple (10), compound (6) or complex repeat patterns (1). Only a limited number of alleles with the same fragment size showed different repeat structures. The allele designation described here was based on the number of repeats, including all variable regions within the amplified fragment.     

Helping decision making for reliable and cost‐effective 2b‐RAD sequencing and genotyping analyses in non‐model species

High‐throughput sequencing has revolutionized population and conservation genetics. RAD sequencing methods, such as 2b‐RAD, can be used on species lacking a reference genome. However, transferring protocols across taxa can potentially lead to poor results. We tested two different IIB enzymes (AlfI and CspCI) on two species with different genome sizes (the loggerhead turtle Caretta caretta and the sharpsnout seabream Diplodus puntazzo) to build a set of guidelines to improve 2b‐RAD protocols on non‐model organisms while optimising costs. Good results were obtained even with degraded samples, showing the value of 2b‐RAD in studies with poor DNA quality. However, library quality was found to be a critical parameter on the number of reads and loci obtained for genotyping. Resampling analyses with different number of reads per individual showed a trade‐off between number of loci and number of reads per sample. The resulting accumulation curves can be used as a tool to calculate the number of sequences per individual needed to reach a mean depth ≥20 reads to acquire good genotyping results. Finally, we demonstrated that selective‐base ligation does not affect genomic differentiation between individuals, indicating that this technique can be used in species with large genome sizes to adjust the number of loci to the study scope, to reduce sequencing costs and to maintain suitable sequencing depth for a reliable genotyping without compromising the results. Here, we provide a set of guidelines to improve 2b‐RAD protocols on non‐model organisms with different genome sizes, helping decision‐making for a reliable and cost‐effective genotyping.     

Seasonal diversity and dynamics of haptophytes in the Skagerrak,Norway, explored by high‐throughput sequencing

Microalgae in the division Haptophyta play key roles in the marine ecosystem and in global biogeochemical processes. Despite their ecological importance, knowledge on seasonal dynamics, community composition and abundance at the species level is limited due to their small cell size and few morphological features visible under the light microscope. Here, we present unique data on haptophyte seasonal diversity and dynamics from two annual cycles, with the taxonomic resolution and sampling depth obtained with high‐throughput sequencing. From outer Oslofjorden, S Norway, nano‐ and picoplanktonic samples were collected monthly for 2 years, and the haptophytes targeted by amplification of RNA/cDNA with Haptophyta‐specific 18S rDNA V4 primers. We obtained 156 operational taxonomic units (OTUs), from c. 400.000 454 pyrosequencing reads, after rigorous bioinformatic filtering and clustering at 99.5%. Most OTUs represented uncultured and/or not yet 18S rDNA‐sequenced species. Haptophyte OTU richness and community composition exhibited high temporal variation and significant yearly periodicity. Richness was highest in September–October (autumn) and lowest in April–May (spring). Some taxa were detected all year, such as Chrysochromulina simplex, Emiliania huxleyi and Phaeocystis cordata, whereas most calcifying coccolithophores only appeared from summer to early winter. We also revealed the seasonal dynamics of OTUs representing putative novel classes (clades HAP‐3–5) or orders (clades D, E, F). Season, light and temperature accounted for 29% of the variation in OTU composition. Residual variation may be related to biotic factors, such as competition and viral infection. This study provides new, in‐depth knowledge on seasonal diversity and dynamics of haptophytes in North Atlantic coastal waters.     

QDD version 3.1: a user‐friendly computer program for microsatellite selection and primer design revisited: experimental validation of variables determining genotyping success rate

Microsatellite marker development has been greatly simplified by the use of high‐throughput sequencing followed by in silico microsatellite detection and primer design. However, the selection of markers designed by the existing pipelines depends either on arbitrary criteria, or older studies on PCR success. Based on wet laboratory experiments, we have identified the following factors that are most likely to influence genotyping success rate: alignment score between the primers and the amplicon; the distance between primers and microsatellites; the length of the PCR product; target region complexity and the number of reads underlying the sequence. The QDD pipeline has been modified to include these most pertinent factors in the output to help the selection of markers. Furthermore, new features are also included in the present version: (i) not only raw sequencing reads are accepted as input, but also contigs, allowing the analysis of assembled high‐coverage data; (ii) input data can be both in fasta and fastq format to facilitate the use of Illumina and IonTorrent reads; (iii) A comparison to known transposable elements allows their detection; (iv) A contamination check can be carried out by BLASTing potential markers against the nucleotide (nt) database of NCBI; (v) QDD3 is now also available imbedded into a virtual machine making installation easier and operating system independent. It can be used both on command‐line version as well as integrated into a Galaxy server, providing a user‐friendly interface, as well as the possibility to utilize a large variety of NGS tools.     

Diversity and abundance of human‐pathogenic fungi associated with pigeon faeces in urban environments

Pathogenic fungi are a growing health concern worldwide, particularly in large, densely populated cities. The dramatic upsurge of pigeon populations in cities has been implicated in the increased incidence of invasive fungal infections. In this study, we used a culture‐independent, high‐throughput sequencing approach to describe the diversity of clinically relevant fungi (CRF) associated with pigeon faeces and map the relative abundance of CRF across Seoul, Korea. In addition, we tested whether certain geographical, sociological and meteorological factors were significantly associated with the diversity and relative abundance of CRF. Finally, we compared the CRF diversity of fresh and old pigeon faeces to identify the source of the fungi and the role of pigeons in dispersal. Our results demonstrated that both the composition and relative abundance of CRF are unevenly distributed across Seoul. The green area ratio and the number of multiplex houses were positively correlated with species diversity, whereas wind speed and number of households were negatively correlated. The number of workers and green area ratio were positively correlated with the relative abundance of CRF, whereas wind speed was negatively correlated. Because many CRF were absent in fresh faeces, we inferred that most species cannot survive the gastrointestinal tract of pigeons and instead are likely transmitted through soil or air and use pigeon faeces as a substrate for proliferation.     

The decline of the Turtle Dove: Dietary associations with body condition and competition with other columbids analysed using high‐throughput sequencing

Dietary changes linked to the availability of anthropogenic food resources can have complex implications for species and ecosystems, especially when species are in decline. Here, we use recently developed primers targeting the ITS2 region of plants to characterize diet from faecal samples of four UK columbids, with particular focus on the European turtle dove (Streptopelia turtur), a rapidly declining obligate granivore. We examine dietary overlap between species (potential competition), associations with body condition in turtle doves and spatiotemporal variation in diet. We identified 143 taxonomic units, of which we classified 55% to species, another 34% to genus and the remaining 11% to family. We found significant dietary overlap between all columbid species, with the highest between turtle doves and stock doves (Columba oenas), then between turtle doves and woodpigeons (Columba palumbus). The lowest overlap was between woodpigeons and collared doves (Streptopelia decaocto). We show considerable change in columbid diets compared to previous studies, probably reflecting opportunistic foraging behaviour by columbids within a highly anthropogenically modified landscape, although our data for nonturtle doves should be considered preliminary. Nestling turtle doves in better condition had a higher dietary proportion of taxonomic units from natural arable plant species and a lower proportion of taxonomic units from anthropogenic food resources such as garden bird seed mixes and brassicas. This suggests that breeding ground conservation strategies for turtle doves should include provision of anthropogenic seeds for adults early in the breeding season, coupled with habitat rich in accessible seeds from arable plants once chicks have hatched.     

A robust,high‐throughput method for computing maize ear,cob, and kernel attributes automatically from images

Grain yield of the maize plant depends on the sizes, shapes, and numbers of ears and the kernels they bear. An automated pipeline that can measure these components of yield from easily‐obtained digital images is needed to advance our understanding of this globally important crop. Here we present three custom algorithms designed to compute such yield components automatically from digital images acquired by a low‐cost platform. One algorithm determines the average space each kernel occupies along the cob axis using a sliding‐window Fourier transform analysis of image intensity features. A second counts individual kernels removed from ears, including those in clusters. A third measures each kernel's major and minor axis after a Bayesian analysis of contour points identifies the kernel tip. Dimensionless ear and kernel shape traits that may interrelate yield components are measured by principal components analysis of contour point sets. Increased objectivity and speed compared to typical manual methods are achieved without loss of accuracy as evidenced by high correlations with ground truth measurements and simulated data. Millimeter‐scale differences among ear, cob, and kernel traits that ranged more than 2.5‐fold across a diverse group of inbred maize lines were resolved. This system for measuring maize ear, cob, and kernel attributes is being used by multiple research groups as an automated Web service running on community high‐throughput computing and distributed data storage infrastructure. Users may create their own workflow using the source code that is staged for download on a public repository.     

A likelihood‐based approach for assessment of extra‐pair paternity and conspecific brood parasitism in natural populations

Genotypes are frequently used to assess alternative reproductive strategies such as extra‐pair paternity and conspecific brood parasitism in wild populations. However, such analyses are vulnerable to genotyping error or molecular artefacts that can bias results. For example, when using multilocus microsatellite data, a mismatch at a single locus, suggesting the offspring was not directly related to its putative parents, can occur quite commonly even when the offspring is truly related. Some recent studies have advocated an ad‐hoc rule that offspring must differ at more than one locus in order to conclude that they are not directly related. While this reduces the frequency with which true offspring are identified as not directly related young, it also introduces bias in the opposite direction, wherein not directly related young are categorized as true offspring. More importantly, it ignores the additional information on allele frequencies which would reduce overall bias. In this study, we present a novel technique for assessing extra‐pair paternity and conspecific brood parasitism using a likelihood‐based approach in a new version of program cervus . We test the suitability of the technique by applying it to a simulated data set and then present an example to demonstrate its influence on the estimation of alternative reproductive strategies.     

Characterization of repeat arrays in ultra‐long nanopore reads reveals frequent origin of satellite DNA from retrotransposon‐derived tandem repeats

Amplification of monomer sequences into long contiguous arrays is the main feature distinguishing satellite DNA from other tandem repeats, yet it is also the main obstacle in its investigation because these arrays are in principle difficult to assemble. Here we explore an alternative, assembly‐free approach that utilizes ultra‐long Oxford Nanopore reads to infer the length distribution of satellite repeat arrays, their association with other repeats and the prevailing sequence periodicities. Using the satellite DNA‐rich legume plant Lathyrus sativus as a model, we demonstrated this approach by analyzing 11 major satellite repeats using a set of nanopore reads ranging from 30 to over 200 kb in length and representing 0.73× genome coverage. We found surprising differences between the analyzed repeats because only two of them were predominantly organized in long arrays typical for satellite DNA. The remaining nine satellites were found to be derived from short tandem arrays located within LTR‐retrotransposons that occasionally expanded in length. While the corresponding LTR‐retrotransposons were dispersed across the genome, this array expansion occurred mainly in the primary constrictions of the L. sativus chromosomes, which suggests that these genome regions are favourable for satellite DNA accumulation.     

Downregulation of the inflammatory network in senescent fibroblasts and aging tissues of the long‐lived and cancer‐resistant subterranean wild rodent,Spalax

The blind mole rat (Spalax) is a wild, long‐lived rodent that has evolved mechanisms to tolerate hypoxia and resist cancer. Previously, we demonstrated high DNA repair capacity and low DNA damage in Spalax fibroblasts following genotoxic stress compared with rats. Since the acquisition of senescence‐associated secretory phenotype (SASP) is a consequence of persistent DNA damage, we investigated whether cellular senescence in Spalax is accompanied by an inflammatory response. Spalax fibroblasts undergo replicative senescence (RS) and etoposide‐induced senescence (EIS), evidenced by an increased activity of senescence‐associated beta‐galactosidase (SA‐β‐Gal), growth arrest, and overexpression of p21, p16, and p53 mRNAs. Yet, unlike mouse and human fibroblasts, RS and EIS Spalax cells showed undetectable or decreased expression of the well‐known SASP factors: interleukin‐6 (IL6), IL8, IL1α, growth‐related oncogene alpha (GROα), SerpinB2, and intercellular adhesion molecule (ICAM‐1). Apparently, due to the efficient DNA repair in Spalax, senescent cells did not accumulate the DNA damage necessary for SASP activation. Conversely, Spalax can maintain DNA integrity during replicative or moderate genotoxic stress and limit pro‐inflammatory secretion. However, exposure to the conditioned medium of breast cancer cells MDA‐MB‐231 resulted in an increase in DNA damage, activation of the nuclear factor κB (NF‐κB) through nuclear translocation, and expression of inflammatory mediators in RS Spalax cells. Evaluation of SASP in aging Spalax brain and intestine confirmed downregulation of inflammatory‐related genes. These findings suggest a natural mechanism for alleviating the inflammatory response during cellular senescence and aging in Spalax, which can prevent age‐related chronic inflammation supporting healthy aging and longevity.     

Simple embryo injection of long single‐stranded donor templates with the CRISPR/Cas9 system leads to homology‐directed repair in Xenopus tropicalis and Xenopus laevis

We report model experiments in which simple microinjection of fertilized eggs has been used to effectively perform homology‐directed repair (HDR)‐mediated gene editing in the two Xenopus species used most frequently for research: X. tropicalis and X. laevis. We have used long single‐stranded DNAs having phosphorothioate modifications as donor templates for HDR at targeted genomic sites using the Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR‐associated protein 9 (CRISPR/Cas9) system. First, X. tropicalis tyr mutant (i.e., albino) embryos were successfully rescued: partially pigmented tadpoles were seen in up to 35% of injected embryos, demonstrating the potential for efficient insertion of targeted point mutations. Second, in order to demonstrate the ability to tag genes with fluorescent proteins (FPs), we targeted the melanocyte‐specific gene slc45a2.L of X. laevis to label it with the Superfolder green FP (sfGFP), seeing mosaic expression of sfGFP in melanophores in up to 20% of injected tadpoles. Tadpoles generated by these two approaches were raised to sexual maturity, and shown to successfully transmit HDR constructs through the germline with precise targeting and seamless recombination. F1 embryos showed rescue of the tyr mutation (X. tropicalis) and tagging in the appropriate pigment cell‐specific manner of slc45a2.L with sfGFP (X. laevis).     

Noninvasive individual and species identification of jaguars (Panthera onca), pumas (Puma concolor) and ocelots (Leopardus pardalis) in Belize,Central America using cross‐species microsatellites and faecal DNA

There is a great need to develop efficient, noninvasive genetic sampling methods to study wild populations of multiple, co‐occurring, threatened felids. This is especially important for molecular scatology studies occurring in challenging tropical environments where DNA degrades quickly and the quality of faecal samples varies greatly. We optimized 14 polymorphic microsatellite loci for jaguars (Panthera onca), pumas (Puma concolor) and ocelots (Leopardus pardalis) and assessed their utility for cross‐species amplification. Additionally, we tested their reliability for species and individual identification using DNA from faeces of wild felids detected by a scat detector dog across Belize in Central America. All microsatellite loci were successfully amplified in the three target species, were polymorphic with average expected heterozygosities of H_E = 0.60 ± 0.18 (SD) for jaguars, H_E = 0.65 ± 0.21 (SD) for pumas and H_E = 0.70 ± 0.13 (SD) for ocelots and had an overall PCR amplification success of 61%. We used this nuclear DNA primer set to successfully identify species and individuals from 49% of 1053 field‐collected scat samples. This set of optimized microsatellite multiplexes represents a powerful tool for future efforts to conduct noninvasive studies on multiple, wild Neotropical felids.     

RAD‐seq linkage mapping and patterns of segregation distortion in sedges: meiosis as a driver of karyotypic evolution in organisms with holocentric chromosomes

Meiotic drive, the class of meiotic mechanisms that drive unequal segregation of alleles among gametes, may be an important force in karyotype evolution. Its role in holocentric organisms, whose chromosomes lack localized centromeres, is poorly understood. We crossed two individuals of Carex scoparia (Cyperaceae) with different chromosome numbers (2n = 33^II = 66 × 2n = 32^II = 64) to obtain F1 individuals, which we then self‐pollinated to obtain second‐generation (F2) crosses. RAD‐seq was performed for 191 individuals (including the parents, five F1 individuals and 184 F2 individuals). Our F2 linkage map based on stringent editing of the RAD‐seq data set yielded 32 linkage groups. In the final map, 865 loci were located on a linkage map of 3966.99 cM (linkage groups ranged from 24.39 to 193.31 cM in length and contained 5–51 loci each). Three linkage groups exhibit more loci under segregation distortion than expected by chance; within linkage groups, loci exhibiting segregation distortion are clustered. This finding implicates meiotic drive in the segregation of chromosome variants, suggesting that selection of chromosome variants in meiosis may contribute to the establishment and fixation of chromosome variants in Carex, which is renowned for high chromosomal and species diversity. This is an important finding as previous studies demonstrate that chromosome divergence may play a key role in differentiation and speciation in Carex.     

What remains from a 454 run: estimation of success rates of microsatellite loci development in selected newt species (Calotriton asper,Lissotriton helveticus,and Triturus cristatus) and comparison with Illumina‐based approaches

The development of microsatellite loci has become more efficient using next‐generation sequencing (NGS) approaches, and many studies imply that the amount of applicable loci is large. However, few studies have sought to quantify the number of loci that are retained for use out of the thousands of sequence reads initially obtained. We analyzed the success rate of microsatellite loci development for three amphibian species using a 454 NGS approach on tetra‐nucleotide motif‐enriched species‐specific libraries. The number of sequence reads obtained differed strongly between species and ranged from 19,562 for Triturus cristatus to 55,626 for Lissotriton helveticus, with 52,075 reads obtained for Calotriton asper. PHOBOS was used to identify sequences with tetra‐nucleotide repeat motifs with a minimum repeat number of ten and high quality primer binding sites. Of 107 sequences for T. cristatus, 316 for C. asper and 319 for L. helveticus, we tested the amplification success, polymorphism, and degree of heterozygosity for 41 primer combinations each for C. asper and T. cristatus, and 22 for L. helveticus. We found 11 polymorphic loci for T. cristatus, 20 loci for C. asper, and 15 loci for L. helveticus. Extrapolated, the number of potentially amplifiable loci (PALs) resulted in estimated species‐specific success rates of 0.15% (T. cristatus), 0.30% (C. asper), and 0.39% (L. helveticus). Compared with representative Illumina NGS approaches, our applied 454‐sequencing approach on specifically enriched sublibraries proved to be quite competitive in terms of success rates and number of finally applicable loci.     

Weight under Steel Wheels: Material Stock and Flow Analysis of High‐Speed Rail in China

The construction of a nation‐wide high‐speed rail (HSR) network has emerged as a hugely expensive and ambitious infrastructure project in China. As of December 2012, some 8,800 kilometers (km) of double‐track HSR lines came into service in the country, accounting for 40% of the total HSR length in the world. The network is expected to expand to 34,000 km or longer in around two decades. As the first HSR system specially built and operated in an economically developing country, it helps integrate the sprawling economy and lift the quality of life of the increasing urban population. China's experiences in HSR are expected to be of value to other countries aiming to adopt bullet train systems, especially those at a similar level of industrialization and urbanization. This work specifically examines material stocks and flows associated with the HSR infrastructure construction in China. A major distinction from the construction of HSR tracks in Europe is that nearly 70% of the HSR tracks in China are laid upon bridges or inside tunnels, which are structures that demand great amounts of raw materials. The entire network, once completed by 2030, will cumulatively require 83 to 137 million tonnes (Mt) of steel and 560 to 920 Mt of cement. This is still a small share of China's use of material resources. Nonetheless, the massive application of the steel‐ and cement‐intensive structures deserves consideration when assessing the environmental performance of HSR over its entire life cycle.