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.

Identification of Protein Biomarkers for Differentiating Listeria monocytogenes Genetic Lineage III.

Listeria monocytogenes is the causative agent of listeriosis, a severe foodborne illness characterized by septicemia, meningitis, encephalitis, abortions, and occasional death in infants and immunocompromised individuals. L. monocytogenes is composed of four genetic lineages (I, II, III, and IV) and fourteen serotypes. The aim of the current study was to identify proteins that can serve as biomarkers for detection of genetic lineage III strains based on simple antibody-based methods. Liquid chromatography (LC) with electrospray ionization tandem mass spectrometry (ESI MS/MS) followed by bioinformatics and computational analysis were performed on three L. monocytogenes strains (NRRL B-33007, NRRL B-33014, and NRRL B-33077), which were used as reference strains for lineages I, II, and III, respectively. Results from ESI MS/MS revealed 42 unique proteins present in NRRL B-33077 and absent in NRRL B-33007 and NRRL B-33014 strains. BLAST analysis of the 42 proteins against a broader panel of >80 sequenced strains from lineages I and II revealed four proteins [TM2 domain-containing protein (NRRL B-33077_2770), DUF3916 domain-containing protein (NRRL B-33077_1897), DNA adenine methylase (NRRL B-33077_1926), and protein RhsA (NRRL B-33077_1129)] that have no homology with any sequenced strains in lineages I and II. The four genes that encode these proteins were expressed in Escherichia coli strain DE3 and purified. Polyclonal antibodies were prepared against purified recombinant proteins. ELISA using the polyclonal antibodies against 12 L. monocytogenes lineage I, II, and III isolates indicated that TM2 protein and DNA adenine methylase (Dam) detected all lineage III strains with no reaction to lineage I and II strains. In conclusion, two proteins including TM2 protein and Dam are potentially useful biomarkers for detection and differentiation of L. monocytogenes lineage III strains in clinical, environmental, and food processing facilities. Furthermore, these results validate the approach of using a combination of proteomics and bioinformatics to identify useful protein biomarkers.

Thermotolerance capabilities, blood metabolomics, and mammary gland hemodynamics and transcriptomic profiles of slick-haired Holstein cattle during mid lactation in Puerto Rico.

Holstein cattle carrying a prolactin receptor gene mutation (SLICK) exhibit short and sleek hair coats (short-haired Holstein [SLK]) enhancing thermotolerance and productivity compared with wild type-haired Holstein (WT) under tropical conditions. The objectives were to unravel the physiological and molecular mechanisms that confer an advantage to this slick genotype in Puerto Rico and determine potential correlations between metabolites and physiological variables. At 160 ± 3 DIM we compared vaginal temperatures (VT) and voluntary solar radiation exposure (VSRE) during 48 h between 9 SLK and 9 WT Holsteins, whereas a subsample of 7 SLK and 7 WT were used to assess udder skin temperature, mammary gland hemodynamics and transcriptomics, and blood plasma untargeted metabolomics at a single time point. The SLK cattle showed lower VT throughout the day and greater VSRE at 1000 h and 1100 h compared with their WT counterparts. Total mammary blood flow (MBF) was greater in SLK Holsteins compared with WT. The metabolite 9-nitrooctadecenoic acid was identified as a potential biomarker for MBF; moreover, SLK cattle had greater amounts of this metabolite in their plasma. Prostaglandin D2 synthase (PTGS) was upregulated in the slick mammary gland, while plasma prostaglandin D2 was positively correlated with milk yield and increased in SLK Holsteins compared with WT. Interestingly, the arachidonic acid metabolism pathway was enriched in the mammary gland transcriptome and perturbed in the blood metabolome in the SLK Holsteins. In conclusion, SLK Holsteins exhibited lower body temperatures, greater VSRE, enhanced blood supply to the mammary gland, and alterations in genes and metabolites involved in arachidonic acid metabolism at the mammary gland and blood plasma. The usage of the SLK Holstein cattle genetics in dairy operations could be a feasible alternative to mitigate the adverse consequences of heat stress.

