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.

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.

Two pathogen loci determine Blumeria graminis f. sp. tritici virulence to wheat resistance gene Pm1a.

Blumeria graminis f. sp. tritici (Bgt) is a globally important fungal pathogen of wheat that can rapidly evolve to defeat wheat powdery mildew (Pm) resistance genes. Despite periodic regional deployment of the Pm1a resistance gene in US wheat production, Bgt strains that overcome Pm1a have been notably nonpersistent in the United States, while on other continents, they are more widely established. A genome-wide association study (GWAS) was conducted to map sequence variants associated with Pm1a virulence in 216 Bgt isolates from six countries, including the United States. A virulence variant apparently unique to Bgt isolates from the United States was detected in the previously mapped gene AvrPm1a (BgtE-5612) on Bgt chromosome 6; an in vitro growth assay suggested no fitness reduction associated with this variant. A gene on Bgt chromosome 8, Bgt-51526, was shown to function as a second determinant of Pm1a virulence, and despite < 30% amino acid identity, BGT-51526 and BGTE-5612 were predicted to share > 85% of their secondary structure. A co-expression study in Nicotiana benthamiana showed that BGTE-5612 and BGT-51526 each produce a PM1A-dependent hypersensitive response. More than one member of a B. graminis effector family can be recognized by a single wheat immune receptor, and a two-gene model is necessary to explain virulence to Pm1a.

Enhancing Upland cotton for drought resilience, productivity, and fiber quality: comparative evaluation and genetic dissection.

To provision the world sustainably, modern society must increase overall crop production, while conserving and preserving natural resources. Producing more with diminishing water resources is an especially daunting endeavor. Toward the goal of genetically improving drought resilience of cultivated Upland cotton (Gossypium hirsutum L.), this study addresses the genetics of differential yield components referred to as productivity and fiber quality traits under regular-water versus low-water (LW) field conditions. We used ten traits to assess water stress deficit, which included six productivity and four fiber quality traits on two recombinant inbred line (RIL) populations from reciprocally crossed cultivars, Phytogen 72 and Stoneville 474. To facilitate genetic inferences, we genotyped RILs with the CottonSNP63K array, assembled high-density linkage maps of over 7000 SNPs and then analyzed quantitative trait variations. Analysis of variance revealed significant differences for all traits (p < 0.05) in these RIL populations. Although the LW irrigation regime significantly reduced all traits, except lint percent, the RILs exhibited a broad phenotypic spectrum of heritable differences across the water regimes. Transgressive segregation occurred among the RILs, suggesting the possibility of genetic gain through phenotypic selection for drought resilience and perhaps through marker-based selection. Analyses revealed more than 150 quantitative trait loci (QTLs) associated with productivity and fiber quality traits (p < 0.005) on different genomic regions of the cotton genome. The multiple-QTL models analysis with LOD > 3.0 detected 21 QTLs associated with productivity and 22 QTLs associated with fiber quality. For fiber traits, strong clustering and QTL associations occurred in c08 and its homolog c24 as well as c10, c14, and c21. Using contemporary genome sequence assemblies and bioinformatically related information, the identification of genomic regions associated with responses to plant stress/drought elevates the possibility of using marker-assisted and omics-based selection to enhance breeding for drought resilient cultivars and identifying candidate genes and networks. RILs with different responses to drought indicated that it is possible to maintain high fiber quality under LW conditions or reduce the of LW impact on quality. The heritable variation among elite bi-parental RILs for productivity and quality under field drought conditions, and their association of QTLs, and thus specific genomic regions, indicate opportunities for breeding-based gains in water resource conservation, i.e., enhancing cotton's agricultural sustainability.

Correction: DNA Sequence Evolution and Rare Homoeologous Conversion in Tetraploid Cotton.

[This corrects the article DOI: 10.1371/journal.pgen.1006012.].

DNA Sequence Evolution and Rare Homoeologous Conversion in Tetraploid Cotton.

