Gene

BackgroundDUXC is a multi-copy transcription factor gene found within a long tandem repeat locus in several Laurasiatherians. It is suggested to be functionally similar to human DUX4 because of its shared C-terminal domain and its close phylogenetic relationship to DUX4. DUX family genes are transiently expressed in preimplantation embryos of placental mammals. However, early embryo-derived cDNA proof for DUXC, which is needed for its further functional characterization, has not been reported so far. ResultsOur study provides a full-length sequence of DUXC mRNA, derived from the 8-cell stage in vitro fertilization (IVF) bovine embryos, containing double homeobox and 9aa transactivation domain (9aaTAD)-encoding sequences. Identified DUXC sequence uncovered a first exon that was not previously annotated. We showed that DUXC mRNA levels are independent of the embryonic transcription at the 2-, 4-, and 8-cell stage, whereas its decline, observed from the 8-cell stage onwards, is minor embryonic genome activation (EGA)-dependent. We also investigated the genomic organisation of the DUXC array in eight different cattle breed assemblies, revealing polymorphic internal repeats flanked by an incomplete distal unit at the telomeric end and a much shorter unit at the proximal end of the DUXC array. Despite the presence of a putative polyadenylation signal downstream the distal unit, we presented evidence for the expression of internal but not distal DUXC in early bovine IVF embryos. ConclusionsDUXC is a potential bovine EGA inducer, supported by its expression at peak levels at pre-EGA stages followed by a decrease with a dependency on minor EGA.

Short Interrupted Repeats Cassettes (SIRC) are recently discovered eukaryotic DNA elements possessing many traits of satellite DNA and mobile genetic elements, and consisted of short direct repeats interspersed with diverse spacer sequences. The SIRC ensemble of individual species is highly heterogenous and cannot be studied using alignment methods. It was found that number of similar SIRC sequences in a given pair of species is in general correlated with their taxonomic distance, and, at the same time, closely related species can possess very diverged SIRC ensembles, which makes SIRC evolutionary pattern closer to mobile genetic element type. The SIRC sequences make up clusters with comparable sequence patterns, that are likely to demonstrate doublet evolutionary model which strongly supports that the SIRC structure is supported by the evolutionary selection. Several SIRC sequences of Arabidopsis were found to be of ancient origin with traceable evolution history as far as to the moss clade. We carried out unbiased detection of SIRC ensembles in 10 plant genomes and found that, despite very high intraspecies heterogeneity, SIRC sets possess strong interspecies phylogenetic signal. Key messageShort Interrupted Repeats Cassettes are elements of ancient origin, and could potentially be used to trace organism history, and to facilitate syntheny and Hi-C analysis.

The cockroach Blattella germanica possesses panoistic ovaries, in which oocytes lack nurse cells and therefore need to rely on their own transcriptional activity to support embryogenesis. Ovarian development in this species involves the development of a single basal ovarian follicle (BOF) per gonadotropic cycle, a process strictly regulated by endocrine signals, primarily juvenile hormone and ecdysone, which act at both the transcriptional and translational levels. In addition, transcriptional activity in these ovaries is necessary for both regulating and genome protection, and at this level, PIWI-interacting RNAs (piRNAs) play an essential role. Although insect ovaries are known to be particularly rich in piRNAs, their function in ovary maturation is still not well defined. For this purpose, we characterize the piRNA expression dynamics across seven key developmental and reproductive stages, ranging from late nymphal instars to post-vitellogenic adults. piRNA expression in B. germanica shows coordinated fluctuations. Expression remains stable in previtellogenic ovaries, whereas vitellogenic ovaries show pronounced changes. Moreover, vitellogenic ovaries exhibit reduced piRNA diversity due to strong enrichment of a subset of highly expressed piRNAs. Our data show that although piRNAs predominantly map to transposable elements, particularly LINEs, there is a notable increase in gene-derived piRNAs toward the end of the cycle. Our results suggest regulatory roles of piRNAs in modulating both TEs and mRNAs during BOF maturation, likely related to changes in the follicular cell program.