In silico prediction and expression analysis of vaccine candidate genes of Campylobacter jejuni.

Campylobacter jejuni (C. jejuni) is the most common food-borne pathogen that causes human gastroenteritis in the United States. Consumption of contaminated poultry products is considered as the major source of human Campylobacter infection. An effective vaccine would be a promising alternative to antibiotic supplements to curb C. jejuni colonization in poultry gastrointestinal (GI) tract. However, the genetic diversity among the C. jejuni isolates makes vaccine production more challenging. Despite many attempts, an effective Campylobacter vaccine is not yet available. This study aimed to identify suitable candidates to develop a subunit vaccine against C. jejuni, which could reduce colonization in the GI tract of the poultry. In the current study, 4 C. jejuni strains were isolated from retail chicken meat and poultry litter samples and their genomes were sequenced utilizing next-generation sequencing technology. The genomic sequences of C. jejuni strains were screened to identify potential antigens utilizing the reverse vaccinology approach. In silico genome analysis predicted 3 conserved potential vaccine candidates (phospholipase A [PldA], TonB dependent vitamin B12 transporter [BtuB], and cytolethal distending toxin subunit B [CdtB]) suitable for the development of a vaccine. Furthermore, the expression of predicted genes during host-pathogen interaction was analyzed by an infection study using an avian macrophage-like immortalized cell line (HD11). The HD11 was infected with C. jejuni strains, and the RT-qPCR assay was performed to determine the expression of the predicted genes. The expression difference was analyzed using ΔΔCt methods. The results indicate that all 3 predicted genes, PldA, BtuB, and CdtB, were upregulated in 4 tested C. jejuni strains irrespective of their sources of isolation. In conclusion, in silico prediction and gene expression analysis during host-pathogen interactions identified 3 potential vaccine candidates for C. jejuni.

A high-quality chromosome-level genome assembly of rohu carp, Labeo rohita, and its utilization in SNP-based exploration of gene flow and sex determination.

Labeo rohita (rohu) is a carp important to aquaculture in South Asia, with a production volume close to Atlantic salmon. While genetic improvements to rohu are ongoing, the genomic methods commonly used in other aquaculture improvement programs have historically been precluded in rohu, partially due to the lack of a high-quality reference genome. Here we present a high-quality de novo genome produced using a combination of next-generation sequencing technologies, resulting in a 946 Mb genome consisting of 25 chromosomes and 2,844 unplaced scaffolds. Notably, while approximately half the size of the existing genome sequence, our genome represents 97.9% of the genome size newly estimated here using flow cytometry. Sequencing from 120 individuals was used in conjunction with this genome to predict the population structure, diversity, and divergence in three major rivers (Jamuna, Padma, and Halda), in addition to infer a likely sex determination mechism in rohu. These results demonstrate the utility of the new rohu genome in modernizing some aspects of rohu genetic improvement programs.

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.

Innovations in double digest restriction-site associated DNA sequencing (ddRAD-Seq) method for more efficient SNP identification.

We present an improved ddRAD-Seq protocol for identifying single nucleotide polymorphisms (SNPs). It utilizes selected restriction enzyme digestion fragments, quick acting ligases that are neutral with the restriction enzyme buffer eliminating buffer exchange steps, and adapters designed to be compatible with Illumina index primers. Library amplification and barcoding are completed in one PCR step, and magnetic beads are used to purify the genomic fragments from the ligation and library generation steps. Our protocol increases the efficiency and decreases the time to complete a ddRAD-Seq experiment. To demonstrate its utility, we compared SNPs from our protocol with those from whole genome resequencing data from Gossypium herbaceum and Gossypium arboreum. Principal component analysis demonstrated that the variability of the combined data was explained by the genotype (PC1) and methodology applied (PC2). Phylogenetic analysis showed that the SNPs from our method clustered with SNPs from the resequencing data of the corresponding genotype. Sequence alignments illustrated that for homozygous loci, more than 90% of the SNPs from the resequencing data were discovered by our method. Our analyses suggest that our ddRAD-Seq method is reliable in identifying SNPs suitable for phylogenetic and association genetic studies while reducing cost and time over known methods.

