Hypothalamic Transcriptome Response To Simulated Diel Earthen Pond Hypoxia Cycles In Channel Catfish (Ictalurus punctatus).

Commercial culture of channel catfish (Ictalurus punctatus) occurs in earthen ponds that are characterized by diel swings in dissolved oxygen concentration that can fall to severe levels of hypoxia which can suppress appetite and lead to suboptimal growth. Given the significance of the hypothalamus in regulating these processes in other fishes, an investigation into the hypothalamus transcriptome was conducted to identify specific genes and expression patterns responding to hypoxia. Channel catfish in normoxic water were compared to catfish subjected to 12 hours of hypoxia (20% oxygen saturation; 1.8 mg O2/L; 27 °C) followed by 12 hours of recovery in normoxia to mimic 24-hours in a catfish aquaculture pond. Fish were sampled at 0-, 6-, 12-, 18-, and 24-hour time points, with the 6- and 12-hour samplings occurring during hypoxia. A total of 190 genes were differentially expressed during the experiment, with most occurring during hypoxia and returning to baseline values within 6 hours of normoxia. Differentially expressed genes were sorted by function into Gene Ontology biological processes and revealed that most were categorized as "response to hypoxia", "sprouting angiogenesis", and "cellular response to xenobiotic stimulus". The patterns of gene expression reported here suggest that transcriptome responses to hypoxia are broad and quickly reversibly with the onset of normoxia. Although no genes commonly reported to modulate appetite were found to be differentially expressed in this experiment, several candidates were identified for future studies investigating the interplay between hypoxia and appetite in channel catfish, including adm, igfbp1a, igfbp7, and stc2b.

Gene disruption by structural mutations drives selection in US rice breeding over the last century.

The genetic basis of general plant vigor is of major interest to food producers, yet the trait is recalcitrant to genetic mapping because of the number of loci involved, their small effects, and linkage. Observations of heterosis in many crops suggests that recessive, malfunctioning versions of genes are a major cause of poor performance, yet we have little information on the mutational spectrum underlying these disruptions. To address this question, we generated a long-read assembly of a tropical japonica rice (Oryza sativa) variety, Carolina Gold, which allowed us to identify structural mutations (>50 bp) and orient them with respect to their ancestral state using the outgroup, Oryza glaberrima. Supporting prior work, we find substantial genome expansion in the sativa branch. While transposable elements (TEs) account for the largest share of size variation, the majority of events are not directly TE-mediated. Tandem duplications are the most common source of insertions and are highly enriched among 50-200bp mutations. To explore the relative impact of various mutational classes on crop fitness, we then track these structural events over the last century of US rice improvement using 101 resequenced varieties. Within this material, a pattern of temporary hybridization between medium and long-grain varieties was followed by recent divergence. During this long-term selection, structural mutations that impact gene exons have been removed at a greater rate than intronic indels and single-nucleotide mutations. These results support the use of ab initio estimates of mutational burden, based on structural data, as an orthogonal predictor in genomic selection.

Legacy genetics of Arachis cardenasii in the peanut crop shows the profound benefits of international seed exchange.

The narrow genetics of most crops is a fundamental vulnerability to food security. This makes wild crop relatives a strategic resource of genetic diversity that can be used for crop improvement and adaptation to new agricultural challenges. Here, we uncover the contribution of one wild species accession, Arachis cardenasii GKP 10017, to the peanut crop (Arachis hypogaea) that was initiated by complex hybridizations in the 1960s and propagated by international seed exchange. However, until this study, the global scale of the dispersal of genetic contributions from this wild accession had been obscured by the multiple germplasm transfers, breeding cycles, and unrecorded genetic mixing between lineages that had occurred over the years. By genetic analysis and pedigree research, we identified A. cardenasii-enhanced, disease-resistant cultivars in Africa, Asia, Oceania, and the Americas. These cultivars provide widespread improved food security and environmental and economic benefits. This study emphasizes the importance of wild species and collaborative networks of international expertise for crop improvement. However, it also highlights the consequences of the implementation of a patchwork of restrictive national laws and sea changes in attitudes regarding germplasm that followed in the wake of the Convention on Biological Diversity. Today, the botanical collections and multiple seed exchanges which enable benefits such as those revealed by this study are drastically reduced. The research reported here underscores the vital importance of ready access to germplasm in ensuring long-term world food security.

