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.

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.

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.

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.

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.

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.

Genomic Characterization of Interspecific Hybrids and an Admixture Population Derived from Panicum amarum × P. virgatum.

Switchgrass (Panicum virgatum L.) and its relatives are regarded as top bioenergy crop candidates; however, one critical barrier is the introduction of useful genetic diversity and the development of new cultivars and hybrids. Combining genomes from related cultivars and species provides an opportunity to introduce new traits. In switchgrass, a breeding advantage would be achieved by combining the genomes of intervarietal ecotypes or interspecific hybrids. The recovery of wide crosses, however, is often tedious and may involve complicated embryo rescue and numerous backcrosses. Here, we demonstrate a straightforward approach to wide crosses involving the use of a selectable transgene for recovery of interspecific [P. virgatum cv. Alamo × Panicum amarum Ell. var amarulum or Atlantic Coastal Panicgrass (ACP)] F1 hybrids followed by backcrossing to generate a nontransgenic admixture population. A nontransgenic herbicide-sensitive (HbS) admixture population of 83 F1 BC1 progeny was analyzed by genotyping-by-sequencing (GBS) to characterize local ancestry, parental contribution, and patterns of recombination. These results demonstrate a widely applicable breeding strategy that makes use of transgenic selectable resistance to identify and recover true hybrids.

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.

Identification of the maize gravitropism gene lazy plant1 by a transposon-tagging genome resequencing strategy.

Since their initial discovery, transposons have been widely used as mutagens for forward and reverse genetic screens in a range of organisms. The problems of high copy number and sequence divergence among related transposons have often limited the efficiency at which tagged genes can be identified. A method was developed to identity the locations of Mutator (Mu) transposons in the Zea mays genome using a simple enrichment method combined with genome resequencing to identify transposon junction fragments. The sequencing library was prepared from genomic DNA by digesting with a restriction enzyme that cuts within a perfectly conserved motif of the Mu terminal inverted repeats (TIR). Paired-end reads containing Mu TIR sequences were computationally identified and chromosomal sequences flanking the transposon were mapped to the maize reference genome. This method has been used to identify Mu insertions in a number of alleles and to isolate the previously unidentified lazy plant1 (la1) gene. The la1 gene is required for the negatively gravitropic response of shoots and mutant plants lack the ability to sense gravity. Using bioinformatic and fluorescence microscopy approaches, we show that the la1 gene encodes a cell membrane and nuclear localized protein. Our Mu-Taq method is readily adaptable to identify the genomic locations of any insertion of a known sequence in any organism using any sequencing platform.

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.

Genetic Background and Allorecognition Phenotype in Hydractinia symbiolongicarpus.

The Hydractinia allorecognition complex (ARC) was initially identified as a single chromosomal interval using inbred and congenic lines. The production of defined lines necessarily homogenizes genetic background and thus may be expected to obscure the effects of unlinked allorecognition loci should they exist. Here, we report the results of crosses in which inbred lines were out-crossed to wild-type animals in an attempt to identify dominant, codominant, or incompletely dominant modifiers of allorecognition. A claim for the existence of modifiers unlinked to ARC was rejected for three different genetic backgrounds. Estimates of the genetic map distance of ARC in two wild-type haplotypes differed markedly from one another and from that measured in congenic lines. These results suggest that additional allodeterminants exist in the Hydractinia ARC.

Hydractinia allodeterminant alr1 resides in an immunoglobulin superfamily-like gene complex.

Allorecognition, the ability to discriminate between self and nonself, is ubiquitous among colonial metazoans and widespread among aclonal taxa. Genetic models for the study of allorecognition have been developed in the jawed vertebrates, invertebrate chordate Botryllus, and cnidarian Hydractinia. In Botryllus, two genes contribute to the histocompatibility response, FuHC and fester. In the cnidarian Hydractinia, one of the two known allorecognition loci, alr2, has been isolated, and a second linked locus, alr1, has been mapped to the same chromosomal region, called the allorecognition complex (ARC). Here we isolate alr1 by positional cloning and report it to encode a transmembrane receptor protein with two hypervariable extracellular regions similar to immunoglobulin (Ig)-like domains. Variation in the extracellular domain largely predicts fusibility within and between laboratory strains and wild-type isolates. alr1 was found embedded in a family of immunoglobulin superfamily (IgSF)-like genes, thus establishing that the ARC histocompatibility complex is an invertebrate IgSF-like gene complex.