Vitamin D signaling has recognized roles in female reproductive physiology, but its effects at the chromatin level in endometrial stromal cells are still unclear. Here, we investigated how the active form of vitamin D, 1,25-dihydroxyvitamin D3, or calcitriol, influences the accessible chromatin landscape of human endometrial stromal cells. Assay for transposase-accessible chromatin using sequencing (ATAC-seq) was performed on T-HESCs treated with either a vehicle or 1,25(OH)2D3. Ligand treatment increased overall chromatin accessibility, shown by higher ATAC-seq signal intensity, while causing only minor changes in the total number of called peaks. Peak annotation revealed that accessible regions were spread across both promoter-proximal and distal genomic areas. Integrating this data with CUT&RUN and RNA sequencing showed that most vitamin D-responsive cistromic modifications and transcripts were linked to nearby open chromatin, though fewer were associated with regions that were significantly differentially accessible. These results suggest that 1,25(OH)2D3-dependent transcription mainly occurs within a permissive, pre-accessible chromatin environment. This study offers new evidence that active vitamin D influences the epigenomic landscape of human endometrial stromal cells, establishing the chromatin-based molecular response to a chemically-defined VDR ligand, 1,25(OH)2D3, relevant to stromal differentiation and preparation for decidualization. HighlightsO_LIFirst evidence suggesting the direct impact of active vitamin D, 1,25-dihydroxyvitamin D3, 1,25(OH)2D3, enhanced the signal intensity of chromatin accessibility in human endometrial stromal cells C_LIO_LIMost accessible chromatin regions were shared between vehicle and ligand-treated human endometrial stromal cells C_LIO_LI1,25(OH)2D3-responsive transcription occurs largely within pre-accessible chromatin in human endometrial stromal cells C_LIO_LIAssay for transposase-accessible chromatin sequencing (ATAC-seq) defines a chromatin-level pharmacologic response to a chemically defined VDR ligand in human endometrial stromal cells C_LI

Cancer-related genomic instability (GI) may cause genetic alterations in spermatozoa, implying health issues not only in cancer survivors, but also in their children [1, 2]. We therefore studied Loss of Y chromosome (LOY), considered as hallmark of GI [3-15], in spermatozoa and blood from survivors of childhood and testicular cancer (CC, TC), and controls (CTRL). We found that LOY was statistically significantly more frequent in spermatozoa from cancer survivors than in controls (Odds Ratio [OR]=2.2 for CC vs. CTRL and OR=2.4 for TC vs. CTRL). Furthermore, LOY was about an order of magnitude more prevalent in spermatozoa than in blood among 18-53-year-old males within all cohorts. Our findings suggest that LOY in spermatozoa might be a clinically useful marker of GI, reduced fertility and disease predisposition in males. Introducing LOY in spermatozoa as a biomarker opens a new research avenue into disease prevention and the causes and consequences of LOY.

Photoperiod sensitivity (PS) is a key biological response in plants as they adapt to specific environments. Rice (Oryza sativa L.) exhibits a clear PS, as it implements critical phase transition decisions based on PS signals. In this study, we identified a novel PS gene, JMJ706, that is expected to deliver photoperiod-related signals to the flowering-time regulatory network in a day-length-dependent manner. The JMJ706 mutants exhibit early flowering under LD and later flowering under SD compared to WT plants. The gene encodes an H3K9me2 demethylase, and under long-day (LD) conditions, its demethylase activity facilitates the expression of Grain number, Plant height, and Heading-date7 (Ghd7). Since Ghd7 is a floral repressor in LD, it promotes the vegetative phase by delaying flowering. Under short-day conditions (SD), H3K9me2 demethylase activity facilitates Early heading-date 1 (Ehd1) expression, and it acts as a floral accelerator by inducing Heading date 3 (Hd3a) and RICE FLOWERING LOCUS T 1 (RFT1). Furthermore, we propose that the daylength-dependent promotion of target genes (Ghd7 and Ehd1) occurs through demethylation of specific promoter regions at a crucial time window. In addition, JMJ706 may play an important role in regulating plant architecture, including plant height. The natural variation in JMJ706 alleles shows high frequencies across major rice subpopulations, suggesting that JMJ706 could play an important role in the geographical distribution and adaptation of rice cultivars. Our results may add a new layer to the rice flowering-time regulatory pathway, supporting regional adaptation and potential for future breeding.

