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

\u3cem\u3eIn situ\u3c/em\u3e embryo rescue for generation of wide intra- and interspecific hybrids of \u3cem\u3ePanicum virgatum\u3c/em\u3e 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-bysequencing 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

Genomic Characterization of Interspecific Hybrids and an Admixture Population Derived from Panicum amarum × \u3cem\u3eP. virgatum\u3c/em\u3e

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 F1BC1 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

Control of sexuality by the \u3cem\u3esk1\u3c/em\u3e-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

Constructing linkage maps in the genomics era with MapDisto 2.0

International audienceMotivation: Genotyping by sequencing (GBS) generates datasets that are challenging to handle by current genetic mapping software with graphical interface. Geneticists need new user-friendly computer programs that can analyze GBS data on desktop computers. This requires improvements in computation efficiency, both in terms of speed and use of random-access memory (RAM).Results: MapDisto v.2.0 is a user-friendly computer program for construction of genetic linkage maps. It includes several new major features: (i) handling of very large genotyping datasets like the ones generated by GBS; (ii) direct importation and conversion of Variant Call Format (VCF) files; (iii) detection of linkage, i.e. construction of linkage groups in case of segregation distortion; (iv) data imputation on VCF files using a new approach, called LB-Impute. Features i to iv operate through inclusion of new Java modules that are used transparently by MapDisto; (v) QTL detection via a new R/qtl graphical interface

A high throughput embryonic stem cell screen identifies Oct-2 as a bifunctional regulator of neuronal differentiation

Neuronal differentiation is a complex process that involves a plethora of regulatory steps. To identify transcription factors that influence neuronal differentiation we developed a high throughput screen using embryonic stem (ES) cells. Seven-hundred human transcription factor clones were stably introduced into mouse ES (mES) cells and screened for their ability to induce neuronal differentiation of mES cells. Twenty-four factors that are capable of inducing neuronal differentiation were identified, including four known effectors of neuronal differentiation, 11 factors with limited evidence of involvement in regulating neuronal differentiation, and nine novel factors. One transcription factor, Oct-2, was studied in detail and found to be a bifunctional regulator: It can either repress or induce neuronal differentiation, depending on the particular isoform. Ectopic expression experiments demonstrate that isoform Oct-2.4 represses neuronal differentiation, whereas Oct-2.2 activates neuron formation. Consistent with a role in neuronal differentiation, Oct-2.2 expression is induced during differentiation, and cells depleted of Oct-2 and its homolog Oct-1 have a reduced capacity to differentiate into neurons. Our results reveal a number of transcription factors potentially important for mammalian neuronal differentiation, and indicate that Oct-2 may serve as a binary switch to repress differentiation in precursor cells and induce neuronal differentiation later during neuronal development

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

A parthenogenesis gene candidate and evidence for segmental allopolyploidy in apomictic brachiaria decumbens

Apomixis, asexual reproduction through seed, enables breeders to identify and faithfully propagate superior heterozygous genotypes by seed without the disadvantages of vegetative propagation or the expense and complexity of hybrid seed production. The availability of new tools such as genotyping by sequencing and bioinformatics pipelines for species lacking reference genomes now makes the construction of dense maps possible in apomictic species, despite complications including polyploidy, multisomic inheritance, self-incompatibility, and high levels of heterozygosity. In this study, we developed saturated linkage maps for the maternal and paternal genomes of an interspecific Brachiaria ruziziensis (R. Germ. and C. M. Evrard) × B. decumbens Stapf. F1 mapping population in order to identify markers linked to apomixis. High-resolution molecular karyotyping and comparative genomics with Setaria italica (L.) P. Beauv provided conclusive evidence for segmental allopolyploidy in B. decumbens, with strong preferential pairing of homologs across the genome and multisomic segregation relatively more common in chromosome 8. The apospory-specific genomic region (ASGR) was mapped to a region of reduced recombination on B. decumbens chromosome 5. The Pennisetum squamulatum (L.) R.Br. PsASGR-BABY BOOM-like (psASGR–BBML)-specific primer pair p779/p780 was in perfect linkage with the ASGR in the F1 mapping population and diagnostic for reproductive mode in a diversity panel of known sexual and apomict Brachiaria (Trin.) Griseb. and P. maximum Jacq. germplasm accessions and cultivars. These findings indicate that ASGR–BBML gene sequences are highly conserved across the Paniceae and add further support for the postulation of the ASGR–BBML as candidate genes for the apomictic function of parthenogenesis

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

Correlation of Global MicroRNA Expression With Basal Cell Carcinoma Subtype

Basal cell carcinomas (BCCs) are the most common cancers in the United States. The histologic appearance distinguishes several subtypes, each of which can have a different biologic behavior. In this study, global miRNA expression was quantified by high-throughput sequencing in nodular BCCs, a subtype that is slow growing, and infiltrative BCCs, aggressive tumors that extend through the dermis and invade structures such as cutaneous nerves. Principal components analysis correctly classified seven of eight infiltrative tumors on the basis of miRNA expression. The remaining tumor, on pathology review, contained a mixture of nodular and infiltrative elements. Nodular tumors did not cluster tightly, likely reflecting broader histopathologic diversity in this class, but trended toward forming a group separate from infiltrative BCCs. Quantitative polymerase chain reaction assays were developed for six of the miRNAs that showed significant differences between the BCC subtypes, and five of these six were validated in a replication set of four infiltrative and three nodular tumors. The expression level of miR-183, a miRNA that inhibits invasion and metastasis in several types of malignancies, was consistently lower in infiltrative than nodular tumors and could be one element underlying the difference in invasiveness. These results represent the first miRNA profiling study in BCCs and demonstrate that miRNA gene expression may be involved in tumor pathogenesis and particularly in determining the aggressiveness of these malignancies