Benchmarking variant identification tools for plant diversity discovery.

BACKGROUND: The ability to accurately and comprehensively identify genomic variations is critical for plant studies utilizing high-throughput sequencing. Most bioinformatics tools for processing next-generation sequencing data were originally developed and tested in human studies, raising questions as to their efficacy for plant research. A detailed evaluation of the entire variant calling pipeline, including alignment, variant calling, variant filtering, and imputation was performed on different programs using both simulated and real plant genomic datasets. RESULTS: A comparison of SOAP2, Bowtie2, and BWA-MEM found that BWA-MEM was consistently able to align the most reads with high accuracy, whereas Bowtie2 had the highest overall accuracy. Comparative results of GATK HaplotypCaller versus SAMtools mpileup indicated that the choice of variant caller affected precision and recall differentially depending on the levels of diversity, sequence coverage and genome complexity. A cross-reference experiment of S. lycopersicum and S. pennellii reference genomes revealed the inadequacy of single reference genome for variant discovery that includes distantly-related plant individuals. Machine-learning-based variant filtering strategy outperformed the traditional hard-cutoff strategy resulting in higher number of true positive variants and fewer false positive variants. A 2-step imputation method, which utilized a set of high-confidence SNPs as the reference panel, showed up to 60% higher accuracy than direct LD-based imputation. CONCLUSIONS: Programs in the variant discovery pipeline have different performance on plant genomic dataset. Choice of the programs is subjected to the goal of the study and available resources. This study serves as an important guiding information for plant biologists utilizing next-generation sequencing data for diversity characterization and crop improvement.

In situ embryo rescue for generation of wide intra- and interspecific hybrids of Panicum virgatum L.

Wide crosses have been used for decades as a method for transferring novel genetic material and traits in plant breeding. Historically, many products of wide crosses require tedious and inefficient surgical embryo rescue prior to embryo abortion to recover single plantlets. We have utilized transgenic switchgrass (Panicum virgatum L. cv Alamo) as a pollen donor in conjunction with antibiotic or herbicide selection for recovery of intra-and interspecific F1 crosses by using developing ovules from the female parent and selecting for embryogenic cultures derived from the in situ immature embryo. Using this approach, several intravarietial crosses were generated between transgenic Alamo and the switchgrass varieties Kanlow, Blackwell and Cave-in-Rock as well as an interspecific cross with Atlantic coastal panicgrass. This procedure selected F1 embryogenic callus produced from the developing embryo contained within isolated immature ovules. Several clonal plants were successfully regenerated from each cross. Southern blot, PCR, phenotypic analyses and genomic analysis confirmed F1 hybrids. Using genotyping-by-sequencing shows the hybridization of the recovered plants by determining the ratio of transgressive markers to total compared markers between parents and their potential offspring. The ratio of transgressive markers to total compared markers was significantly lower between parents and their predicted offspring than between parents and offspring unrelated to them. This approach provides the possibility to move useful transgenes into varieties that are recalcitrant to direct transformation which can be optionally segregated thus useful to create new hybrids, as well as recovery of wide crosses that are either difficult or impossible using traditional techniques.

Flexible and scalable genotyping-by-sequencing strategies for population studies.

BACKGROUND: Many areas critical to agricultural production and research, such as the breeding and trait mapping in plants and livestock, require robust and scalable genotyping platforms. Genotyping-by-sequencing (GBS) is a one such method highly suited to non-human organisms. In the GBS protocol, genomic DNA is fractionated via restriction digest, then reduced representation is achieved through size selection. Since many restriction sites are conserved across a species, the sequenced portion of the genome is highly consistent within a population. This makes the GBS protocol highly suited for experiments that require surveying large numbers of markers within a population, such as those involving genetic mapping, breeding, and population genomics. We have modified the GBS technology in a number of ways. Custom, enzyme specific adaptors have been replaced with standard Illumina adaptors compatible with blunt-end restriction enzymes. Multiplexing is achieved through a dual barcoding system, and bead-based library preparation protocols allows for in-solution size selection and eliminates the need for columns and gels. RESULTS: A panel of eight restriction enzymes was selected for testing on B73 maize and Nipponbare rice genomic DNA. Quality of the data was demonstrated by identifying that the vast majority of reads from each enzyme aligned to restriction sites predicted in silico. The link between enzyme parameters and experimental outcome was demonstrated by showing that the sequenced portion of the genome was adaptable by selecting enzymes based on motif length, complexity, and methylation sensitivity. The utility of the new GBS protocol was demonstrated by correctly mapping several in a maize F2 population resulting from a B73×Country Gentleman test cross. CONCLUSIONS: This technology is readily adaptable to different genomes, highly amenable to multiplexing and compatible with over forty commercially available restriction enzymes. These advancements represent a major improvement in genotyping technology by providing a highly flexible and scalable GBS that is readily implemented for studies on genome-wide variation.