Long-read, chromosome-scale assembly of Vitis rotundifolia cv. Carlos and its unique resistance to Xylella fastidiosa subsp. fastidiosa.

BACKGROUND: Muscadine grape (Vitis rotundifolia) is resistant to many of the pathogens that negatively impact the production of common grape (V. vinifera), including the bacterial pathogen Xylella fastidiosa subsp. fastidiosa (Xfsf), which causes Pierce's Disease (PD). Previous studies in common grape have indicated Xfsf delays host immune response with a complex O-chain antigen produced by the wzy gene. Muscadine cultivars range from tolerant to completely resistant to Xfsf, but the mechanism is unknown. RESULTS: We assembled and annotated a new, long-read genome assembly for 'Carlos', a cultivar of muscadine that exhibits tolerance, to build upon the existing genetic resources available for muscadine. We used these resources to construct an initial pan-genome for three cultivars of muscadine and one cultivar of common grape. This pan-genome contains a total of 34,970 synteny-constrained entries containing genes of similar structure. Comparison of resistance gene content between the 'Carlos' and common grape genomes indicates an expansion of resistance (R) genes in 'Carlos.' We further identified genes involved in Xfsf response by transcriptome sequencing 'Carlos' plants inoculated with Xfsf. We observed 234 differentially expressed genes with functions related to lipid catabolism, oxidation-reduction signaling, and abscisic acid (ABA) signaling as well as seven R genes. Leveraging public data from previous experiments of common grape inoculated with Xfsf, we determined that most differentially expressed genes in the muscadine response were not found in common grape, and three of the R genes identified as differentially expressed in muscadine do not have an ortholog in the common grape genome. CONCLUSIONS: Our results support the utility of a pan-genome approach to identify candidate genes for traits of interest, particularly disease resistance to Xfsf, within and between muscadine and common grape.

Microscopic and Transcriptomic Analyses of Dalbergoid Legume Peanut Reveal a Divergent Evolution Leading to Nod-Factor-Dependent Epidermal Crack-Entry and Terminal Bacteroid Differentiation.

Root nodule symbiosis (RNS) is the pillar behind sustainable agriculture and plays a pivotal role in the environmental nitrogen cycle. Most of the genetic, molecular, and cell-biological knowledge on RNS comes from model legumes that exhibit a root-hair mode of bacterial infection, in contrast to the Dalbergoid legumes exhibiting crack-entry of rhizobia. As a step toward understanding this important group of legumes, we have combined microscopic analysis and temporal transcriptome to obtain a dynamic view of plant gene expression during Arachis hypogaea (peanut) nodule development. We generated comprehensive transcriptome data by mapping the reads to A. hypogaea, and two diploid progenitor genomes. Additionally, we performed BLAST searches to identify nodule-induced yet-to-be annotated peanut genes. Comparison between peanut, Medicago truncatula, Lotus japonicus, and Glycine max showed upregulation of 61 peanut orthologs among 111 tested known RNS-related genes, indicating conservation in mechanisms of nodule development among members of the Papilionoid family. Unlike model legumes, recruitment of class 1 phytoglobin-derived symbiotic hemoglobin (SymH) in peanut indicates diversification of oxygen-scavenging mechanisms in the Papilionoid family. Finally, the absence of cysteine-rich motif-1-containing nodule-specific cysteine-rich peptide (NCR) genes but the recruitment of defensin-like NCRs suggest a diverse molecular mechanism of terminal bacteroid differentiation. In summary, our work describes genetic conservation and diversification in legume-rhizobia symbiosis in the Papilionoid family, as well as among members of the Dalbergoid legumes.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Pyrrocidine, a molecular off switch for fumonisin biosynthesis.