A greater genetic susceptibility has been proposed as an explanation of the greater rates of cardiovascular and metabolic disease in South Asian relative to European populations. We first demonstrate that after accounting for technical artefacts the genetic effects for related traits are largely consistent between ancestral groups, which downplays the role of GxG or GxE interactions driving differential prevalence. If higher genetic susceptibility in South Asians is due to selective pressures acting through adiposity-related traits in the evolutionary past, signatures of selection should be evident at loci associated with cardiometabolic disease and other causally related traits (e.g. fat distribution). We tested for enrichment of several selection statistics (FST, XP-EHH and XP-nSL) at loci associated with a range of traits related to cardiometabolic disease, in comparison to a null distribution of linkage disequilibrium (LD) score and minor allele frequency (MAF) matched SNPs. Loci associated with a subset of these traits (Type 2 diabetes mellitus, trunk fat percentage, body fat percentage and trunk fat mass) exhibited enrichment for FST, consistent with a moderate adaptive explanation for their cross-population differentiation. In contrast, none of the studied traits were enriched for haplotype-based statistics, indicative that cross population genetic divergence is unlikely to have been driven by recent selective sweeps but has rather likely arisen from either ancient selection or recent polygenic selection acting on standing variation.

In wheat, weak seed dormancy (SD) is related to an increased tendency for pre-harvest sprouting (PHS), which reduces yield and quality. However, the molecular mechanism underlying SD remains elusive. Here, we identified a wheat R2R3-MYB transcription factor (TaMYB83-7B) related to SD. Expression analysis showed that TaMYB83-7B was highly expressed in wheat seeds, and was more highly expressed in strong-dormancy varieties than in weak-dormancy varieties. Sequence and association analysis indicated that T/C mutations at -907 bp and -1133 bp in the TaMYB83-7B promoter were significantly associated with wheat SD, with C at both sites related to strong dormancy. Dual-luciferase reporter assays demonstrated that the transcriptional activity of the TaMYB83-7B promoter was significantly higher in strong-dormancy varieties than in weak-dormancy varieties. Further analyses indicated that TaMYB83-7B functions as a transcriptional inhibitor. Germination experiments revealed that overexpression of TaMYB83-7B significantly enhanced SD, while its loss-of-function reduced SD. Finally, TaMYB83-7B was found to regulate SD by influencing the balance between abscisic acid (ABA) and gibberellin (GA) in wheat seeds. Overall, the results of this study enhance our understanding of the complex regulatory mechanism underlying SD, and provide gene targets and molecular markers for the genetic improvement of PHS resistance in wheat.

The only autonomously active transposable elements in the human genome are Long interspersed nuclear element-1 (LINE-1) elements. These elements are known to play an important role in changing the transcriptome. LINE-1 sequences affect gene regulation during post-transcription processing, along with their established role in retrotransposition. Exonization is one mechanism where the LINE-1 integrated genome undergoes alternative splicing to produce new isoforms of transcripts. Our work mainly highlights the effect of LINE-1 associated exonization, focusing on the formation of isoforms of transcripts. Using Non-small cell lung cancer (NSCLC) as a model, we conducted a detailed transcriptome study that combines splice junction profiling with gene expression data. Our results show that LINE-1 sequences are often included as exons in host transcripts, leading to the formation of new exons and their various isoforms. The events are validated by solid splice junction evidence that proves the reliability and reproducibility. In particular, it was observed that repetitive analyses revealed certain LINE-1 exonization events that were consistent. The finding indicates that LINE-1 act as recurrent sources of splice ready sequences. Though exonizations do not necessarily affect the total expression levels of genes, our study reveals that they certainly contribute to transcript diversity. The diversity of isoforms generated potentially contributes to the effects of gene function. This study is limited to NSCLC, but it is likely that the exonizations events play a crucial role in the altering RNA diversity in cancers. Therefore the study elucidates new insights into how transposable elements modify gene structure and function during cancer development.

