Genome-wide genetic marker discovery and genotyping using next-generation sequencing.

The advent of next-generation sequencing (NGS) has revolutionized genomic and transcriptomic approaches to biology. These new sequencing tools are also valuable for the discovery, validation and assessment of genetic markers in populations. Here we review and discuss best practices for several NGS methods for genome-wide genetic marker development and genotyping that use restriction enzyme digestion of target genomes to reduce the complexity of the target. These new methods -- which include reduced-representation sequencing using reduced-representation libraries (RRLs) or complexity reduction of polymorphic sequences (CRoPS), restriction-site-associated DNA sequencing (RAD-seq) and low coverage genotyping -- are applicable to both model organisms with high-quality reference genome sequences and, excitingly, to non-model species with no existing genomic data.

A role for migration-linked genes and genomic islands in divergence of a songbird.

Next-generation sequencing has made it possible to begin asking questions about the process of divergence at the level of the genome. For example, recently, there has been a debate around the role of 'genomic islands of divergence' (i.e. blocks of outlier loci) in facilitating the process of speciation-with-gene-flow. The Swainson's thrush, Catharus ustulatus, is a migratory songbird with two genetically distinct subspecies that differ in a number of traits known to be involved in reproductive isolation in birds (plumage coloration, song and migratory behaviour), despite contemporary gene flow along a secondary contact zone. Here, we use RAD-PE sequencing to test emerging hypotheses about the process of divergence at the level of the genome and identify genes and gene regions involved in differentiation in this migratory songbird. Our analyses revealed distinct genomic islands on 15 of the 23 chromosomes and an accelerated rate of divergence on the Z chromosome, one of the avian sex chromosomes. Further, an analysis of loci linked to traits known to be involved in reproductive isolation in songbirds showed that genes linked to migration are significantly more differentiated than expected by chance, but that these genes lie primarily outside the genomic islands. Overall, our analysis supports the idea that genes linked to migration play an important role in divergence in migratory songbirds, but we find no compelling evidence that the observed genomic islands are facilitating adaptive divergence in migratory behaviour.

Genetic Analysis of Formosan Subterranean Termite (Blattodea: Rhinotermitidae) Populations in California.

A new infestation of the Formosan subterranean termite, Coptotermes formosanus Shiraki (Blattodea: Rhinotermitidae), was discovered in Canyon Lake, Riverside County, California. We used three mitochondrial DNA (COI, COII, and 16S) and seven polymorphic microsatellite markers to characterize the genetic relationship of the colony with two other colonies that were collected in 1992 and 2018 in La Mesa, San Diego County. Maximum likelihood phylogeny of C. formosanus based on concatenated COI and COII sequences revealed that the two La Mesa populations (CA01 and CA02) and the Canyon Lake population (CA03) were from different maternal lineages. Based on the 14 COII haplotypes of C. formosanus found world-wide, CA01 and CA02 belonged to a haplotype widely distributed across the United States, while CA03 was grouped under a haplotype predominantly found in Asia. Microsatellite allele frequencies across all loci for both La Mesa populations were relatively similar, but significant genetic differences were found between CA02 and CA03 colonies (FST = 0.24; Dest = 0.30; G″ ST = 0.55; P < 0.01).

Prepulse inhibition and fear-potentiated startle are altered in tissue inhibitor of metalloproteinase-2 (TIMP-2) knockout mice.

The ability to discriminate between potential dangers and recall those stimuli is essential for survival. This emotional learning requires the involvement of higher brain structures, including the amygdala, hippocampus and related cortical structures. Long-term changes in synaptic transmission and structure are important for the establishment and consolidation of fear memory. The structural changes associated with this synaptic plasticity likely require alterations in the composition of the extracellular matrix (ECM). ECM integrity is maintained by the opposing action of matrix metalloproteinases (MMPs) and their specific inhibitors, tissue inhibitors of metalloproteinases (TIMPs). To date, no studies have examined the role of MMPs or TIMPs in conditioned fear. Here, we show that neither male nor female mice deficient in TIMP-2 (knockout) exhibit prepulse inhibition of the startle reflex, suggesting deficits in pre-attentional sensorimotor gating. In addition, knockout mice and mice expressing a mutant truncated TIMP-2 (knock-down) show deficits in fear-potentiated startle. This is the first report of a phenotype for the TIMP-2(-/-) mice and suggests that TIMP-2 may play a role in the synaptic plasticity underlying learning and memory.

Restriction site-associated DNA sequencing generates high-quality single nucleotide polymorphisms for assessing hybridization between bighead and silver carp in the United States and China.