Dual Domestication, Diversity, and Differential Introgression in Old World Cotton Diploids.

Domestication in the cotton genus is remarkable in that it has occurred independently four different times at two different ploidy levels. Relatively little is known about genome evolution and domestication in the cultivated diploid species Gossypium herbaceum and Gossypium arboreum, due to the absence of wild representatives for the latter species, their ancient domestication, and their joint history of human-mediated dispersal and interspecific gene flow. Using in-depth resequencing of a broad sampling from both species, we provide support for their independent domestication, as opposed to a progenitor-derivative relationship, showing that diversity (mean π = 6 × 10-3) within species is similar, and that divergence between species is modest (FST = 0.413). Individual accessions were homozygous for ancestral single-nucleotide polymorphisms at over half of variable sites, while fixed, derived sites were at modest frequencies. Notably, two chromosomes with a paucity of fixed, derived sites (i.e., chromosomes 7 and 10) were also strongly implicated as having experienced high levels of introgression. Collectively, these data demonstrate variable permeability to introgression among chromosomes, which we propose is due to divergent selection under domestication and/or the phenomenon of F2 breakdown in interspecific crosses. Our analyses provide insight into the evolutionary forces that shape diversity and divergence in the diploid cultivated species and establish a foundation for understanding the contribution of introgression and/or strong parallel selection to the extensive morphological similarities shared between species.

Complete Genome Sequences of Four Campylobacter jejuni Strains Isolated from Retail Chicken Meat and Broiler Feces.

Campylobacter jejuni is the leading pathogen that causes foodborne infections. Here, we report the complete genome sequences of four C. jejuni strains isolated from retail chicken meat and broiler feces samples. Genes encoding type VI secretion and antibiotic resistance were detected among these isolates.

Transcriptome analysis of the 2,4-dichlorophenoxyacetic acid (2,4-D)-tolerant cotton chromosome substitution line CS-B15sh and its susceptible parental lines G. hirsutum L. cv. Texas Marker-1 and G. barbadense L. cv. Pima 379.

The cotton chromosome substitution line, CS-B15sh, exhibits 41% lower injury from 2,4-D when applied at the field recommended rate of 1.12 kg ae ha-1 (1×) than does Texas Marker-1 (TM-1). CS-B15sh was developed in the genetic background of Gossypium hirsutum L. cv TM-1 and has chromosome introgression on the short arm of chromosome 15 from Gossypium barbadense L. cv. Pima 379. In a previous experiment, we observed reduced translocation of [14C]2,4-D outside the treated leaf tissue in CS-B15sh, which contrasted with an increased translocation of the herbicide in the tissues above and below the treated leaf in TM-1. Our results indicate a potential 2,4-D tolerance mechanism in CS-B15sh involving altered movement of 2,4-D. Here, we used RNA sequencing (RNA-seq) to determine the differential expression of genes between 2,4-D-challenged and control plants of the tolerant (CS-B15sh) and susceptible lines (TM-1 and Pima 379). Several components of the 2,4-D/auxin-response pathway-including ubiquitin E3 ligase, PB1|AUX/IAA, ARF transcription factors, and F-box proteins of the SCFTIR1/AFB complex-were upregulated with at least threefold higher expression in TM-1 compared with CS-B15sh, while both Pima 379 and TM-1 showed the same fold change expression for PB1|AUX/IAA mRNA. Some genes associated with herbicide metabolism, including flavin monooxygenase (Gohir.A01G174100) and FAD-linked oxidase (Gohir.D06G002600), exhibited at least a twofold increase in CS-B15sh than in TM-1 (the gene was not expressed in Pima 379), suggesting a potential relationship between the gene's expression and 2,4-D tolerance. It is interesting to note that glutathione S-transferase was differentially expressed in both CS-B15sh and Pima 379 but not in TM-1, while cytochrome P450 and other genes involved in the oxidation-reduction process were significantly expressed only in CS-B15sh in response to 2,4-D. Gene set enrichment analysis on the union DEGs of the three cotton genotypes revealed the depletion of transcripts involved in photosynthesis and enrichment of transcripts involved in ABA response and signaling.