BackgroundConventional models of human ribosomal DNA (rDNA) array organization have historically depended on transcription-centric boundaries, partitioning the unit into a [~]13 kb rDNA transcription region and a monolithic [~]31 kb intergenic spacer (IGS). While our previous identification of Duplication Segment Units (DSUs) mapped these arrays based on an intuitive analysis of the microsatellite density landscape of the complete reference human genome, our present deep mining of this landscape has revealed a more accurate rDNA Gene Unit Pattern. Methods & ResultsIn this study, we conducted a deep mining analysis of our previously established microsatellite density landscape of the T2T-CHM13 assembly, focusing specifically on nucleolar organizing regions (NORs). We suggest a more accurate rDNA Gene Unit Pattern containing a (CTTT)n microsatellite aggregation ahead of the rDNA gene and a (CT)n microsatellite aggregation behind the gene, rather than a pattern featuring an IGS region inserted between two rDNA genes. ConclusionsA correct rDNA gene pattern of the human genome probably includes a (CTTT)n microsatellite aggregation ahead of the gene and a (CT)n microsatellite aggregation behind it, which possibly constitute cis- and trans-regulating regions; the (CTTT)n and (CT)n microsatellite aggregations may provide two different local stable DNA structures for regulatory protein binding.

Some Anopheles species that act as malaria vectors are members of species complexes, a concept whereby sibling species cannot be differentiated solely on the basis of morphological characters. Therefore, species complexes represent a major problem in malaria vector control, because within an Anopheles complex, vectors cannot be differentiated from non-vector species, unless molecular techniques are used to identify them. The Anopheles tessellatus species complex is an important potential vector in South, East, and Southeast Asia, including certain regions of Indonesia. However, no in-depth studies have been conducted on this species complex in that country. Therefore, this study investigated the taxonomic status of An. tessellatus from diverse populations across five Indonesian islands (Sumatra, Java, West Nusa Tenggara, East Nusa Tenggara, and Sulawesi) and identified interpopulation genetic variation based on molecular data of the ITS2, COI, and COII genes. Phylogenetic relationships were constructed using the Maximum Likelihood method. Haplotype and network analysis were also conducted. The results indicate that An. tessellatus constitutes a monophyletic group comprising three well-defined lineages that exhibit clear intraspecific genetic differentiation. Cluster 1 corresponds to the population of Sumatra, Cluster 2 represents population from Sulawesi, and Cluster 3 encompasses populations from Java, West Nusa Tenggara, and East Nusa Tenggara. These findings demonstrate high haplotype diversity and low nucleotide diversity within the species. Populations from West Sumatra, Manado, Tojo Una - Una, and North Morowali (Sulawesi) have the potential for speciation with a genetic distance of 0.61 - 0.94% for COI, between 0.81 - 0.95% for ITS2, and between 0.62 - 0.71% for COII. These findings underscore the need for further integrative studies to obtain a more comprehensive understanding of the An. tessellatus complex in Indonesia and its role in malaria transmission.