Bighead carp (Hypophthalmichthys nobilis) and silver carp (H. molitrix) are invasive species and listed as US federally injurious species under the Lacy Act. They have established populations in much of the Mississippi River Basin (MRB; Mississippi, Illinois, and Missouri rivers) and are capable of producing fertile hybrids and complex introgression. Characterizing the composition of this admixture requires a large set of high-quality, evolutionarily conserved, diagnostic genetic markers to aid in the identification and management of these species in the midst of morphological ambiguity. Restriction site-associated DNA (RAD) sequencing of 45 barcoded bighead and silver carp from the United States and China produced reads that were aligned to the silver carp transcriptome yielded 261 candidate single nucleotide polymorphisms (SNPs) with fixed allelic differences between the two species. We selected the highest quality 112 SNP loci for validation using 194 putative pure-species and F1 hybrids from the MRB and putative bighead carp and silver carp pure species from China (Amur, Pearl and Yangtze rivers). Fifty SNPs were omitted due to design/amplification failure or lack of diagnostic utility. A total of 57 species-diagnostic SNPs conserved between carp species in US and Chinese rivers were identified; 32 were annotated to functional gene loci. Twenty-seven of the 181 (15%) putative pure species were identified as hybrid backcrosses after validation, including three backcrosses from the Amur River, where hybridization has not been documented previously. The 57 SNPs identified through RAD sequencing provide a diagnostic tool to detect population admixture and to identify hybrid and pure-species Asian carps in the United States and China.

Two different high throughput sequencing approaches identify thousands of de novo genomic markers for the genetically depleted Bornean elephant.

High throughput sequencing technologies are being applied to an increasing number of model species with a high-quality reference genome. The application and analyses of whole-genome sequence data in non-model species with no prior genomic information are currently under way. Recent sequencing technologies provide new opportunities for gathering genomic data in natural populations, laying the empirical foundation for future research in the field of conservation and population genomics. Here we present the case study of the Bornean elephant, which is the most endangered subspecies of Asian elephant and exhibits very low genetic diversity. We used two different sequencing platforms, the Roche 454 FLX (shotgun) and Illumina, GAIIx (Restriction site associated DNA, RAD) to evaluate the feasibility of the two methodologies for the discovery of de novo markers (single nucleotide polymorphism, SNPs and microsatellites) using low coverage data. Approximately, 6,683 (shotgun) and 14,724 (RAD) SNPs were detected within our elephant sequence dataset. Genotyping of a representative sample of 194 SNPs resulted in a SNP validation rate of ~83 to 94% and 17% of the loci were polymorphic with a low diversity (H(o)=0.057). Different numbers of microsatellites were identified through shotgun (27,226) and RAD (868) techniques. Out of all di-, tri-, and tetra-microsatellite loci, 1,706 loci had sufficient flanking regions (shotgun) while only 7 were found with RAD. All microsatellites were monomorphic in the Bornean but polymorphic in another elephant subspecies. Despite using different sample sizes, and the well known differences in the two platforms used regarding sequence length and throughput, the two approaches showed high validation rate. The approaches used here for marker development in a threatened species demonstrate the utility of high throughput sequencing technologies as a starting point for the development of genomic tools in a non-model species and in particular for a species with low genetic diversity.

Drosophila type II neuroblast lineages keep Prospero levels low to generate large clones that contribute to the adult brain central complex.

Tissue homeostasis depends on the ability of stem cells to properly regulate self-renewal versus differentiation. Drosophila neural stem cells (neuroblasts) are a model system to study self-renewal and differentiation. Recent work has identified two types of larval neuroblasts that have different self-renewal/differentiation properties. Type I neuroblasts bud off a series of small basal daughter cells (ganglion mother cells) that each generate two neurons. Type II neuroblasts bud off small basal daughter cells called intermediate progenitors (INPs), with each INP generating 6 to 12 neurons. Type I neuroblasts and INPs have nuclear Asense and cytoplasmic Prospero, whereas type II neuroblasts lack both these transcription factors. Here we test whether Prospero distinguishes type I/II neuroblast identity or proliferation profile, using several newly characterized Gal4 lines. We misexpress prospero using the 19H09-Gal4 line (expressed in type II neuroblasts but no adjacent type I neuroblasts) or 9D11-Gal4 line (expressed in INPs but not type II neuroblasts). We find that differential prospero expression does not distinguish type I and type II neuroblast identities, but Prospero regulates proliferation in both type I and type II neuroblast lineages. In addition, we use 9D11 lineage tracing to show that type II lineages generate both small-field and large-field neurons within the adult central complex, a brain region required for locomotion, flight, and visual pattern memory.

Identification of Drosophila type II neuroblast lineages containing transit amplifying ganglion mother cells.