The Trichoplax genome and the nature of placozoans.

As arguably the simplest free-living animals, placozoans may represent a primitive metazoan form, yet their biology is poorly understood. Here we report the sequencing and analysis of the approximately 98 million base pair nuclear genome of the placozoan Trichoplax adhaerens. Whole-genome phylogenetic analysis suggests that placozoans belong to a 'eumetazoan' clade that includes cnidarians and bilaterians, with sponges as the earliest diverging animals. The compact genome shows conserved gene content, gene structure and synteny in relation to the human and other complex eumetazoan genomes. Despite the apparent cellular and organismal simplicity of Trichoplax, its genome encodes a rich array of transcription factor and signalling pathway genes that are typically associated with diverse cell types and developmental processes in eumetazoans, motivating further searches for cryptic cellular complexity and/or as yet unobserved life history stages.

Evolutionary genetics of the hydroid allodeterminant alr2.

We surveyed genetic variation in alr2, an allodeterminant of the colonial hydroid Hydractinia symbiolongicarpus. We generated cDNA from a sample of 239 Hydractinia colonies collected at Lighthouse Point, Connecticut, and identified 473 alr2 alleles, 198 of which were unique. Rarefaction analysis suggested that the sample was near saturation. Most alleles were rare, with 86% occurring at frequencies of 1% or less. Alleles were highly variable, diverging on average by 18% of the amino acids in a predicted extracellular domain of the molecule. Analysis of 152 full-length alleles confirmed the existence of two structural types, defined by exons 4-8 of the gene. Several residues of the predicted immunoglobulin superfamily-like domains display signatures of positive selection. We also identified 77 unique alr2 pseudogene sequences from 85 colonies. Twenty-seven of these sequences matched expressed alr2 sequences from other colonies. This observation is consistent with pseudogenes contributing to alr2 diversification through sequence donation. A more limited collection of animals was made from a distant, relict population of H. symbiolongicarpus. Sixty percent of the unique sequences identified in this sample were found to match sequences from the Lighthouse Point population. The large number of alr2 alleles, their degree of divergence, the predominance of rare alleles in the population, their persistence over broad spatial and temporal scales, and the signatures of positive selection in multiple residues of the putative recognition domain paint a consistent picture of negative-frequency-dependent selection operating in this system. The genetic diversity observed at alr2 is comparable to that of the most highly polymorphic genetic systems known to date.

Genetic diversity of the allodeterminant alr2 in Hydractinia symbiolongicarpus.

Hydractinia symbiolongicarpus, a colonial cnidarian (class Hydrozoa) epibiont on hermit crab shells, is well established as a model for genetic studies of allorecognition. Recently, two linked loci, allorecognition (alr) 1 and alr2, were identified by positional cloning and shown to be major determinants of histocompatibility. Both genes encode putative transmembrane proteins with hypervariable extracellular domains similar to immunoglobulin (Ig)-like domains. We sought to characterize the naturally occurring variation at the alr2 locus and to understand the origins of this molecular diversity. We examined full-length cDNA coding sequences derived from a sample of 21 field-collected colonies, including 18 chosen haphazardly and two laboratory reference strains. Of the 35 alleles recovered from the 18 unbiased samples, 34 encoded unique gene products. We identified two distinct structural classes of alleles that varied over a large central region of the gene but both possessed highly polymorphic extracellular domains I, similar to an Ig-like V-set domain. The discovery of structurally chimeric alleles provided evidence that interallelic recombination may contribute to alr2 variation. Comparisons of the genomic region encompassing alr2 from two field-derived haplotypes and one laboratory reference sequence revealed a history of structural variation at the haplotype level as well. Maintenance of large numbers of equally rare alleles in a natural population is a hallmark of negative frequency-dependent selection and is expected to produce high levels of heterozygosity. The observed alr2 allelic diversity is comparable with that found in immune recognition molecules such as human leukocyte antigens, B cell Igs, or natural killer cell Ig-like receptors.

Concatenated analysis sheds light on early metazoan evolution and fuels a modern "urmetazoon" hypothesis.