MicroRNAs (miRNAs) play essential roles in the posttranscriptional regulation of gene expression in organisms. In the process of synthesizing mature miRNAs from miRNA precursors, the miRNA precursors are cleaved via Dicer at their loop structure, after which the miRNA precursors become mature and regulate transcription. However, the consequences of altering the loop sequence are not fully understood. The silkworm Bombyx mori is a lepidopteran insect with many genetic strains. We identified a mutant of the miRNA miR-3260 whose the part of the loop structure was lacking in a silkworm strain with translucent larval skin. Here, we aimed to analyze the role of wild-type miR-3260 and the influence of the mutation of the loop structure in B. mori. First, we identified the genomic region responsible for the translucent larval skin phenotype and determined that the mutated miR-3260 nucleotide sequences. Then, we predicted the binding partners of wild-type miR-3260 using the RNA hybrid tool and found two juvenile hormone (JH)-related genes as targets of wild-type miR-3260. Next, we assessed the relationships between miR-3260 and JH and found that miR-3260 was highly expressed in the Corpora allata and its expression responded to JH treatment. Meanwhile, miR-3260 mimic and inhibitor did not induce the typical phenotypes associated with JH in B. mori. Then, we compared the dicing products from wild-type and mutant miR-3260 precursors and observed that neither form underwent Dicer-mediated cleavage when the loop structure was altered. These results suggest that loop mutations in the miR-3260 precursor may not influence dicing activity, consistent with the lack of observable phenotypic effects.

Xenopus is a genus of entirely aquatic frogs found in sub-Saharan Africa. Currently, the complete genomes of two species within the Xenopus genus, Xenopus laevis and Xenopus tropicalis, have been fully sequenced, annotated, and made publicly available. The two species inhabit markedly different environments: X. tropicalis lives in the hot, equatorial regions of Africa, whereas X. laevis resides in the cooler climates of southern Africa. In the present study, mutational profiling, comparative homology modeling, and computational bioinformatics were used to identify the features of adaptive evolution in Xenopus endonuclease G (EndoG) proteins. The multiple characteristics of EndoG isozymes were discovered to vary considerably between the two Xenopus species dwelling in different locations. Most notably, EndoG proteins from the psychrophilic X. laevis exhibit the increased contents of charged and polar residues, elevated pI, higher intramolecular interaction energies, B factors, molecular void volumes, and solvent accessibilities, but the decreased contents of nonpolar and aromatic amino acids, lower hydrophobicity, buried surface area, and molecular packing density compared to those from the thermophilic X. tropicalis. The observed differences strongly suggest that temperature plays a dominant role in EndoG diversification. Evaluation of intramolecular interaction energies appears to be a particularly sensitive and discriminative framework for assessing protein divergence at the structural level. Overall, this study highlights the diversification of homologous proteins in ectothermic vertebrate eukaryotes and provides mechanistic insight into protein adaptation to contrasting environments.

BACKGROUNDDuring spermiogenesis, histones are exchanged by protamines (PRMs) in spermatids, which results in DNA hypercondensation and protection. Rodents and primates express two PRMs (PRM1 and PRM2) in a species-specific ratio. Maintaining this ratio is necessary for functional chromatin reorganization and alteration is associated with sub- or infertility in mice and humans. Prm1 and Prm2 deficient mice are infertile, while Prm1+/- males are subfertile showing a severely altered PRM ratio. Prm2+/- males are fertile and display a protamine ratio comparable to WT. OBJECTIVESHere, we addressed the question whether loss of one allele of Prm1 and one allele of Prm2 affects fertility. MATERIAL AND METHODSDouble heterozygous (dHET) mice lacking one allele of Prm1 and one allele of Prm2 were generated and analyzed RESULTSdHET males were infertile with sperm showing retention of histones and TNPs, high levels of PRM2 precursor and decreased levels of mature PRM2. In mature sperm the PRM ratio and the total PRM content was not altered. However, CMA3 staining revealed incomplete protamination and sperm nuclei appeared more rounded and slightly bigger, suggesting impaired DNA-hypercondensation. In dHET sperm, DNA degradation was apparent, but to a lower level compared to sperm from Prm1 and Prm2 deficient males. Increased 8-OHdG levels suggested oxidative stress in the epididymis of dHET mice. However, a fraction of dHET sperm were capable of fertilization, with embryonic development up to 8-cell stage. DISCUSSION AND CONCLUSIONThese results suggest, that male factor infertility might not be reliably detected by measuring PRM1/PRM2 ratio but rather by determining the level of protamination by e.g. CMA3 analysis and pre-PRM2 retention.