Sarocladium zeae is a fungal endophyte of maize and can be found co-inhabiting a single seed with Fusarium verticillioides, a major mycotoxigenic food safety threat. S. zeae produces pyrrocidines A and B that inhibit the growth of F. verticillioides and may limit its spread within the seed to locations lacking S. zeae. Although coinhabiting single seeds, the fungi are generally segregated in separate tissues. To understand F. verticillioides' protective physiological response to pyrrocidines we sequenced the F. verticillioides transcriptome upon exposure to purified pyrrocidine A or B at sub-inhibitory concentrations. Through this work we identified a F. verticillioides locus FvABC3 (FVEG_11089) encoding a transporter critical for resistance to pyrrocidine. We also identified FvZBD1 (FVEG_00314), a gene directly adjacent to the fumonisin biosynthetic gene cluster that was induced several thousand-fold in response to pyrrocidines. FvZBD1 is postulated to act as a genetic repressor of fumonisin production since deletion of the gene resulted in orders of magnitude increase in fumonisin. Further, pyrrocidine acts, likely through FvZBD1, to shut off fumonisin biosynthesis. This suggests that S. zeae is able to hack the secondary metabolic program of a competitor fungus, perhaps as preemptive self-protection, in this case impacting a mycotoxin of central concern for food safety.

Molecular insight into cotton leaf curl geminivirus disease resistance in cultivated cotton (Gossypium hirsutum).

Cultivated cotton (Gossypium hirsutum) is the most important fibre crop in the world. Cotton leaf curl disease (CLCuD) is the major limiting factor and a threat to textile industry in India and Pakistan. All the local cotton cultivars exhibit moderate to no resistance against CLCuD. In this study, we evaluated an exotic cotton accession Mac7 as a resistance source to CLCuD by challenging it with viruliferous whiteflies and performing qPCR to evaluate the presence/absence and relative titre of CLCuD-associated geminiviruses/betasatellites. The results indicated that replication of pathogenicity determinant betasatellite is significantly attenuated in Mac7 and probably responsible for resistance phenotype. Afterwards, to decipher the genetic basis of CLCuD resistance in Mac7, we performed RNA sequencing on CLCuD-infested Mac7 and validated RNA-Seq data with qPCR on 24 independent genes. We performed co-expression network and pathway analysis for regulation of geminivirus/betasatellite-interacting genes. We identified nine novel modules with 52 hubs of highly connected genes in network topology within the co-expression network. Analysis of these hubs indicated the differential regulation of auxin stimulus and cellular localization pathways in response to CLCuD. We also analysed the differential regulation of geminivirus/betasatellite-interacting genes in Mac7. We further performed the functional validation of selected candidate genes via virus-induced gene silencing (VIGS). Finally, we evaluated the genomic context of resistance responsive genes and found that these genes are not specific to A or D sub-genomes of G. hirsutum. These results have important implications in understanding CLCuD resistance mechanism and developing a durable resistance in cultivated cotton.

Whole genome sequencing of Asia II 1 species of whitefly reveals that genes involved in virus transmission and insecticide resistance have genetic variances between Asia II 1 and MEAM1 species.

BACKGROUND: Whiteflies (Bemisia tabaci) are phloem sap-sucking pests that because of their broad host range and ability to transmit viruses damage crop plants worldwide. B. tabaci are now known to be a complex of cryptic species that differ from each other in many characteristics such as mode of interaction with viruses, invasiveness, and resistance to insecticides. Asia II 1 is an indigenous species found on the Indian sub-continent and south-east Asia while the species named as Middle East Asia Minor 1 (MEAM1), likely originated from the Middle-East and has spread worldwide in recent decades. The purpose of this study is to find genomic differences between these two species. RESULTS: Sequencing of the nuclear genome of Asia II 1 with Illumina HiSeq and MiSeq generated 198.90 million reads that covers 88% of the reference genome. The sequence comparison with MEAM1 identified 2,327,972 SNPs and 202,479 INDELs. In Total, 1294 genes were detected with high impact variants. The functional analysis revealed that some of the genes are involved in virus transmission including 4 genes in Tomato yellow leaf curl virus (TYLCV) transmission, 96 in Tomato crinivirus (ToCV) transmission, and 14 genes in insecticide resistance. CONCLUSIONS: These genetic differences between Asia II 1 and MEAM1 may underlie the major biological differences between the two species such as virus transmission, insecticide resistance, and range of host plants. The present study provides new genomic data and information resources for Asia II 1 that will not only contribute to the species delimitation of whitefly, but also help in conceiving future research studies to develop more targeted management strategies against whitefly.