For more than a century, the origin of metazoan animals has been debated. One aspect of this debate has been centered on what the hypothetical "urmetazoon" bauplan might have been. The morphologically most simply organized metazoan animal, the placozoan Trichoplax adhaerens, resembles an intriguing model for one of several "urmetazoon" hypotheses: the placula hypothesis. Clear support for a basal position of Placozoa would aid in resolving several key issues of metazoan-specific inventions (including, for example, head-foot axis, symmetry, and coelom) and would determine a root for unraveling their evolution. Unfortunately, the phylogenetic relationships at the base of Metazoa have been controversial because of conflicting phylogenetic scenarios generated while addressing the question. Here, we analyze the sum of morphological evidence, the secondary structure of mitochondrial ribosomal genes, and molecular sequence data from mitochondrial and nuclear genes that amass over 9,400 phylogenetically informative characters from 24 to 73 taxa. Together with mitochondrial DNA genome structure and sequence analyses and Hox-like gene expression patterns, these data (1) provide evidence that Placozoa are basal relative to all other diploblast phyla and (2) spark a modernized "urmetazoon" hypothesis.

Social network analysis of the genealogy of strawberry: retracing the wild roots of heirloom and modern cultivars.

The widely recounted story of the origin of cultivated strawberry (Fragaria × ananassa) oversimplifies the complex interspecific hybrid ancestry of the highly admixed populations from which heirloom and modern cultivars have emerged. To develop deeper insights into the three-century-long domestication history of strawberry, we reconstructed the genealogy as deeply as possible-pedigree records were assembled for 8,851 individuals, including 2,656 cultivars developed since 1775. The parents of individuals with unverified or missing pedigree records were accurately identified by applying an exclusion analysis to array-genotyped single-nucleotide polymorphisms. We identified 187 wild octoploid and 1,171 F. × ananassa founders in the genealogy, from the earliest hybrids to modern cultivars. The pedigree networks for cultivated strawberry are exceedingly complex labyrinths of ancestral interconnections formed by diverse hybrid ancestry, directional selection, migration, admixture, bottlenecks, overlapping generations, and recurrent hybridization with common ancestors that have unequally contributed allelic diversity to heirloom and modern cultivars. Fifteen to 333 ancestors were predicted to have transmitted 90% of the alleles found in country-, region-, and continent-specific populations. Using parent-offspring edges in the global pedigree network, we found that selection cycle lengths over the past 200 years of breeding have been extraordinarily long (16.0-16.9 years/generation), but decreased to a present-day range of 6.0-10.0 years/generation. Our analyses uncovered conspicuous differences in the ancestry and structure of North American and European populations, and shed light on forces that have shaped phenotypic diversity in F. × ananassa.

Comparative genomics of large mitochondria in placozoans.

The first sequenced mitochondrial genome of a placozoan, Trichoplax adhaerens, challenged the conventional wisdom that a compact mitochondrial genome is a common feature among all animals. Three additional placozoan mitochondrial genomes representing highly divergent clades have been sequenced to determine whether the large Trichoplax mtDNA is a shared feature among members of the phylum Placozoa or a uniquely derived condition. All three mitochondrial genomes were found to be very large, 32- to 37-kb, circular molecules, having the typical 12 respiratory chain genes, 24 tRNAs, rnS, and rnL. They share with the Trichoplax mitochondrial genome the absence of atp8, atp9, and all ribosomal protein genes, the presence of several cox1 introns, and a large open reading frame containing an intron group I LAGLIDADG endonuclease domain. The differences in mtDNA size within Placozoa are due to variation in intergenic spacer regions and the presence or absence of long open reading frames of unknown function. Phylogenetic analyses of the 12 respiratory chain genes support the monophyly of Placozoa. The similarities in composition and structure between the three mitochondrial genomes reported here and that of Trichoplax's mtDNA suggest that their uncompacted state is a shared ancestral feature to other nonmetazoans while their gene content is a derived feature shared only among the Metazoa.

Axial patterning and diversification in the cnidaria predate the Hox system.

Across the animal kingdom, Hox genes are organized in clusters whose genomic organization reflects their central roles in patterning along the anterior/posterior (A/P) axis . While a cluster of Hox genes was present in the bilaterian common ancestor, the origins of this system remain unclear (cf. ). With new data for two representatives of the closest extant phylum to the Bilateria, the sea anemone Nematostella and the hydromedusa Eleutheria, we argue here that the Cnidaria predate the evolution of the Hox system. Although Hox-like genes are present in a range of cnidarians, many of these are paralogs and in neither Nematostella nor Eleutheria is an equivalent of the Hox cluster present. With the exception of independently duplicated genes, the cnidarian genes are unlinked and in several cases are flanked by non-Hox genes. Furthermore, the cnidarian genes are expressed in patterns that are inconsistent with the Hox paradigm. We conclude that the Cnidaria/Bilateria split occurred before a definitive Hox system developed. The spectacular variety in morphological and developmental characteristics shown by extant cnidarians demonstrates that there is no obligate link between the Hox system and morphological diversity in the animal kingdom and that a canonical Hox system is not mandatory for axial patterning.