The paper presents the results of a comparative analysis of gene networks activated by water stress and low temperatures in extensive (Saratovskaya 29, S29) and intensive (Yanetskis Probat, YP) wheat varieties during the seedling development stage. It is concluded that the creation of the S29 variety, which occurred through complex stepwise hybridization and selection for morphological traits, productivity, and grain quality traits, resulted in the emergence and inheritance of regulatory gene networks involving proteins with the CC domain, as well as the BTB/POZ and TAZ domains, which have an increased affinity for transcription factors involved in the transcriptomic response to changing external conditions. It was established that, at the transcriptomic level, the S29 variety is characterized by a transition to an energy saving mode to maintain the activity of the Calvin-Benson cycle under the water deficit conditions and the inhibition of proteolytic processes at low temperatures. The transcriptional response of the high-yielding YP variety to 24-hour low-temperature treatment was more active and involved a larger number of genes compared to the S29 variety. Identifying varietal variability in molecular genetic mechanisms of resistance to abiotic stressors facilitates the development of marker-assisted and genomic selection technologies for common wheat. Key messageThe extensive S29 variety was characterized by its transition to energy-saving mode to maintain the Calvin-Benson cycle under water deficit and a reduction in proteolytic processes under low temperature.

Nucleobindin 1 (NUCB1) is a multifunctional conserved protein located in Golgi luminal, nucleus, extracellular and cytosolic pools. NUCB1 is multidomain protein comprised of a signal peptide, a DNA-binding domain, a leucine zipper and Ca2+ -binding domain. The multiple domains and localization of NUCB1 potentiates its interactions with various partners, such as DNA, Gi3 protein, cyclooxygenase 2, LRP10 and RNA suggests its importance in the regulation of many cellular events. We revealed that NUCB1 contains three RNA-binding regions and able to interact with two RNA fragments. It was suggested possible variants of the participation of NUCB1 in the interaction of the two partially complementary RNAs. The RNA-binding properties of the NUCB1 were also confirmed in vivo experiments.

In mammals, silent state of one of the X-chromosomes in female balance the X-dosage between sexes. In parallel, X to autosome imbalance due to monoallelic expression of X-linked genes relative to biallelic autosomal genes, is primarily compensated through the X-chromosome upregulation (XCU). It has been demonstrated that X-chromosome inactivation (XCI) and XCU coincides during embryogenesis, however, XCU is not global and it occurs in a gene-specific manner. The underlying mechanistic aspects of such specificity of XCU remain unknown. Here, we provide systematic and comparative analysis across eutherians (mouse, human and pig) and metatherian (opossum) embryogenesis to determine if evolutionary origins shape the XCU. Intriguingly, we show that while evolutionary older X-linked genes (predating mammalian divergence) undergo robust XCU consistently across developmental stages, younger mammalianLJorigin genes do not. Similarly, the eutherian X-linked genes conserved in metatherian X (X-ortholog) undergo robust XCU, whereas genes orthologus to metatherian autosome (Auto-ortholog) exhibit weaker pattern of XCU. Further, strata-wise comparison revealed that genes in older XLJchromosome strata (1-2) consistently undergo upregulation, whereas strata 3-4 genes do not. Importantly, we show that different evolutionary classes of X-linked genes, which undergo robust XCU, are often enriched with active chromatin marks (H3K36me3, H3K27ac and H3K4me1) relative to the autosome, suggesting that chromatin state mediate the XCU. Moreover, we show that often active-marks enrichment correlates with differential XCU dynamics of different class of genes. Taken together, our study provides significant insight into the evolutionary dynamics of XCU and underlying mechanistic framework.