Mammalian neural stem cells generate transit amplifying progenitors that expand the neuronal population, but these type of progenitors have not been studied in Drosophila. The Drosophila larval brain contains approximately 100 neural stem cells (neuroblasts) per brain lobe, which are thought to bud off smaller ganglion mother cells (GMCs) that each produce two post-mitotic neurons. Here, we use molecular markers and clonal analysis to identify a novel neuroblast cell lineage containing "transit amplifying GMCs" (TA-GMCs). TA-GMCs differ from canonical GMCs in several ways: each TA-GMC has nuclear Deadpan, cytoplasmic Prospero, forms Prospero crescents at mitosis, and generates up to 10 neurons; canonical GMCs lack Deadpan, have nuclear Prospero, lack Prospero crescents at mitosis, and generate two neurons. We conclude that there are at least two types of neuroblast lineages: a Type I lineage where GMCs generate two neurons, and a type II lineage where TA-GMCs have longer lineages. Type II lineages allow more neurons to be produced faster than Type I lineages, which may be advantageous in a rapidly developing organism like Drosophila.

Drosophila Activin- and the Activin-like product Dawdle function redundantly to regulate proliferation in the larval brain.

The Drosophila Activin-like ligands Activin-beta and Dawdle control several aspects of neuronal morphogenesis, including mushroom body remodeling, dorsal neuron morphogenesis and motoneuron axon guidance. Here we show that the same two ligands act redundantly through the Activin receptor Babo and its transcriptional mediator Smad2 (Smox), to regulate neuroblast numbers and proliferation rates in the developing larval brain. Blocking this pathway results in the development of larvae with small brains and aberrant photoreceptor axon targeting, and restoring babo function in neuroblasts rescued these mutant phenotypes. These results suggest that the Activin signaling pathway is required for producing the proper number of neurons to enable normal connection of incoming photoreceptor axons to their targets. Furthermore, as the Activin pathway plays a key role in regulating propagation of mouse and human embryonic stem cells, our observation that it also regulates neuroblast numbers and proliferation in Drosophila suggests that involvement of Activins in controlling stem cell propagation may be a common regulatory feature of this family of TGF-beta-type ligands.

Nkx6 proteins specify one zebrafish primary motoneuron subtype by regulating late islet1 expression.

The ability of animals to carry out their normal behavioral repertoires requires exquisitely precise matching between specific motoneuron subtypes and the muscles they innervate. However, the molecular mechanisms that regulate motoneuron subtype specification remain unclear. Here, we use individually identified zebrafish primary motoneurons to describe a novel role for Nkx6 and Islet1 proteins in the specification of vertebrate motoneuron subtypes. We show that zebrafish primary motoneurons express two related Nkx6 transcription factors. In the absence of both Nkx6 proteins, the CaP motoneuron subtype develops normally, whereas the MiP motoneuron subtype develops a more interneuron-like morphology. In the absence of Nkx6 function, MiPs exhibit normal early expression of islet1, which is required for motoneuron formation; however, they fail to maintain islet1 expression. Misexpression of islet1 RNA can compensate for loss of Nkx6 function, providing evidence that Islet1 acts downstream of Nkx6. We suggest that Nkx6 proteins regulate MiP development at least in part by maintaining the islet1 expression that is required both to promote the MiP subtype and to suppress interneuron development.

Regulation of spindle orientation and neural stem cell fate in the Drosophila optic lobe.

BACKGROUND: The choice of a stem cell to divide symmetrically or asymmetrically has profound consequences for development and disease. Unregulated symmetric division promotes tumor formation, whereas inappropriate asymmetric division affects organ morphogenesis. Despite its importance, little is known about how spindle positioning is regulated. In some tissues cell fate appears to dictate the type of cell division, whereas in other tissues it is thought that stochastic variation in spindle position dictates subsequent sibling cell fate. RESULTS: Here we investigate the relationship between neural progenitor identity and spindle positioning in the Drosophila optic lobe. We use molecular markers and live imaging to show that there are two populations of progenitors in the optic lobe: symmetrically dividing neuroepithelial cells and asymmetrically dividing neuroblasts. We use genetically marked single cell clones to show that neuroepithelial cells give rise to neuroblasts. To determine if a change in spindle orientation can trigger a neuroepithelial to neuroblast transition, we force neuroepithelial cells to divide along their apical/basal axis by misexpressing Inscuteable. We find that this does not induce neuroblasts, nor does it promote premature neuronal differentiation. CONCLUSION: We show that symmetrically dividing neuroepithelial cells give rise to asymmetrically dividing neuroblasts in the optic lobe, and that regulation of spindle orientation and division symmetry is a consequence of cell type specification, rather than a mechanism for generating cell type diversity.