Allotetraploid cotton species are a vital source of spinnable fiber for textiles. The polyploid nature of the cotton genome raises many evolutionary questions as to the relationships between duplicated genomes. We describe the evolution of the cotton genome (SNPs and structural variants) with the greatly improved resolution of 34 deeply re-sequenced genomes. We also explore the evolution of homoeologous regions in the AT- and DT-genomes and especially the phenomenon of conversion between genomes. We did not find any compelling evidence for homoeologous conversion between genomes. These findings are very different from other recent reports of frequent conversion events between genomes. We also identified several distinct regions of the genome that have been introgressed between G. hirsutum and G. barbadense, which presumably resulted from breeding efforts targeting associated beneficial alleles. Finally, the genotypic data resulting from this study provides access to a wealth of diversity sorely needed in the narrow germplasm of cotton cultivars.

Comparative transcriptomics and genomic patterns of discordance in Capsiceae (Solanaceae).

The integration of genomics and phylogenetics allows new insight into the structure of gene tree discordance, the relationships among gene position, gene history, and rate of evolution, as well as the correspondence of gene function, positive selection, and gene ontology enrichment across lineages. We explore these issues using the tribe Capsiceae (Solanaceae), which is comprised of the genera Lycianthes and Capsicum (peppers). In combining the annotated genomes of Capsicum with newly sequenced transcriptomes of four species of Lycianthes and Capsicum, we develop phylogenies for 6747 genes, and construct a backbone species tree using both concordance and explicit phylogenetic network approaches. We quantify phylogenetic discordance among individual gene trees, measure their rates of synonymous and nonsynonymous substitution, and test whether they were positively selected along any branch of the phylogeny. We then map these genes onto the annotated Capsicum genome and test whether rates of evolution, gene history, and gene ontology vary significantly with gene position. We observed substantial discordance among gene trees. A bifurcating species tree placing Capsicum within a paraphyletic Lycianthes was supported over all phylogenetic networks. Rates of synonymous and nonsynonymous substitution varied 41-fold and 130-fold among genes, respectively, and were significantly lower in pericentromeric regions. We found that results of concordance tree analyses vary depending on the subset of genes used, and that genes within the pericentromeric regions only capture a portion of the observed discordance. We identified 787 genes that have been positively selected throughout the diversification history of Capsiceae, and discuss the importance of these genes as targets for investigation of economically important traits in the domesticated peppers.

Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array.

BACKGROUND: Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits. RESULTS: The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content. CONCLUSIONS: Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton.

Development and bin mapping of gene-associated interspecific SNPs for cotton (Gossypium hirsutum L.) introgression breeding efforts.

BACKGROUND: Cotton (Gossypium spp.) is the largest producer of natural fibers for textile and is an important crop worldwide. Crop production is comprised primarily of G. hirsutum L., an allotetraploid. However, elite cultivars express very small amounts of variation due to the species monophyletic origin, domestication and further bottlenecks due to selection. Conversely, wild cotton species harbor extensive genetic diversity of prospective utility to improve many beneficial agronomic traits, fiber characteristics, and resistance to disease and drought. Introgression of traits from wild species can provide a natural way to incorporate advantageous traits through breeding to generate higher-producing cotton cultivars and more sustainable production systems. Interspecific introgression efforts by conventional methods are very time-consuming and costly, but can be expedited using marker-assisted selection. RESULTS: Using transcriptome sequencing we have developed the first gene-associated single nucleotide polymorphism (SNP) markers for wild cotton species G. tomentosum, G. mustelinum, G. armourianum and G. longicalyx. Markers were also developed for a secondary cultivated species G. barbadense cv. 3-79. A total of 62,832 non-redundant SNP markers were developed from the five wild species which can be utilized for interspecific germplasm introgression into cultivated G. hirsutum and are directly associated with genes. Over 500 of the G. barbadense markers have been validated by whole-genome radiation hybrid mapping. Overall 1,060 SNPs from the five different species have been screened and shown to produce acceptable genotyping assays. CONCLUSIONS: This large set of 62,832 SNPs relative to cultivated G. hirsutum will allow for the first high-density mapping of genes from five wild species that affect traits of interest, including beneficial agronomic and fiber characteristics. Upon mapping, the markers can be utilized for marker-assisted introgression of new germplasm into cultivated cotton and in subsequent breeding of agronomically adapted types, including cultivar development.