The mitochondrial genome of the hexactinellid sponge Aphrocallistes vastus: evidence for programmed translational frameshifting.

BACKGROUND: Mitochondrial genomes (mtDNA) of numerous sponges have been sequenced as part of an ongoing effort to resolve the class-level phylogeny of the Porifera, as well as to place the various lower metazoan groups on the animal-kingdom tree. Most recently, the partial mtDNA of two glass sponges, class Hexactinellida, were reported. While previous phylogenetic estimations based on these data remain uncertain due to insufficient taxon sampling and accelerated rates of evolution, the mtDNA molecules themselves reveal interesting traits that may be unique to hexactinellids. Here we determined the first complete mitochondrial genome of a hexactinellid sponge, Aphrocallistes vastus, and compared it to published poriferan mtDNAs to further describe characteristics specific to hexactinellid and other sponge mitochondrial genomes. RESULTS: The A. vastus mtDNA consisted of a 17,427 base pair circular molecule containing thirteen protein-coding genes, divergent large and small subunit ribosomal RNAs, and a reduced set of 18 tRNAs. The A. vastus mtDNA showed a typical hexactinellid nucleotide composition and shared a large synteny with the other sequenced glass sponge mtDNAs. It also contained an unidentified open reading frame and large intergenic space region. Two frameshifts, in the cox3 and nad6 genes, were not corrected by RNA editing, but rather possessed identical shift sites marked by the extremely rare tryptophan codon (UGG) followed by the common glycine codon (GGA) in the +1 frame. CONCLUSION: Hexactinellid mtDNAs have shown similar trends in gene content, nucleotide composition, and codon usage, and have retained a large gene syntenty. Analysis of the mtDNA of A. vastus has provided evidence diagnostic for +1 programmed translational frameshifting, a phenomenon disparately reported throughout the animal kingdom, but present in the hexactinellid mtDNAs that have been sequenced to date.

A guard-cell-specific MYB transcription factor regulates stomatal movements and plant drought tolerance.

Stomatal pores located on the plant epidermis regulate CO(2) uptake for photosynthesis and the loss of water by transpiration. The opening and closing of the pore is mediated by turgor-driven volume changes of two surrounding guard cells. These highly specialized cells integrate internal signals and environmental stimuli to modulate stomatal aperture for plant survival under diverse conditions. Modulation of transcription and mRNA processing play important roles in controlling guard-cell activity, although the details of these levels of regulation remain mostly unknown. Here we report the characterization of AtMYB60, a R2R3-MYB gene of Arabidopsis, as the first transcription factor involved in the regulation of stomatal movements. AtMYB60 is specifically expressed in guard cells, and its expression is negatively modulated during drought. A null mutation in AtMYB60 results in the constitutive reduction of stomatal opening and in decreased wilting under water stress conditions. Transcript levels of a limited number of genes are altered in the mutant, and many of these genes are involved in the plant response to stress. Our data indicate that AtMYB60 is a transcriptional modulator of physiological responses in guard cells and open new possibilities to engineering stomatal activity to help plants survive desiccation.

Differential effect of allorecognition loci on phenotype in Hydractinia symbiolongicarpus (Cnidaria: Hydrozoa).

The allorecognition complex of Hydractinia symbiolongicarpus is a chromosomal interval containing two loci, alr1 and alr2, that controls fusion between genetically distinct colonies. Recombination between these two loci has been associated with a heterogeneous class of phenotypes called transitory fusion. A large-scale backcross was performed to generate a population of colonies (N = 106) with recombination breakpoints within the allorecognition complex. Two distinct forms of transitory fusion were correlated with reciprocal recombination products, suggesting that alr1 and alr2 contributed differentially to the allorecognition response. Specifically, type I transitory fusion is associated with rapid and persistent separation of allogeneic tissues, whereas type II transitory fusion generates a patchwork of continuously fusing and separating tissues.

Cell cycle arrest of stamen initials in maize sex determination.

The maize sex determination pathway results in the arrest of stamen in ear spikelets and the abortion of pistils in both the tassel spikelets and in the secondary florets of ear spikelets. Arrested stamen cells showed no signs of DNA fragmentation, an absence of CYCLIN B expression, and an accumulation of the negative cell cycle regulator WEE1 RNA.