Male fathead minnow transcriptomes and associated chemical analytes in the Milwaukee estuary system.

Contaminants of Emerging Concern (CECs) can be measured in waters across the United States, including the tributaries of the Great Lakes. The extent to which these contaminants affect gene expression in aquatic wildlife is unclear. This dataset presents the full hepatic transcriptomes of laboratory-reared fathead minnows (Pimephales promelas) caged at multiple sites within the Milwaukee Estuary Area of Concern and control sites. Following 4 days of in situ exposure, liver tissue was removed from males at each site for RNA extraction and sequencing, yielding a total of 116 samples from which libraries were prepared, pooled, and sequenced. For each exposure site, 179 chemical analytes were also assessed. These data were created with the intention of inviting research on possible transcriptomic changes observed in aquatic species exposed to CECs. Access to both full sequencing reads of animal samples as well as water contaminant data across multiple Great Lakes sites will allow others to explore the health of these ecosystems in support of the aims of the Great Lakes Restoration Initiative.

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.

Transcriptomic Analysis of Listeria monocytogenes in Response to Bile Under Aerobic and Anaerobic Conditions.

Listeria monocytogenes is a gram-positive facultative anaerobic bacterium that causes the foodborne illness listeriosis. The pathogenesis of this bacterium depends on its survival in anaerobic, acidic, and bile conditions encountered throughout the gastrointestinal (GI) tract. This transcriptomics study was conducted to analyze the differences in transcript levels produced under conditions mimicking the GI tract. Changes in transcript levels were analyzed using RNA isolated from L. monocytogenes strain F2365 at both aerobic and anaerobic conditions, upon exposure to 0 and 1% bile at acidic and neutral pH. Transcripts corresponding to genes responsible for pathogenesis, cell wall associated proteins, DNA repair, transcription factors, and stress responses had variations in levels under the conditions tested. Upon exposure to anaerobiosis in acidic conditions, there were variations in the transcript levels for the virulence factors internalins, listeriolysin O, etc., as well as many histidine sensory kinases. These data indicate that the response to anaerobiosis differentially influences the transcription of several genes related to the survival of L. monocytogenes under acidic and bile conditions. Though further research is needed to decipher the role of oxygen in pathogenesis of L. monocytogenes, these data provide comprehensive information on how this pathogen responds to the GI tract.

Genome assembly of two nematode-resistant cotton lines (Gossypium hirsutum L.).

Upland cotton (Gossypium hirsutum L.) is susceptible to damage by the root-knot and the reniform nematodes, causing yield losses greater than 4% annually in the United States. In addition, these nematodes are synergistic with seeding disease and root rot pathogens that exacerbate diseases and subsequent yield losses. Production practices to minimize nematode damage include crop rotation and nematicides, but these techniques need to be repeated and are expensive. The use of resistant cultivars is deemed the most effective and economical approach for managing nematodes in cotton. Here, we describe the genomes of two nematode-resistant lines of cotton, BARBREN-713 and BAR 32-30. These genomes may expedite the development of DNA markers that can be used to efficiently introduce nematode resistance into commercially valuable Upland lines.

The Gossypium anomalum genome as a resource for cotton improvement and evolutionary analysis of hybrid incompatibility.