A hypervariable invertebrate allodeterminant.

Colonial marine invertebrates, such as sponges, corals, bryozoans, and ascidians, often live in densely populated communities where they encounter other members of their species as they grow over their substratum. Such encounters typically lead to a natural histocompatibility response in which colonies either fuse to become a single, chimeric colony or reject and aggressively compete for space. These allorecognition phenomena mediate intraspecific competition, support allotypic diversity, control the level at which selection acts, and resemble allogeneic interactions in pregnancy and transplantation. Despite the ubiquity of allorecognition in colonial phyla, however, its molecular basis has not been identified beyond what is currently known about histocompatibility in vertebrates and protochordates. We positionally cloned an allorecognition gene by using inbred strains of the cnidarian, Hydractinia symbiolongicarpus, which is a model system for the study of invertebrate allorecognition. The gene identified encodes a putative transmembrane receptor expressed in all tissues capable of allorecognition that is highly polymorphic and predicts allorecognition responses in laboratory and field-derived strains. This study reveals that a previously undescribed hypervariable molecule bearing three extracellular domains with greatest sequence similarity to the immunoglobulin superfamily is an allodeterminant in a lower metazoan.

tasselseed1 is a lipoxygenase affecting jasmonic acid signaling in sex determination of maize.

Sex determination in maize is controlled by a developmental cascade leading to the formation of unisexual florets derived from an initially bisexual floral meristem. Abortion of pistil primordia in staminate florets is controlled by a tasselseed-mediated cell death process. We positionally cloned and characterized the function of the sex determination gene tasselseed1 (ts1). The TS1 protein encodes a plastid-targeted lipoxygenase with predicted 13-lipoxygenase specificity, which suggests that TS1 may be involved in the biosynthesis of the plant hormone jasmonic acid. In the absence of a functional ts1 gene, lipoxygenase activity was missing and endogenous jasmonic acid concentrations were reduced in developing inflorescences. Application of jasmonic acid to developing inflorescences rescued stamen development in mutant ts1 and ts2 inflorescences, revealing a role for jasmonic acid in male flower development in maize.

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.

The SDR (short-chain dehydrogenase/reductase and related enzymes) nomenclature initiative.

Short-chain dehydrogenases/reductases (SDR) constitute one of the largest enzyme superfamilies with presently over 46,000 members. In phylogenetic comparisons, members of this superfamily show early divergence where the majority have only low pairwise sequence identity, although sharing common structural properties. The SDR enzymes are present in virtually all genomes investigated, and in humans over 70 SDR genes have been identified. In humans, these enzymes are involved in the metabolism of a large variety of compounds, including steroid hormones, prostaglandins, retinoids, lipids and xenobiotics. It is now clear that SDRs represent one of the oldest protein families and contribute to essential functions and interactions of all forms of life. As this field continues to grow rapidly, a systematic nomenclature is essential for future annotation and reference purposes. A functional subdivision of the SDR superfamily into at least 200 SDR families based upon hidden Markov models forms a suitable foundation for such a nomenclature system, which we present in this paper using human SDRs as examples.

The early ANTP gene repertoire: insights from the placozoan genome.

The evolution of ANTP genes in the Metazoa has been the subject of conflicting hypotheses derived from full or partial gene sequences and genomic organization in higher animals. Whole genome sequences have recently filled in some crucial gaps for the basal metazoan phyla Cnidaria and Porifera. Here we analyze the complete genome of Trichoplax adhaerens, representing the basal metazoan phylum Placozoa, for its set of ANTP class genes. The Trichoplax genome encodes representatives of Hox/ParaHox-like, NKL, and extended Hox genes. This repertoire possibly mirrors the condition of a hypothetical cnidarian-bilaterian ancestor. The evolution of the cnidarian and bilaterian ANTP gene repertoires can be deduced by a limited number of cis-duplications of NKL and "extended Hox" genes and the presence of a single ancestral "ProtoHox" gene.

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.