New tools to screen wild peanut species for aflatoxin accumulation and genetic fingerprinting.

BACKGROUND: Aflatoxin contamination in peanut seeds is still a serious problem for the industry and human health. No stable aflatoxin resistant cultivars have yet been produced, and given the narrow genetic background of cultivated peanuts, wild species became an important source of genetic diversity. Wild peanut seeds, however, are not abundant, thus, an effective method of screening for aflatoxin accumulation using minimal seeds is highly desirable. In addition, keeping record of genetic fingerprinting of each accession would be very useful for breeding programs and for the identification of accessions within germplasm collections. RESULTS: In this study, we report a method of screening for aflatoxin accumulation that is applicable to the small-size seeds of wild peanuts, increases the reliability by testing seed viability, and records the genetic fingerprinting of the samples. Aflatoxin levels observed among 20 wild peanut species varied from zero to 19000 ng.g-1 and 155 ng.g-1 of aflatoxin B1 and B2, respectively. We report the screening of 373 molecular markers, including 288 novel SSRs, tested on 20 wild peanut species. Multivariate analysis by Neighbor-Joining, Principal Component Analysis and 3D-Principal Coordinate Analysis using 134 (36 %) transferable markers, in general grouped the samples according to their reported genomes. The best 88 markers, those with high fluorescence, good scorability and transferability, are reported with BLAST results. High quality markers (total 98) that discriminated genomes are reported. A high quality marker with UPIC score 16 (16 out of 20 species discriminated) had significant hits on BLAST2GO to a pentatricopeptide-repeat protein, another marker with score 5 had hits on UDP-D-apiose synthase, and a third one with score 12 had BLASTn hits on La-RP 1B protein. Together, these three markers discriminated all 20 species tested. CONCLUSIONS: This study provides a reliable method to screen wild species of peanut for aflatoxin resistance using minimal seeds. In addition we report 288 new SSRs for peanut, and a cost-effective combination of markers sufficient to discriminate all 20 species tested. These tools can be used for the systematic search of aflatoxin resistant germplasm keeping record of the genetic fingerprinting of the accessions tested for breeding purpose.

Study of the genetic diversity of the aflatoxin biosynthesis cluster in Aspergillus section Flavi using insertion/deletion markers in peanut seeds from Georgia, USA.

Aflatoxins are among the most powerful carcinogens in nature. The major aflatoxin-producing fungi are Aspergillus flavus and A. parasiticus. Numerous crops, including peanut, are susceptible to aflatoxin contamination by these fungi. There has been an increased use of RNA interference (RNAi) technology to control phytopathogenic fungi in recent years. In order to develop molecular tools targeting specific genes of these fungi for the control of aflatoxins, it is necessary to obtain their genome sequences. Although high-throughput sequencing is readily available, it is still impractical to sequence the genome of every isolate. Thus, in this work, the authors proposed a workflow that allowed prescreening of 238 Aspergillus section Flavi isolates from peanut seeds from Georgia, USA. The aflatoxin biosynthesis cluster (ABC) of the isolates was fingerprinted at 25 InDel (insertion/deletion) loci using capillary electrophoresis. All isolates were tested for aflatoxins using ultra-high-performance liquid chromatography. The neighbor-joining, three-dimension (3D) principal coordinate, and Structure analyses revealed that the Aspergillus isolates sampled consisted of three main groups determined by their capability to produce aflatoxins. Group I comprised 10 non-aflatoxigenic A. flavus; Group II included A. parasiticus; and Group III included mostly aflatoxigenic A. flavus and the three non-aflatoxigenic A. caelatus. Whole genomes of 10 representative isolates from different groups were sequenced. Although InDels in Aspergillus have been used by other research groups, this is the first time that the cluster analysis resulting from fingerprinting was followed by whole-genome sequencing of representative isolates. In our study, cluster analysis of ABC sequences validated the results obtained with fingerprinting. This shows that InDels used here can predict similarities at the genome level. Our results also revealed a relationship between groups and their capability to produce aflatoxins. The database generated of Aspergillus spp. can be used to select target genes and assess the effectiveness of RNAi technology to reduce aflatoxin contamination in peanut.

