Seminal plasma uterine priming alters uterine transcriptomics and negatively impacts embryo growth and uterine artery resistance but not offspring liver transcriptomics in beef cattle.

Seminal plasma uterine priming is important for pregnancy and offspring phenotype in mice and swine; however, impacts on the uterus of the dam and her offspring in cattle are unknown. We sought to determine the effects of seminal plasma uterine priming at estrus on uterine transcriptomics, early gestation (d 35, 40, and 45) embryo morphometrics, mid- to late-gestation (d 140 to 220) uterine artery hemodynamics, birth morphometrics, and liver transcriptomics in offspring at 30 d of age. Multiparous Angus-based commercial beef cows were randomly assigned to receive treatment at estrus: 0.5 mL pooled seminal plasma in the uterine body (n = 31, seminal plasma primed) or no treatment (n = 31, control). Seven d later a subset of cows (n = 4/treatment) underwent uterine biopsies, and the remaining cows underwent embryo transfer. Embryo crown-rump length and uterine artery hemodynamics were measured during gestation using ultrasonography. Morphometrics of the calf were collected within 24 h of parturition. Liver biopsies were collected at 30 d of age. Data were analyzed by ANOVA in a completely randomized design for the effect of treatment. Myosin heavy chain I (JSP.1) was downregulated [Benjamin-Hochberg adj P (BH) <= 0.05] and ABO alpha 1-3-N-acetylgalactosaminyltransferase and alpha 1-3-galactosyltransferase (ABO) was upregulated (BH adj P <= 0.05) in the uterus of seminal plasma primed cows 7 d after treatment. Embryo crown-rump length was less (P < 0.05) in seminal plasma primed cows. Mid- to late-gestation (d 140 to 220) uterine artery resistance was increased (P < 0.05) in seminal plasma primed cows. Seminal plasma priming did not alter birth weights or curve-crown-rump length, but heart girth was increased (P < 0.05) in offspring from seminal plasma primed cows. There were no differentially expressed genes (BH adj P <= 0.05) in offspring liver at 30 d of age; however, myosin light chain, phosphorylatable, fast skeletal muscle (MYLPF) was absent in all liver samples from calves from seminal plasma primed cows. In contrast, vomeronasal 1 receptor bosTauV1R414 (BOSTAUV1R414) was present in 6 of the 7 liver samples from calves from seminal plasma primed cows. Seminal plasma uterine priming alters uterine transcriptomics, negatively impacts early gestation embryo growth and mid- to late-gestation uterine artery resistance suggesting downstream vascular anomalies. However, these in utero conditions did not impact offspring from birth to 30 d of age.

Origin and diversity of the wild cottons (Gossypium hirsutum) of Mound Key, Florida.

Elucidating genetic diversity within wild forms of modern crops is essential for understanding domestication and the possibilities of wild germplasm utilization. Gossypium hirsutum is a predominant source of natural plant fibers and the most widely cultivated cotton species. Wild forms of G. hirsutum are challenging to distinguish from feral derivatives, and truly wild populations are uncommon. Here we characterize a population from Mound Key Archaeological State Park, Florida using genome-wide SNPs extracted from 25 individuals over three sites. Our results reveal that this population is genetically dissimilar from other known wild, landrace, and domesticated cottons, and likely represents a pocket of previously unrecognized wild genetic diversity. The unexpected level of divergence between the Mound Key population and other wild cotton populations suggests that the species may harbor other remnant and genetically distinct populations that are geographically scattered in suitable habitats throughout the Caribbean. Our work thus has broader conservation genetic implications and suggests that further exploration of natural diversity in this species is warranted.

The Gossypium herbaceum L. Wagad genome as a resource for understanding cotton domestication.

Gossypium herbaceum is a species of cotton native to Africa and Asia that is one of the 2 domesticated diploids. Together with its sister-species G. arboreum, these A-genome taxa represent models of the extinct A-genome donor of modern polyploid cotton, which provide about 95% of cotton grown worldwide. As part of a larger effort to characterize variation and improve resources among diverse diploid and polyploid cotton genomes, we sequenced and assembled the genome of G. herbaceum cultivar (cv.) Wagad, representing the first domesticated accession for this species. This chromosome-level genome was generated using a combination of PacBio long-read technology, HiC, and Bionano optical mapping and compared to existing genome sequences in cotton. We compare the genome of this cultivar to the existing genome of wild G. herbaceum subspecies africanum to elucidate changes in the G. herbaceum genome concomitant with domestication and extend these analyses to gene expression using available RNA-seq. Our results demonstrate the utility of the G. herbaceum cv. Wagad genome in understanding domestication in the diploid species, which could inform modern breeding programs.