Biochemical characterization of TASSELSEED 2, an essential plant short-chain dehydrogenase/reductase with broad spectrum activities.

The development of unisexual flowers in maize and other plants proceeds through selective elimination of floral organs in an initially bisexual floral meristem. The essential character of the tasselseed 2 gene (TS2) in this cell-death pathway has been established previously. Molecular cloning of TS2 reveals membership to the evolutionarily conserved superfamily of short-chain dehydrogenases/reductases, but its substrate specificity remained unknown. Recombinant TS2 protein was produced in Escherichia coli, and purified to apparent homogeneity. Analytical ultracentrifugation and gel filtration experiments show that TS2 is a tetrameric enzyme. Thermal denaturation followed by circular dichroism spectroscopy reveals that TS2 binds NAD(H) and NAD(P)(H). Substrate screening demonstrates that TS2 converts steroids with specificities found at positions 3 and 17, and several dicarbonyl and quinone compounds, thus establishing TS2 as a plant 3beta/17beta-hydroxysteroid dehydrogenase and carbonyl/quinone reductase. Taken together, the genetic data and the substrate specificities determined suggest that TS2 converts specific plant compounds and acts as a prereceptor control mechanism, in a manner similar to that of mammalian hydroxysteroid dehydrogenases.

Genetic Architecture of a Rice Nested Association Mapping Population.

Describing the genetic diversity in the gene pool of crops will provide breeders with novel resources for varietal improvement. Nested Association Mapping (NAM) populations are uniquely suited for characterizing parental diversity through the shuffling and fixation of parental haplotypes. Here, we describe a set of 1879 rice NAM lines created through the selfing and single-seed descent of F1 hybrids derived from elite IR64 indica crossed with 10 diverse tropical japonica lines. Genotyping data indicated tropical japonica alleles were captured at every queried locus despite the presence of segregation distortion factors. Several distortion loci were mapped, both shared and unique, among the 10 populations. Using two-point and multi-point genetic map calculations, our datasets achieved the ∼1500 cM expected map size in rice. Finally, we highlighted the utility of the NAM lines for QTL mapping, including joint analysis across the 10 populations, by confirming known QTL locations for the trait days to heading.

Caribbean placozoan phylogeography.

We here address placozoan distribution and phylogeography in five locations in the Caribbean Sea. We performed a coarse-resolution presence/absence survey of placozoans in Belize, Bermuda, Grenada, Jamaica, and Panama and a fine-resolution study of the distribution of placozoans in Twin Cays, Belize. Placozoans were recovered in every country sampled. Animals were sequenced at the mitochondrial 16S rDNA locus, and our analysis identified four of the five previously identified clades present in the Caribbean. In addition, we discovered two new haplotypes within one of these clades, and we found sympatric clades in Belize, Bermuda, Jamaica, and Panama. These studies provide further molecular evidence for species diversity within the Phylum Placozoa.

Imputing Genotypes in Biallelic Populations from Low-Coverage Sequence Data.

Low-coverage next-generation sequencing methodologies are routinely employed to genotype large populations. Missing data in these populations manifest both as missing markers and markers with incomplete allele recovery. False homozygous calls at heterozygous sites resulting from incomplete allele recovery confound many existing imputation algorithms. These types of systematic errors can be minimized by incorporating depth-of-sequencing read coverage into the imputation algorithm. Accordingly, we developed Low-Coverage Biallelic Impute (LB-Impute) to resolve missing data issues. LB-Impute uses a hidden Markov model that incorporates marker read coverage to determine variable emission probabilities. Robust, highly accurate imputation results were reliably obtained with LB-Impute, even at extremely low (<1×) average per-marker coverage. This finding will have implications for the design of genotype imputation algorithms in the future. LB-Impute is publicly available on GitHub at https://github.com/dellaporta-laboratory/LB-Impute.

Control of sexuality by the sk1-encoded UDP-glycosyltransferase of maize.

Sex determination in maize involves the production of staminate and pistillate florets from an initially bisexual floral meristem. Pistil elimination in staminate florets requires jasmonic acid signaling, and functional pistils are protected by the action of the silkless 1 (sk1) gene. The sk1 gene was identified and found to encode a previously uncharacterized family 1 uridine diphosphate glycosyltransferase that localized to the plant peroxisomes. Constitutive expression of an sk1 transgene protected all pistils in the plant, causing complete feminization, a gain-of-function phenotype that operates by blocking the accumulation of jasmonates. The segregation of an sk1 transgene was used to effectively control the production of pistillate and staminate inflorescences in maize plants.