Genetic Diversity of the Two Commercial Tetraploid Cotton Species in the Gossypium Diversity Reference Set.

A diversity reference set has been constructed for the Gossypium accessions in the US National Cotton Germplasm Collection to facilitate more extensive evaluation and utilization of accessions held in the Collection. A set of 105 mapped simple sequence repeat markers was used to study the allelic diversity of 1933 tetraploid Gossypium accessions representative of the range of diversity of the improved and wild accessions of G. hirsutum and G. barbadense. The reference set contained 410 G. barbadense accessions and 1523 G. hirsutum accessions. Observed numbers of polymorphic and private bands indicated a greater diversity in G. hirsutum as compared to G. barbadense as well as in wild-type accessions as compared to improved accessions in both species. The markers clearly differentiated the 2 species. Patterns of diversity within species were observed but not clearly delineated, with much overlap occurring between races and regions of origin for wild accessions and between historical and geographic breeding pools for cultivated accessions. Although the percentage of accessions showing introgression was higher among wild accessions than cultivars in both species, the average level of introgression within individual accessions, as indicated by species-specific bands, was much higher in wild accessions of G. hirsutum than in wild accessions of G. barbadense. The average level of introgression within individual accessions was higher in improved G. barbadense cultivars than in G. hirsutum cultivars. This molecular characterization reveals the levels and distributions of genetic diversity that will allow for better exploration and utilization of cotton genetic resources.

Metabolomic and transcriptomic insights into how cotton fiber transitions to secondary wall synthesis, represses lignification, and prolongs elongation.

BACKGROUND: The morphogenesis of single-celled cotton fiber includes extreme elongation and staged cell wall differentiation. Designing strategies for improving cotton fiber for textiles and other uses relies on uncovering the related regulatory mechanisms. In this research we compared the transcriptomes and metabolomes of two Gossypium genotypes, Gossypium barbadense cv Phytogen 800 and G. hirsutum cv Deltapine 90. When grown in parallel, the two types of fiber developed similarly except for prolonged fiber elongation in the G. barbadense cultivar. The data were collected from isolated fibers between 10 to 28 days post anthesis (DPA) representing: primary wall synthesis to support elongation; transitional cell wall remodeling; and secondary wall cellulose synthesis, which was accompanied by continuing elongation only in G. barbadense fiber. RESULTS: Of 206 identified fiber metabolites, 205 were held in common between the two genotypes. Approximately 38,000 transcripts were expressed in the fiber of each genotype, and these were mapped to the reference set and interpreted by homology to known genes. The developmental changes in the transcriptomes and the metabolomes were compared within and across genotypes with several novel implications. Transitional cell wall remodeling is a distinct stable developmental stage lasting at least four days (18 to 21 DPA). Expression of selected cell wall related transcripts was similar between genotypes, but cellulose synthase gene expression patterns were more complex than expected. Lignification was transcriptionally repressed in both genotypes. Oxidative stress was lower in the fiber of G. barbadense cv Phytogen 800 as compared to G. hirsutum cv Deltapine 90. Correspondingly, the G. barbadense cultivar had enhanced capacity for management of reactive oxygen species during its prolonged elongation period, as indicated by a 138-fold increase in ascorbate concentration at 28 DPA. CONCLUSIONS: The parallel data on deep-sequencing transcriptomics and non-targeted metabolomics for two genotypes of single-celled cotton fiber showed that a discrete developmental stage of transitional cell wall remodeling occurs before secondary wall cellulose synthesis begins. The data showed how lignification can be transcriptionally repressed during secondary cell wall synthesis, and they implicated enhanced capacity to manage reactive oxygen species through the ascorbate-glutathione cycle as a positive contributor to fiber length.