Variation in cytonuclear expression accommodation among allopolyploid plants.

Cytonuclear coevolution is a common feature among plants, which coordinates gene expression and protein products between the nucleus and organelles. Consequently, lineage-specific differences may result in incompatibilities between the nucleus and cytoplasm in hybrid taxa. Allopolyploidy is also a common phenomenon in plant evolution. The hybrid nature of allopolyploids may result in cytonuclear incompatibilities, but the massive nuclear redundancy created during polyploidy affords additional avenues for resolving cytonuclear conflict (i.e. cytonuclear accommodation). Here we evaluate expression changes in organelle-targeted nuclear genes for 6 allopolyploid lineages that represent 4 genera (i.e. Arabidopsis, Arachis, Chenopodium, and Gossypium) and encompass a range in polyploid ages. Because incompatibilities between the nucleus and cytoplasm could potentially result in biases toward the maternal homoeolog and/or maternal expression level, we evaluate patterns of homoeolog usage, expression bias, and expression-level dominance in cytonuclear genes relative to the background of noncytonuclear expression changes and to the diploid parents. Although we find subsets of cytonuclear genes in most lineages that match our expectations of maternal preference, these observations are not consistent among either allopolyploids or categories of organelle-targeted genes. Our results indicate that cytonuclear expression evolution may be subtle and variable among genera and genes, likely reflecting a diversity of mechanisms to resolve nuclear-cytoplasmic incompatibilities in allopolyploid species.

Organellar transcripts dominate the cellular mRNA pool across plants of varying ploidy levels.

Mitochondrial and plastid functions depend on coordinated expression of proteins encoded by genomic compartments that have radical differences in copy number of organellar and nuclear genomes. In polyploids, doubling of the nuclear genome may add challenges to maintaining balanced expression of proteins involved in cytonuclear interactions. Here, we use ribo-depleted RNA sequencing (RNA-seq) to analyze transcript abundance for nuclear and organellar genomes in leaf tissue from four different polyploid angiosperms and their close diploid relatives. We find that even though plastid genomes contain <1% of the number of genes in the nuclear genome, they generate the majority (69.9 to 82.3%) of messenger RNA (mRNA) transcripts in the cell. Mitochondrial genes are responsible for a much smaller percentage (1.3 to 3.7%) of the leaf mRNA pool but still produce much higher transcript abundances per gene compared to nuclear genome. Nuclear genes encoding proteins that functionally interact with mitochondrial or plastid gene products exhibit mRNA expression levels that are consistently more than 10-fold lower than their organellar counterparts, indicating an extreme cytonuclear imbalance at the RNA level despite the predominance of equimolar interactions at the protein level. Nevertheless, interacting nuclear and organellar genes show strongly correlated transcript abundances across functional categories, suggesting that the observed mRNA stoichiometric imbalance does not preclude coordination of cytonuclear expression. Finally, we show that nuclear genome doubling does not alter the cytonuclear expression ratios observed in diploid relatives in consistent or systematic ways, indicating that successful polyploid plants are able to compensate for cytonuclear perturbations associated with nuclear genome doubling.

Homoeologous gene expression and co-expression network analyses and evolutionary inference in allopolyploids.