Gossypium hirsutum gene of unknown function Gohir.A02G131900 encodes a potential plant-specific, dual-domain exo-1,3-ß-glucosidase.

A gene of unknown function, Gohir.A02G131900.1, identified in Gossypium hirsutum was studied using computational sequence and structure bioinformatic tools. The encoded protein GhGH5BG-A0A1U8NW40 (UniProt A0A1U8NW40) is predicted to be secreted and localized to the cell wall. Homology and conserved residues indicate it belongs to a plant-specific subgroup of the glycoside hydrolase family 5 and likely has exo-1,3-ß-glucosidase activity. This subgroup is unique in containing a fascin-like domain which may have evolved a unique glucan binding site of interest for further research.

Gossypium hirsutum gene of unknown function Gohir.A02G161000 encodes a potential transmembrane Root UVB Sensitive 4 Protein with a putative protein-protein interaction interface.

A gene of unknown function, Gohir.A02G161000.1, identified in Gossypium hirsutum was studied using computational sequence and structure bioinformatics tools. The associated protein GhRUS4-A0A1U8JPV7 (UniProt A0A1U8JPV7) is predicted to be a plastid-localized, transmembrane root UVB-sensitive 4 (RUS4) protein with a newly identified potential dimerization surface. Evidence from homology and sequence conservation suggest involvement in auxin transport and pollen maturation.

Gossypium hirsutum gene of unknown function, Gohir.A02G039501.1, encodes a potential DNA-binding ALOG protein involved in gene regulation.

A protein of unknown function encoded by gene Gohir.A02G039501.1 in Gossypium hirsutum , was studied using sequence and structure bioinformatic tools leading to its proposed function as a nuclear, DNA-binding ALOG protein involved in gene regulation during organ boundary specification and maintenance. The encoded protein contains a predicted nuclear localization sequence, an ALOG domain with conserved residues in the modeled DNA-binding regions and nearly identical sequence identity to Arabidopsis homologs involved in development of organ boundaries at the shoot apical meristem. The protein was modeled by AlphaFold2 to have a four-helix bundle that is structurally analogous to DNA-binding domains of XerC/D-like recombinases.

Gossypium hirsutum gene of unknown function Gohir.A03G007700.1 encodes a potential VAN3-binding protein with a phosphoinositide-binding site.

A gene of unknown function, Gohir.A03G007700.1 (gene ID: Gohir.A03G007700_UTX-TM1_v2.1; transcript ID: Gohir.A03G007700.1_UTX-TM1_v2.1), identified in Gossypium hirsutum was studied using bioinformatic analyses of the sequence and structure of its encoded protein. Results from domain prediction, conserved residues and structure comparison indicate the encoded plant-specific protein (UniProt A0A1U8N485) is part of the VAN3-binding protein family with a conserved phosphoinositide-binding site. Homology comparison suggests functional similarity with Arabidopsis FORKED-like FL5 and 6, which localize to the Golgi apparatus and are linked to vein development and leaf size phenotypes.

A first look at the ability to use genomic prediction for improving the ratooning ability of sugarcane.

The sugarcane ratooning ability (RA) is the most important target trait for breeders seeking to enhance the profitability of sugarcane production by reducing the planting cost. Understanding the genetics governing the RA could help breeders by identifying molecular markers that could be used for genomics-assisted breeding (GAB). A replicated field trial was conducted for three crop cycles (plant cane, first ratoon, and second ratoon) using 432 sugarcane clones and used for conducting genome-wide association and genomic prediction of five sugar and yield component traits of the RA. The RA traits for economic index (EI), stalk population (SP), stalk weight (SW), tonns of cane per hectare (TCH), and tonns of sucrose per hectare (TSH) were estimated from the yield and sugar data. A total of six putative quantitative trait loci and eight nonredundant single-nucleotide polymorphism (SNP) markers were associated with all five tested RA traits and appear to be unique. Seven putative candidate genes were colocated with significant SNPs associated with the five RA traits. The genomic prediction accuracies for those tested traits were moderate and ranged from 0.21 to 0.36. However, the models fitting fixed effects for the most significant associated markers for each respective trait did not give any advantages over the standard models without fixed effects. As a result of this study, more robust markers could be used in the future for clone selection in sugarcane, potentially helping resolve the genetic control of the RA in sugarcane.