Molecular characterization of the Gossypium Diversity Reference Set of the US National Cotton Germplasm Collection.

KEY MESSAGE: A core marker set containing markers developed to be informative within a single commercial cotton species can elucidate diversity structure within a multi-species subset of the Gossypium germplasm collection. An understanding of the genetic diversity of cotton (Gossypium spp.) as represented in the US National Cotton Germplasm Collection is essential to develop strategies for collecting, conserving, and utilizing these germplasm resources. The US collection is one of the largest world collections and includes not only accessions with improved yield and fiber quality within cultivated species, but also accessions possessing sources of abiotic and biotic stress resistance often found in wild species. We evaluated the genetic diversity of a subset of 272 diploid and 1,984 tetraploid accessions in the collection (designated the Gossypium Diversity Reference Set) using a core set of 105 microsatellite markers. Utility of the core set of markers in differentiating intra-genome variation was much greater in commercial tetraploid genomes (99.7 % polymorphic bands) than in wild diploid genomes (72.7 % polymorphic bands), and may have been influenced by pre-selection of markers for effectiveness in the commercial species. Principal coordinate analyses revealed that the marker set differentiated interspecific variation among tetraploid species, but was only capable of partially differentiating among species and genomes of the wild diploids. Putative species-specific marker bands in G. hirsutum (73) and G. barbadense (81) were identified that could be used for qualitative identification of misclassifications, redundancies, and introgression within commercial tetraploid species. The results of this broad-scale molecular characterization are essential to the management and conservation of the collection and provide insight and guidance in the use of the collection by the cotton research community in their cotton improvement efforts.

Development of a Large Set of Microsatellite Markers in Zapote Mamey (Pouteria sapota (Jacq.) H.E. Moore & Stearn) and Their Potential Use in the Study of the Species.

Pouteria sapota is known for its edible fruits that contain unique carotenoids, as well as for its fungitoxic, anti-inflammatory and anti-oxidant activity. However, its genetics is mostly unknown, including aspects about its genetic diversity and domestication process. We did high-throughput sequencing of microsatellite-enriched libraries of P. sapota, generated 5223 contig DNA sequences, 1.8 Mbp, developed 368 microsatellites markers and tested them on 29 individuals from 10 populations (seven wild, three cultivated) from Mexico, its putative domestication center. Gene ontology BLAST analysis of the DNA sequences containing microsatellites showed potential association to physiological functions. Genetic diversity was slightly higher in cultivated than in the wild gene pool (HE = 0.41 and HE = 0.35, respectively), although modified Garza-Williamson Index and Bottleneck software showed evidence for a reduction in genetic diversity for the cultivated one. Neighbor Joining, 3D Principal Coordinates Analysis and assignment tests grouped most individuals according to their geographic origin but no clear separation was observed between wild or cultivated gene pools due to, perhaps, the existence of several admixed populations. The developed microsatellites have a great potential in genetic population and domestication studies of P. sapota but additional sampling will be necessary to better understand how the domestication process has impacted the genetic diversity of this fruit crop.

Candidate gene association mapping of Sclerotinia stalk rot resistance in sunflower (Helianthus annuus L.) uncovers the importance of COI1 homologs.

KEY MESSAGE: Functional markers for Sclerotinia basal stalk rot resistance in sunflower were obtained using gene-level information from the model species Arabidopsis thaliana. Sclerotinia stalk rot, caused by Sclerotinia sclerotiorum, is one of the most destructive diseases of sunflower (Helianthus annuus L.) worldwide. Markers for genes controlling resistance to S. sclerotiorum will enable efficient marker-assisted selection (MAS). We sequenced eight candidate genes homologous to Arabidopsis thaliana defense genes known to be associated with Sclerotinia disease resistance in a sunflower association mapping population evaluated for Sclerotinia stalk rot resistance. The total candidate gene sequence regions covered a concatenated length of 3,791 bp per individual. A total of 187 polymorphic sites were detected for all candidate gene sequences, 149 of which were single nucleotide polymorphisms (SNPs) and 38 were insertions/deletions. Eight SNPs in the coding regions led to changes in amino acid codons. Linkage disequilibrium decay throughout the candidate gene regions declined on average to an r (2) = 0.2 for genetic intervals of 120 bp, but extended up to 350 bp with r (2) = 0.1. A general linear model with modification to account for population structure was found the best fitting model for this population and was used for association mapping. Both HaCOI1-1 and HaCOI1-2 were found to be strongly associated with Sclerotinia stalk rot resistance and explained 7.4 % of phenotypic variation in this population. These SNP markers associated with Sclerotinia stalk rot resistance can potentially be applied to the selection of favorable genotypes, which will significantly improve the efficiency of MAS during the development of stalk rot resistant cultivars.