Polyploidy is a widespread phenomenon throughout eukaryotes. Due to the coexistence of duplicated genomes, polyploids offer unique challenges for estimating gene expression levels, which is essential for understanding the massive and various forms of transcriptomic responses accompanying polyploidy. Although previous studies have explored the bioinformatics of polyploid transcriptomic profiling, the causes and consequences of inaccurate quantification of transcripts from duplicated gene copies have not been addressed. Using transcriptomic data from the cotton genus (Gossypium) as an example, we present an analytical workflow to evaluate a variety of bioinformatic method choices at different stages of RNA-seq analysis, from homoeolog expression quantification to downstream analysis used to infer key phenomena of polyploid expression evolution. In general, EAGLE-RC and GSNAP-PolyCat outperform other quantification pipelines tested, and their derived expression dataset best represents the expected homoeolog expression and co-expression divergence. The performance of co-expression network analysis was less affected by homoeolog quantification than by network construction methods, where weighted networks outperformed binary networks. By examining the extent and consequences of homoeolog read ambiguity, we illuminate the potential artifacts that may affect our understanding of duplicate gene expression, including an overestimation of homoeolog co-regulation and the incorrect inference of subgenome asymmetry in network topology. Taken together, our work points to a set of reasonable practices that we hope are broadly applicable to the evolutionary exploration of polyploids.

The Gossypium longicalyx Genome as a Resource for Cotton Breeding and Evolution.

Cotton is an important crop that has made significant gains in production over the last century. Emerging pests such as the reniform nematode have threatened cotton production. The rare African diploid species Gossypium longicalyx is a wild species that has been used as an important source of reniform nematode immunity. While mapping and breeding efforts have made some strides in transferring this immunity to the cultivated polyploid species, the complexities of interploidal transfer combined with substantial linkage drag have inhibited progress in this area. Moreover, this species shares its most recent common ancestor with the cultivated A-genome diploid cottons, thereby providing insight into the evolution of long, spinnable fiber. Here we report a newly generated de novo genome assembly of G. longicalyx This high-quality genome leveraged a combination of PacBio long-read technology, Hi-C chromatin conformation capture, and BioNano optical mapping to achieve a chromosome level assembly. The utility of the G. longicalyx genome for understanding reniform immunity and fiber evolution is discussed.

Therapeutic potential of N-acetylcysteine in acrylamide acute neurotoxicity in adult zebrafish.

Two essential key events in acrylamide (ACR) acute neurotoxicity are the formation of adducts with nucleophilic sulfhydryl groups on cysteine residues of selected proteins in the synaptic terminals and the depletion of the glutathione (GSx) stores in neural tissue. The use of N-acetylcysteine (NAC) has been recently proposed as a potential antidote against ACR neurotoxicity, as this chemical is not only a well-known precursor of the reduced form of glutathione (GSH), but also is an scavenger of soft electrophiles such as ACR. In this study, the suitability of 0.3 and 0.75 mM NAC to protect against the neurotoxic effect of 0.75 mM ACR has been tested in vivo in adult zebrafish. NAC provided only a mild to negligible protection against the changes induced by ACR in the motor function, behavior, transcriptome and proteome. The permeability of NAC to cross blood-brain barrier (BBB) was assessed, as well as the ACR-scavenging activity and the gamma-glutamyl-cysteine ligase (γ-GCL) and acylase I activities. The results show that ACR not only depletes GSx levels but also inhibits it synthesis from NAC/cysteine, having a dramatic effect over the glutathione system. Moreover, results indicate a very low NAC uptake to the brain, probably by a combination of low BBB permeability and high deacylation of NAC during the intestinal absorption. These results strongly suggest that the use of NAC is not indicated in ACR acute neurotoxicity treatment.

Independent Domestication of Two Old World Cotton Species.

Domesticated cotton species provide raw material for the majority of the world's textile industry. Two independent domestication events have been identified in allopolyploid cotton, one in Upland cotton (Gossypium hirsutum L.) and the other to Egyptian cotton (Gossypium barbadense L.). However, two diploid cotton species, Gossypium arboreum L. and Gossypium herbaceum L., have been cultivated for several millennia, but their status as independent domesticates has long been in question. Using genome resequencing data, we estimated the global abundance of various repetitive DNAs. We demonstrate that, despite negligible divergence in genome size, the two domesticated diploid cotton species contain different, but compensatory, repeat content and have thus experienced cryptic alterations in repeat abundance despite equivalence in genome size. Evidence of independent origin is bolstered by estimates of divergence times based on molecular evolutionary analysis of f7,000 orthologous genes, for which synonymous substitution rates suggest that G. arboreum and G. herbaceum last shared a common ancestor approximately 0.4-2.5 Ma. These data are incompatible with a shared domestication history during the emergence of agriculture and lead to the conclusion that G. arboreum and G. herbaceum were each domesticated independently.