Gossypium hirsutum gene of unknown function Gohir.A03G0737001 encodes a potential Chaperone-like Protein of protochlorophyllide oxidoreductase (CPP1).

A gene of unknown function identified in Gossypium hirsutum , Gohir.A03G0737001.1, was studied using sequence and bioinformatic tools. The encoded protein (referred to here as GhCPP1-A0A1U8HKT6) was predicted to function as a Chaperone-like protein of protochlorophyllide oxidoreductase (CPP1), which is involved with initiation of photochemical reactions of chlorophyll biosynthesis. Sequence analysis indicates it is embedded in the chloroplast envelope membrane through four transmembrane regions and contains a J-like domain that is structurally similar to the J domain of DnaJ/Hsp40 "holdase" chaperone proteins.

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.

Detecting Cotton Leaf Curl Virus Resistance Quantitative Trait Loci in Gossypium hirsutum and iCottonQTL a New R/Shiny App to Streamline Genetic Mapping.

Cotton leaf curl virus (CLCuV) causes devastating losses to fiber production in Central Asia. Viral spread across Asia in the last decade is causing concern that the virus will spread further before resistant varieties can be bred. Current development depends on screening each generation under disease pressure in a country where the disease is endemic. We utilized quantitative trait loci (QTL) mapping in four crosses with different sources of resistance to identify single nucleotide polymorphism (SNP) markers associated with the resistance trait to allow development of varieties without the need for field screening every generation. To assist in the analysis of multiple populations, a new publicly available R/Shiny App was developed to streamline genetic mapping using SNP arrays and to also provide an easy method to convert and deposit genetic data into the CottonGen database. Results identified several QTL from each cross, indicating possible multiple modes of resistance. Multiple sources of resistance would provide several genetic routes to combat the virus as it evolves over time. Kompetitive allele specific PCR (KASP) markers were developed and validated for a subset of QTL, which can be used in further development of CLCuV-resistant cotton lines.

Gossypium hirsutum gene of unknown function, Gohir.A02G044702.1, encodes a potential B3 Transcription Factor of the REM subfamily.

A gene of unknown function, Gohir.A02G044702.1, identified in Gossypium hirsutum was studied using sequence and structure bioinformatic tools. The encoded protein (UniProt A0A1U8MGX4) was predicted to localize to the nucleus, was found to retain the B3 transcription factor domain with conserved DNA-binding residues and to most closely cluster with REM subfamily members of B3-domain containing proteins.

Outlook for Implementation of Genomics-Based Selection in Public Cotton Breeding Programs.

Researchers have used quantitative genetics to map cotton fiber quality and agronomic performance loci, but many alleles may be population or environment-specific, limiting their usefulness in a pedigree selection, inbreeding-based system. Here, we utilized genotypic and phenotypic data on a panel of 80 important historical Upland cotton (Gossypium hirsutum L.) lines to investigate the potential for genomics-based selection within a cotton breeding program's relatively closed gene pool. We performed a genome-wide association study (GWAS) to identify alleles correlated to 20 fiber quality, seed composition, and yield traits and looked for a consistent detection of GWAS hits across 14 individual field trials. We also explored the potential for genomic prediction to capture genotypic variation for these quantitative traits and tested the incorporation of GWAS hits into the prediction model. Overall, we found that genomic selection programs for fiber quality can begin immediately, and the prediction ability for most other traits is lower but commensurate with heritability. Stably detected GWAS hits can improve prediction accuracy, although a significance threshold must be carefully chosen to include a marker as a fixed effect. We place these results in the context of modern public cotton line-breeding and highlight the need for a community-based approach to amass the data and expertise necessary to launch US public-sector cotton breeders into the genomics-based selection era.