Camelina seed transcriptome: a tool for meal and oil improvement and translational research.

Camelina (Camelina sativa), a Brassicaceae oilseed, has received recent interest as a biofuel crop and production platform for industrial oils. Limiting wider production of camelina for these uses is the need to improve the quality and content of the seed protein-rich meal and oil, which is enriched in oxidatively unstable polyunsaturated fatty acids that are deleterious for biodiesel. To identify candidate genes for meal and oil quality improvement, a transcriptome reference was built from 2047 Sanger ESTs and more than 2 million 454-derived sequence reads, representing genes expressed in developing camelina seeds. The transcriptome of approximately 60K transcripts from 22 597 putative genes includes camelina homologues of nearly all known seed-expressed genes, suggesting a high level of completeness and usefulness of the reference. These sequences included candidates for 12S (cruciferins) and 2S (napins) seed storage proteins (SSPs) and nearly all known lipid genes, which have been compiled into an accessible database. To demonstrate the utility of the transcriptome for seed quality modification, seed-specific RNAi lines deficient in napins were generated by targeting 2S SSP genes, and high oleic acid oil lines were obtained by targeting FATTY ACID DESATURASE 2 (FAD2) and FATTY ACID ELONGASE 1 (FAE1). The high sequence identity between Arabidopsis thaliana and camelina genes was also exploited to engineer high oleic lines by RNAi with Arabidopsis FAD2 and FAE1 sequences. It is expected that these transcriptomic data will be useful for breeding and engineering of additional camelina seed traits and for translating findings from the model Arabidopsis to an oilseed crop.

Development and characterization of simple sequence repeats for Bipolaris sorokiniana and cross transferability to related species.

Simple sequence repeats (SSR) markers were developed from a small insert genomic library for Bipolaris sorokiniana, a mitosporic fungal pathogen that causes spot blotch and root rot in switchgrass. About 59% of sequenced clones (n = 384) harbored SSR motifs. After eliminating redundant sequences, 196 SSR loci were identified, of which 84.7% were dinucleotide repeats and 9.7% and 5.6% were tri- and tetra-nucleotide repeats, respectively. Primer pairs were designed for 105 loci and 85 successfully amplified loci. Sixteen polymorphic loci were characterized with 15 B. sorokiniana isolates obtained from infected switchgrass plant materials collected from five states in USA. These loci successfully cross-amplified isolates from at least one related species, including Bipolaris oryzae, Bipolaris spicifera and Bipolaris victoriae, that causes leaf spot on switchgrass. Haploid gene diversity per locus across all isolates studied varied 0.633-0.861. Principal component analysis of SSR data clustered isolates according to their respective species. These SSR markers will be a valuable tool for genetic variability and population studies of B. sorokiniana and related species that are pathogenic on switchgrass and other host plants. In addition, these markers are potential diagnostic tools for species in the genus Bipolaris.

Chromosome-scale genome assembly of the pink bollworm, Pectinophora gossypiella, a global pest of cotton.