In many oviparous animals, egg yolk is the sole source of nutrition until feeding begins, and carbohydrates are present in only small amounts in the yolk. Glucose plays an important role in the developmental processes of various animals. In addition, gluconeogenesis has been reported to occur in the yolk syncytial layer (YSL) of cartilaginous fish and teleosts. In contrast, the role of gluconeogenesis in tetrapods remains unclear. In this study, we used Xenopus tropicalis, an anuran amphibian, which lacks YSL, and therefore provide an opportunity to examine the evolutionary conservation of gluconeogenic mechanisms among vertebrates. In X. tropicalis, liquid chromatography/mass spectrometry revealed that glucose levels increased before liver formation. Subsequent tracer experiments using 13C-labeled metabolic substrates detected gluconeogenesis activity from glycerol and lactate. Expression analyses showed that gluconeogenic genes are expressed in the epidermis and endoderm. Consistently, G0 knockout of fbp1, a key gluconeogenic gene, resulted in a significant reduction in glucose levels, affecting brain development. These findings first demonstrate that gluconeogenesis supports development of X. tropicalis. To the best of our knowledge, gluconeogenesis in developing epidermis has not been reported, highlighting previously unrecognized diversity in tissue-specific metabolism during vertebrate development. Comparative analyses across species will provide further insights into the evolution and functional significance of embryonic gluconeogenesis and nutrient metabolism.

Characterizing genomic properties such as genome size, ploidy level, heterozygosity, and repetitive DNA proportion and composition without relying on genome assembly is crucial for profiling the genomes of non-model species. Little is known about the nuclear genome of the large neotropical orchid genus Epidendrum. This study compares genome profiles of Epidendrum anisatum and Epidendrum marmoratum, using flow cytometry and k-mer analysis approaches, as well as bioinformatics ploidy level estimation and repeatome characterization. Multiple depths of coverage, k values, and k-mer-based tools for genome size estimation were explored and contrasted with cytometry genome size estimations. Cytometry and k-mer analyses yielded a consistently higher genome size for E. anisatum (mean 1C genome size = 2.59 Gb) than E. marmoratum (mean 1C genome size = 1.13 Gb), which represents a 2.3-fold genome size difference. Both species were identified as diploid with no evidence of strict partial endoreplication. The most important aspects to be taken into account to improve genome size estimation were heterozygosity, depth of coverage, and the maximum k-mer coverage. The genomes of both species were found to be highly repetitive (63-73%) and heavily dominated by Ty3-gypsy retrotransposons, particularly those of the Ogre family. Additionally, the genome of E. anisatum was characterized by the presence of a 172 bp satellite (AniS1), which represented 11% of the genome size. Together, both Ty3-gypsy transposons and AniS1 shape the genome size difference between the two genomes. This study provides the first genome profiling for species in the genus Epidendrum, but also highlights the importance of using flow cytometry, cytogenetic approaches and bioinformatics techniques in combination for genome profiling.

We examined the overlap in the genes associated with daily rhythms and with behavioral plasticity in ants. We first investigated the daily rhythms of gene expression in the harvester ant, Pogonomyrmex barbatus, and how the rhythmic genes overlap with others previously shown to be associated with plasticity of foraging behavior. Then, to consider whether the overlap is conserved across ant species, we compared rhythms of gene expression in the diurnal, desert harvester ants with those previously reported for a distantly related nocturnal, subtropical carpenter ant, Camponotus floridanus. First, daily transcriptomes in P. barbatus showed that most genes were expressed in light-dark (LD) and constantly dark (DD) conditions at about the same levels; only 11 genes showed at least a two-fold change in expression. Network analysis identified eleven modules of P. barbatus genes under LD conditions. Of these 11 clusters, modules C1 and C2 seem to be central nodes of the gene expression network, because they are the most highly connected in LD, and show the strongest preservation in DD vs. LD, and contain core clock gene Period. Only one module, C2, showed significant overlap with P. barbatus genes that have 24h-rhythmic expression in both LD and DD. There was significant overlap between modules C1, C2, C10, C11, and P. barbatus genes found previously to be associated with plasticity in the regulation of foraging activity to manage water loss. A comparison of the daily transcriptome of P. barbatus with that of C. floridanus showed significant overlap of 24h-rhythmic genes in LD. Modules C1 and C2 of P. barbatus also overlap with C. floridanus genes previously shown to differ in expression rhythms in nurses and foragers. In combination, these results indicate that genes linking plasticity of the circadian clock and of behavior may be broadly conserved in ants.