Cotton is an important crop that has been the beneficiary of multiple genome sequencing efforts, including diverse representatives of wild species for germplasm development. Gossypium anomalum is a wild African diploid species that harbors stress-resistance and fiber-related traits with potential application to modern breeding efforts. In addition, this species is a natural source of cytoplasmic male sterility and a resource for understanding hybrid lethality in the genus. Here, we report a high-quality de novo genome assembly for G. anomalum and characterize this genome relative to existing genome sequences in cotton. In addition, we use the synthetic allopolyploids 2(A2D1) and 2(A2D3) to discover regions in the G. anomalum genome potentially involved in hybrid lethality, a possibility enabled by introgression of regions homologous to the D3 (Gossypium davidsonii) lethality loci into the synthetic 2(A2D3) allopolyploid.

TALEN-Based HvMPK3 Knock-Out Attenuates Proteome and Root Hair Phenotypic Responses to flg22 in Barley.

Mitogen activated protein kinases (MAPKs) integrate elicitor perception with both early and late responses associated with plant defense and innate immunity. Much of the existing knowledge on the role of plant MAPKs in defense mechanisms against microbes stems from extensive research in the model plant Arabidopsis thaliana. In the present study, we investigated the involvement of barley (Hordeum vulgare) MPK3 in response to flagellin peptide flg22, a well-known bacterial elicitor. Using differential proteomic analysis we show that TALEN-induced MPK3 knock-out lines of barley (HvMPK3 KO) exhibit constitutive downregulation of defense related proteins such as PR proteins belonging to thaumatin family and chitinases. Further analyses showed that the same protein families were less prone to flg22 elicitation in HvMPK3 KO plants compared to wild types. These results were supported and validated by chitinase activity analyses and immunoblotting for HSP70. In addition, differential proteomes correlated with root hair phenotypes and suggested tolerance of HvMPK3 KO lines to flg22. In conclusion, our study points to the specific role of HvMPK3 in molecular and root hair phenotypic responses of barley to flg22.

The Gossypium stocksii genome as a novel resource for cotton improvement.

Cotton is an important textile crop whose gains in production over the last century have been challenged by various diseases. Because many modern cultivars are susceptible to several pests and pathogens, breeding efforts have included attempts to introgress wild, naturally resistant germplasm into elite lines. Gossypium stocksii is a wild cotton species native to Africa, which is part of a clade of vastly understudied species. Most of what is known about this species comes from pest resistance surveys and/or breeding efforts, which suggests that G. stocksii could be a valuable reservoir of natural pest resistance. Here, we present a high-quality de novo genome sequence for G. stocksii. We compare the G. stocksii genome with resequencing data from a closely related, understudied species (Gossypium somalense) to generate insight into the relatedness of these cotton species. Finally, we discuss the utility of the G. stocksii genome for understanding pest resistance in cotton, particularly resistance to cotton leaf curl virus.

Transcriptomic response patterns of hornyhead turbot (Pleuronichthys verticalis) dosed with polychlorinated biphenyls and polybrominated diphenyl ethers.

To evaluate the impact of environmental contaminants on aquatic health, extensive surveys of fish populations have been conducted using bioaccumulation as an indicator of impairment. While these studies have reported mixtures of chemicals in fish tissues, the relationship between specific contaminants and observed adverse impacts remains poorly understood. The present study aimed to characterize the toxicological responses induced by persistent organic pollutants in wild-caught hornyhead turbot (P. verticalis). To do so, hornyhead turbot were interperitoneally injected with a single dose of PCB or PBDE congeners prepared using environmentally realistic mixture proportions. After 96-hour exposure, the livers were excised and analyzed using transcriptomic approaches and analytical chemistry. Concentrations of PCBs and PBDEs measured in the livers indicated clear differences across treatments, and congener profiles closely mirrored our expectations. Distinct gene profiles were characterized for PCB and PBDE exposed fish, with significant differences observed in the expression of genes associated with immune responses, endocrine-related functions, and lipid metabolism. Our findings highlight the key role that transcriptomics can play in monitoring programs to assess chemical-induced toxicity in heterogeneous group of fish (mixed gender and life stage) as is typically found during field surveys. Altogether, the present study provides further evidence of the potential of transcriptomic tools to improve aquatic health assessment and identify causative agents.

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.