The pink bollworm, Pectinophora gossypiella (Saunders) (Lepidoptera: Gelechiidae), is a major global pest of cotton. Current management practices include chemical insecticides, cultural strategies, sterile insect releases, and transgenic cotton producing crystalline (Cry) protein toxins of the bacterium Bacillus thuringiensis (Bt). These strategies have contributed to the eradication of P. gossypiella from the cotton-growing areas of the United States and northern Mexico. However, this pest has evolved resistance to Bt cotton in Asia, where it remains a critical pest, and the benefits of using transgenic Bt crops have been lost. A complete annotated reference genome is needed to improve global Bt resistance management of the pink bollworm. We generated the first chromosome-level genome assembly for pink bollworm from a Bt-susceptible laboratory strain (APHIS-S) using PacBio continuous long reads for contig generation, Illumina Hi-C for scaffolding, and Illumina whole-genome re-sequencing for error correction. The pseudo-haploid assembly consists of 29 autosomes and the Z sex chromosome. The assembly exceeds the minimum Earth BioGenome Project quality standards, has a low error rate, is highly contiguous at both the contig and scaffold levels (L/N50 of 18/8.26 MB and 14/16.44 MB, respectively), and is complete, with 98.6% of lepidopteran single-copy orthologs represented without duplication. The genome was annotated with 50% repeat content and 14,107 protein-coding genes, further assigned to 41,666 functional annotations. This assembly represents the first publicly available complete annotated genome of pink bollworm and will serve as the foundation for advancing molecular genetics of this important pest species.

Representing true plant genomes: haplotype-resolved hybrid pepper genome with trio-binning.

As sequencing costs decrease and availability of high fidelity long-read sequencing increases, generating experiment specific de novo genome assemblies becomes feasible. In many crop species, obtaining the genome of a hybrid or heterozygous individual is necessary for systems that do not tolerate inbreeding or for investigating important biological questions, such as hybrid vigor. However, most genome assembly methods that have been used in plants result in a merged single sequence representation that is not a true biologically accurate representation of either haplotype within a diploid individual. The resulting genome assembly is often fragmented and exhibits a mosaic of the two haplotypes, referred to as haplotype-switching. Important haplotype level information, such as causal mutations and structural variation is therefore lost causing difficulties in interpreting downstream analyses. To overcome this challenge, we have applied a method developed for animal genome assembly called trio-binning to an intra-specific hybrid of chili pepper (Capsicum annuum L. cv. HDA149 x Capsicum annuum L. cv. HDA330). We tested all currently available softwares for performing trio-binning, combined with multiple scaffolding technologies including Bionano to determine the optimal method of producing the best haplotype-resolved assembly. Ultimately, we produced highly contiguous biologically true haplotype-resolved genome assemblies for each parent, with scaffold N50s of 266.0 Mb and 281.3 Mb, with 99.6% and 99.8% positioned into chromosomes respectively. The assemblies captured 3.10 Gb and 3.12 Gb of the estimated 3.5 Gb chili pepper genome size. These assemblies represent the complete genome structure of the intraspecific hybrid, as well as the two parental genomes, and show measurable improvements over the currently available reference genomes. Our manuscript provides a valuable guide on how to apply trio-binning to other plant genomes.

Two haplotype-resolved genomes reveal important flower traits in bigleaf hydrangea (Hydrangea macrophylla) and insights into Asterid evolution.

The Hydrangea genus belongs to the Hydrangeaceae family, in the Cornales order of flowering plants, which early diverged among the Asterids, and includes several species that are commonly used ornamental plants. Of them, Hydrangea macrophylla is one of the most valuable species in the nursery trade, yet few genomic resources are available for this crop or closely related Asterid species. Two high-quality haplotype-resolved reference genomes of hydrangea cultivars 'Veitchii' and 'Endless Summer' [highest quality at 2.22 gigabase pairs (Gb), 396 contigs, N50 22.8 megabase pairs (Mb)] were assembled and scaffolded into the expected 18 pseudochromosomes. Utilizing the newly developed high-quality reference genomes along with high-quality genomes of other related flowering plants, nuclear data were found to support a single divergence point in the Asterids clade where both the Cornales and Ericales diverged from the euasterids. Genetic mapping with an F1 hybrid population demonstrated the power of linkage mapping combined with the new genomic resources to identify the gene for inflorescence shape, CYP78A5 located on chromosome 4, and a novel gene, BAM3 located on chromosome 17, for causing double flower. Resources developed in this study will not only help to accelerate hydrangea genetic improvement but also contribute to understanding the largest group of flowering plants, the Asterids.