The first ITS1 profiling of honey samples from the Southeast Asian region Lombok, Bali and Banggi Island.

Southern Asian flowers offer honeybees a diversity of nectar. Based on its geographical origin, honey quality varies. Traditional methods are less authentic than DNA-based identification. The origin of honey is determined by pollen, polyphenolic, and macro-microorganisms. In this study, amplicon sequencing targets macro-microorganisms in eDNA using the ITS1 region to explore honey's geographical location and authentication. The variety of honey samples was investigated using ITS1 with Illumina sequencing. For all four honey samples, raw sequence reads showed 979,380 raw ITS1 amplicon reads and 375 ASVs up to the phylum level. The highest total number of 202 ASVs up to phylum level identified Bali honey with 211,189 reads, followed by Banggi honey with 309,207 a total number of 111 ASVs, and Lombok represents only 63 ASVs up to phylum level with several read 458,984. Based on Shannon and Chao1, honey samples from Bali (B2) and (B3) exhibited higher diversity than honey from Lombok (B1) and green honey from Sabah (B4), while the Simpson index showed that Banggi honey (B4) had higher diversity. Honey samples had significant variance in mycobiome taxonomic composition and abundance. Zygosaccharomyces and Aspergillus were the main genera found in Lombok honey, with percentages of 68.81% and 29.76% respectively. Bali honey samples (B2 and B3) were identified as having a significant amount of the genus Aureobasidium, accounting for 40.81% and 25% of the readings, respectively. The microbiome composition of Banggi honey (B4) showed a high presence of Zygosaccharomyces 45.17% and Aureobasidium 35.24%. The ITS1 analysis effectively distinguishes between honey samples of different origins and its potential as a discriminatory tool for honey origin and authentication purposes.

The internal transcribed spacer 1 sequence polymorphism brings updates to tsetse species distribution in the northern Cameroon: Importance in planning efficient vector control.

Vector control remains one of the best strategies to prevent the transmission of trypanosome infections in humans and livestock and, thus, a good way to achieve the elimination of human African trypanosomiasis and animal African trypanosomiasis. A key prerequisite for the success of any vector control strategy is the accurate identification and correct mapping of tsetse species. In this work, we updated the tsetse fly species identification and distribution in many geographical areas in Cameroon. Tsetse flies were captured from six localities in Cameroon, and their species were morphologically identified. Thereafter, DNA was extracted from legs of each tsetse fly and the length polymorphism of internal transcribed spacer-1 (ITS1) region of each fly was investigated using PCR. ITS1 DNA fragments of each tsetse species were sequenced. The sequences obtained were analysed and compared to those available in GenBank. This enabled to confirm/infirm results of the morphologic identification and then, to establish the phylogenetic relationships between tsetse species. Morphologic features allowed to clearly distinguish all the tsetse species captured in the South Region of Cameroon, that is, Glossina palpalis palpalis, G. pallicera, G. caliginea and G. nigrofusca. In the northern area, G. morsitans submorsitans could also be distinguished from G. palpalis palpalis, G. tachinoides and G. fuscipes, but these three later could not be distinguished with routine morphological characters. The ITS1 length polymorphism was high among most of the studied species and allowed to identify the following similar species with a single PCR, that is, G. palpalis palpalis with 241 or 242 bp and G. tachinoides with 221 or 222 bp, G. fuscipes with 236 or 237 bp. We also updated the old distribution of tsetse species in the areas assessed, highlighting the presence of G. palpalis palpalis instead of G. fuscipes in Mbakaou, or in sympatry with G. morsitans submorsitans in Dodeo (northern Cameroon). This study confirms the presence of G. palpalis palpalis in the Adamawa Region of Cameroon. It highlights the limits of using morphological criteria to differentiate some tsetse species. Molecular tools based on the polymorphism of ITS1 of tsetse flies can differentiate tsetse species through a simple PCR before downstream analyses or vector control planning.

Genome composition and GC content influence loci distribution in reduced representation genomic studies.

BACKGROUND: Genomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored. RESULTS: Here, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected. CONCLUSIONS: Our results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.

Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428.

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

The discovery of novel noncoding RNAs in 50 bacterial genomes.

Structured noncoding RNAs (ncRNAs) contribute to many important cellular processes involving chemical catalysis, molecular recognition and gene regulation. Few ncRNA classes are broadly distributed among organisms from all three domains of life, but the list of rarer classes that exhibit surprisingly diverse functions is growing. We previously developed a computational pipeline that enables the near-comprehensive identification of structured ncRNAs expressed from individual bacterial genomes. The regions between protein coding genes are first sorted based on length and the fraction of guanosine and cytidine nucleotides. Long, GC-rich intergenic regions are then examined for sequence and structural similarity to other bacterial genomes. Herein, we describe the implementation of this pipeline on 50 bacterial genomes from varied phyla. More than 4700 candidate intergenic regions with the desired characteristics were identified, which yielded 44 novel riboswitch candidates and numerous other putative ncRNA motifs. Although experimental validation studies have yet to be conducted, this rate of riboswitch candidate discovery is consistent with predictions that many hundreds of novel riboswitch classes remain to be discovered among the bacterial species whose genomes have already been sequenced. Thus, many thousands of additional novel ncRNA classes likely remain to be discovered in the bacterial domain of life.

Intergenic sequences harboring potential enhancer elements contribute to Axenfeld-Rieger syndrome by regulating PITX2.

Recent studies have uncovered that noncoding sequence variants may relate to Axenfeld-Rieger syndrome (ARS), a rare developmental anomaly with genetic heterogeneity. However, how these genomic regions are functionally and structurally associated with ARS is still unclear. In this study, we performed genome-wide linkage analysis and whole-genome sequencing in a Chinese family with ARS and identified a heterozygous deletion of about 570 kb (termed LOH-1) in the intergenic sequence between paired-like homeodomain transcription factor 2 (PITX2) and family with sequence similarity 241 member A. Knockout of LOH-1 homologous sequences caused ARS phenotypes in mice. RNA-Seq and real-time quantitative PCR revealed a significant reduction in Pitx2 gene expression in LOH-1-/- mice, while forkhead box C1 expression remained unchanged. ChIP-Seq and bioinformatics analysis identified a potential enhancer region (LOH-E1) within LOH-1. Deletion of LOH-E1 led to a substantial downregulation of the PITX2 gene. Mechanistically, we found a sequence (hg38 chr4:111,399,594-111,399,691) that is on LOH-E1 could regulate PITX2 by binding to RAD21, a critical component of the cohesin complex. Knockdown of RAD21 resulted in reduced PITX2 expression. Collectively, our findings indicate that a potential enhancer sequence that is within LOH-1 may regulate PITX2 expression remotely through cohesin-mediated loop domains, leading to ARS when absent.

Bartonella in bat flies from the Egyptian fruit bat in the Middle East.

In the family of fruit bats, Pteropodidae Gray, 1821, as in the third most diverse group of bats (Chiroptera), the bacterium of the genus Bartonella was detected in several species as well as in a few species of their insect ectoparasites in some tropical and sub-tropical regions of the Old World. The Egyptian fruit bat, Rousettus aegyptiacus (Geoffroy, 1810), is one of the most widespread fruit bats, occurring between South Africa, Senegal, and Pakistan. In this bat species, Candidatus Bartonella rousetti has been detected in three African populations in Nigeria, Kenya, and Zambia. This fruit bat, however, also occurs in the Palaearctic, an area isolating the species geographically and phylogenetically from the Afrotropical part of its distribution range. We screened the blood-sucking bat flies (family Nycteribiidae) from R. aegyptiacus for the presence of the Bartonella bacteria. A rich material of bat fly Eucampsipoda aegyptia (Macquart, 1850), a monoxenous ectoparasite of the Egyptian fruit bats, was collected at 26 localities in seven countries (Egypt, Iran, Jordan, Lebanon, Oman, United Arab Emirates, and Yemen) of the Middle East in 2007-2013. The DNA isolates from the bat flies were subjected to a three-marker (gltA, ssrA, and intergenic spacer region, ITS) multilocus sequence analysis. Based on the amplification of the fragment of ssrA gene by a real-time PCR, 65 E. aegyptia samples from 19 localities in all seven countries were positive for the bacteria. One to five Bartonella-positive individuals of E. aegyptia were collected per one individual of R. aegyptiacus. An analysis of the ITS and gltA genes indicated the presence of an uncultured Bartonella sp., belonging to the Cand. B. rousetti genogroup, identified from populations of the Egyptian fruit bat in Africa. These results support the hypothesis that Bartonella's diversity corresponds to its host's diversity (and phylogenetic structure). Specific lineages of pathogens are present in specific phylogenetic groups of bats.

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila.

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

Spacer-Programmed Two-Dimensional DNA Origami Assembly.

Two-dimensional (2D) DNA origami assembly represents a powerful approach to the programmable design and construction of advanced 2D materials. Within the context of hybridization-mediated 2D DNA origami assembly, DNA spacers play a pivotal role as essential connectors between sticky-end regions and DNA origami units. Here, we demonstrated that programming the spacer length, which determines the binding radius of DNA origami units, could effectively tune sticky-end hybridization reactions to produce distinct 2D DNA origami arrays. Using DNA-PAINT super-resolution imaging, we unveiled the significant impact of spacer length on the hybridization efficiency of sticky ends for assembling square DNA origami (SDO) units. We also found that the assembly efficiency and pattern diversity of 2D DNA origami assemblies were critically dependent on the spacer length. Remarkably, we realized a near-unity yield of ∼98% for the assembly of SDO trimers and tetramers via this spacer-programmed strategy. At last, we revealed that spacer lengths and thermodynamic fluctuations of SDO are positively correlated, using molecular dynamics simulations. Our study thus paves the way for the precision assembly of DNA nanostructures toward higher complexity.

Identification of 1600 replication origins in S. cerevisiae.

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, 'a site that binds Mcm in G1' might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here, we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least three orders of magnitude. Published data quantifying single-stranded DNA (ssDNA) during S phase revealed replication initiation among the most abundant 1600 of these sites, with replication activity decreasing with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5500 sites. Specifically, these sites: (1) appeared in intergenic nucleosome-free regions flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. We conclude that, if sites at which Mcm double hexamers are loaded can function as replication origins, then DNA replication origins are at least threefold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than widely assumed.

Anopheles arabiensis continues to be the primary vector of Plasmodium falciparum after decades of malaria control in southwestern Ethiopia.

BACKGROUND: Investigating the species distribution and their role in malaria transmission is important as it varies from place to place and is highly needed to design interventions appropriate to the site. The current study aimed to investigate the Anopheles mosquito species distribution and their infection rate in southwestern Ethiopia. METHODS: The study was conducted in 14 malaria-endemic kebeles (the smallest administrative unit), which were situated in eight different malaria-endemic districts and four zones in southwestern Ethiopia. Ten per cent of households in each village were visited to collect adult mosquitoes using Centers for Disease Control and Prevention (CDC) light traps. The larval and pupal collection was done from breeding sites within the villages, and reared to adults. Female mosquitoes were morphologically identified. The head and thorax of adult Anopheles mosquitoes were tested for circumsporozoite proteins (CSPs) using ELISA. At the same time, legs, wings, and abdomen were used to identify sibling species using PCR targeting the rDNA intergenic spacers region for species typing of the Anopheles funestus group and the internal transcribed spacer 2 region genes for Anopheles gambiae complex. RESULTS: A total of 1445 Anopheles mosquitoes comprising eight species were collected. Of 813 An. gambiae complex tested by PCR, 785 (97%) were Anopheles arabiensis, and the remaining 28 (3%) were not amplified. There were 133 An. funestus group captured and tested to identify the species, of which 117 (88%) were positive for Anopheles parensis, and 15 (11%) were not amplified. A single specimen (1%) showed a band with a different base pair length from the known An. funestus group species. Sequencing revealed this was Anopheles sergentii. Among 1399 Anopheles tested for CSPs by ELISA, 5 (0.4%) An. arabiensis were positive for Plasmodium falciparum and a single (0.07%) was positive for Plasmodium vivax. CONCLUSIONS: Anopheles arabiensis continues to play the principal role in malaria transmission despite implementing indoor-based interventions for decades. Sequencing results suggest that An. sergentii was amplified by the An. funestus group primer, producing PCR amplicon size of different length. Therefore, relying solely on amplifying a specific gene of interest in grouping species could be misleading, as different species may share the same gene.

The pan-plastome of Hemerocallis citrina reveals new insights into the genetic diversity and cultivation history of an economically important food plant.

BACKGROUND: Hemerocallis citrina Baroni (Huang hua cai in Chinese) is a perennial herbaceous plant grown for its flower buds that are eaten fresh or dried and is known as the vegetarian three treasures. The nuclear genome of H. citrina has been reported, but the intraspecific variation of the plastome (plastid genome) has not yet been studied. Therefore, the panplastome of this species collected from diverse locations is reported here for the first time. RESULTS: In this study, 65 H. citrina samples were resequenced, de novo assembled, and aligned with the published plastome of H. citrina to resolve the H. citrina panplastome. The sizes of the 65 newly assembled complete plastomes of H. citrina ranged from 156,048 bp to 156,263 bp, and the total GC content ranged from 37.31 to 37.34%. The structure of the complete plastomes showed a typical tetrameric structure, including a large single copy (LSC), a small single copy (SSC), and a pair of inverted repeat regions (IRA and IRB). Many nucleotide variants were identified between plastomes, among which the variants in the intergenic spacer region were the most abundant, with the highest number of variants concentrated in the LSC region. Based on the phylogenetic tree constructed using the ML method, population structure analysis, and principal component analysis (PCA), the panplastome data were subdivided into five genetic clusters. The C5 genetic cluster was mostly represented by samples from Qidong, Hunan Province, while samples from Shanxi and Shaanxi Provinces were classified into the C4 genetic cluster. The greatest genetic diversity was found in the C1 genetic cluster, and the greatest genetic distance between any two clusters was found between the C4 and C5 clusters. CONCLUSION: The resolution of the panplastome and the analysis of the population structure of H. citrina plastomes provide important data for future breeding projects and germplasm preservation.

Evolution of the gene regulatory network of body axis by enhancer hijacking in amphioxus.

A central goal of evolutionary developmental biology is to decipher the evolutionary pattern of gene regulatory networks (GRNs) that control embryonic development, and the mechanism underlying GRNs evolution. The Nodal signaling that governs the body axes of deuterostomes exhibits a conserved GRN orchestrated principally by Nodal, Gdf1/3, and Lefty. Here we show that this GRN has been rewired in cephalochordate amphioxus. We found that while the amphioxus Gdf1/3 ortholog exhibited nearly no embryonic expression, its duplicate Gdf1/3-like, linked to Lefty, was zygotically expressed in a similar pattern as Lefty. Consistent with this, while Gdf1/3-like mutants showed defects in axial development, Gdf1/3 mutants did not. Further transgenic analyses showed that the intergenic region between Gdf1/3-like and Lefty could drive reporter gene expression as that of the two genes. These results indicated that Gdf1/3-like has taken over the axial development role of Gdf1/3 in amphioxus, possibly through hijacking Lefty enhancers. We finally demonstrated that, to compensate for the loss of maternal Gdf1/3 expression, Nodal has become an indispensable maternal factor in amphioxus and its maternal mutants caused axial defects as Gdf1/3-like mutants. We therefore demonstrated a case that the evolution of GRNs could be triggered by enhancer hijacking events. This pivotal event has allowed the emergence of a new GRN in extant amphioxus, presumably through a stepwise process. In addition, the co-expression of Gdf1/3-like and Lefty achieved by a shared regulatory region may have provided robustness during body axis formation, which provides a selection-based hypothesis for the phenomena called developmental system drift.

Mitochondrial genome of Hancornia speciosa gomes: intergenic regions containing retrotransposons and predicted genes.

BACKGROUND: Plant mitochondrial genomes are characterized by high homologous recombination, extensive intergenic spacers, conservation in DNA sequences, and gene content. The Hancornia genus belongs to the Apocynaceae family, with H. speciosa Gomes being the sole species in the genus. It is an siganificant commercial fruit crop; however, only a number of studies have been conducted. In this study, we present the mitochondrial genome of H. speciosa and compare it with other mitochondrial genomes within the Apocynaceae family. METHODS AND RESULTS: A total of 2.8 Gb of Illumina paired-end reads were used to obtain the mitogenome, resulting in 22 contigs that were merged using 6.1 Gb of Illumina mate-pair reads to obtain a circular chromosome. The mitochondrial genome of H. speciosa is circular, containing 63 predicted functional genes, spanning a length of 741,811 bp, with a CG content of 44%. Within the mitogenome, 50 chloroplast DNA sequences, equivalent to 1.72% of the genome, were detected. However, intergenic spaces accounted for 703,139 bp (94.79% of the genome), and 287 genes were predicted, totaling 173,721 bp. CONCLUSION: This suggests the incorporation of nuclear DNA into the mitogenome of H. speciosa and self duplication. Comparative analysis among the mitogenomes in the Apocynaceae family revealed a diversity in the structure mediated by recombination, with similar gene content and large intergenic spaces.

Genomic and cytogenetic analyses reveal satellite repeat signature in allotetraploid okra (Abelmoschus esculentus).

BACKGROUND: Satellite repeats are one of the most rapidly evolving components in eukaryotic genomes and play vital roles in genome regulation, genome evolution, and speciation. As a consequence, the composition, abundance and chromosome distribution of satellite repeats often exhibit variability across various species, genome, and even individual chromosomes. However, we know little about the satellite repeat evolution in allopolyploid genomes. RESULTS: In this study, we investigated the satellite repeat signature in five okra (Abelmoschus esculentus) accessions using genomic and cytogenetic methods. In each of the five accessions, we identified eight satellite repeats, which exhibited a significant level of intraspecific conservation. Through fluorescence in situ hybridization (FISH) experiments, we observed that the satellite repeats generated multiple signals and exhibited variations in copy number across chromosomes. Intriguingly, we found that five satellite repeats were interspersed with centromeric retrotransposons, signifying their involvement in centromeric satellite repeat identity. We confirmed subgenome-biased amplification patterns of these satellite repeats through existing genome assemblies or dual-color FISH, indicating their distinct dynamic evolution in the allotetraploid okra subgenome. Moreover, we observed the presence of multiple chromosomes harboring the 35 S rDNA loci, alongside another chromosomal pair carrying the 5 S rDNA loci in okra using FISH assay. Remarkably, the intensity of 35 S rDNA hybridization signals varied among chromosomes, with the signals predominantly localized within regions of relatively weak DAPI staining, associated with GC-rich heterochromatin regions. Finally, we observed a similar localization pattern between 35 S rDNA and three satellite repeats with high GC content and confirmed their origin in the intergenic spacer region of the 35 S rDNA. CONCLUSIONS: Our findings uncover a unique satellite repeat signature in the allotetraploid okra, contributing to our understanding of the composition, abundance, and chromosomal distribution of satellite repeats in allopolyploid genomes, further enriching our understanding of their evolutionary dynamics in complex allopolyploid genomes.

Whole-genome resequencing facilitates the development of a 50K single nucleotide polymorphism genotyping array for Scots pine (Pinus sylvestris L.) and its transferability to other pine species.

Scots pine (Pinus sylvestris L.) is one of the most widespread and economically important conifer species in the world. Applications like genomic selection and association studies, which could help accelerate breeding cycles, are challenging in Scots pine because of its large and repetitive genome. For this reason, genotyping tools for conifer species, and in particular for Scots pine, are commonly based on transcribed regions of the genome. In this article, we present the Axiom Psyl50K array, the first single nucleotide polymorphism (SNP) genotyping array for Scots pine based on whole-genome resequencing, that represents both genic and intergenic regions. This array was designed following a two-step procedure: first, 192 trees were sequenced, and a 430K SNP screening array was constructed. Then, 480 samples, including haploid megagametophytes, full-sib family trios, breeding population, and range-wide individuals from across Eurasia were genotyped with the screening array. The best 50K SNPs were selected based on quality, replicability, distribution across the draft genome assembly, balance between genic and intergenic regions, and genotype-environment and genotype-phenotype associations. Of the final 49 877 probes tiled in the array, 20 372 (40.84%) occur inside gene models, while the rest lie in intergenic regions. We also show that the Psyl50K array can yield enough high-confidence SNPs for genetic studies in pine species from North America and Eurasia. This new genotyping tool will be a valuable resource for high-throughput fundamental and applied research of Scots pine and other pine species.

Histone H4 asymmetrically dimethylated at arginine 3 (H4R3me2a), a mark of super-enhancers.

Histone H4 asymmetrically dimethylated at arginine 3 (H4R3me2a) is an active histone mark catalyzed by protein arginine methyltransferase 1 (PRMT1), a major arginine methyltransferase in vertebrates catalyzing asymmetric dimethylation of arginine. H4R3me2a stimulates the activity of lysine acetyltransferases such as CBP/p300, which catalyze the acetylation of H3K27, a mark of active enhancers, super-enhancers, and promoters. There are a few studies on the genomic location of H4R3me2a. In chicken polychromatic erythrocytes, H4R3me2a is found in introns and intergenic regions and binds to the globin locus control region (a super-enhancer) and globin regulatory regions. In this report, we analyzed chromatin immunoprecipitation sequencing data for the genomic location of H4R3me2a in the breast cancer cell line MCF7. As in avian cells, MCF7 H4R3me2a is present in intronic and intergenic regions. Nucleosomes with H4R3me2a and H3K27ac next to nucleosome-free regions are found at super-enhancers, enhancers, and promoter regions of expressed genes. Genes with critical roles in breast cancer cells have broad domains of nucleosomes with H4R3me2a, H3K27ac, and H3K4me3. Our results are consistent with PRMT1-mediated H4R3me2a playing a key role in the chromatin organization of regulatory regions of vertebrate genomes.

Whole genome sequencing reveals candidate genes involving in PAS resistance in M. Tuberculosis isolated from patients in Thailand.

Para-amino salicylic acid (PAS) was first reported by Lehmann in 1946 and used for tuberculosis treatment. However, due to its adverse effects, it is now used only as a second line anti-tuberculosis drug for treatment of multidrug resistant or extensively drug resistant M. tuberculosis. The structure of PAS is similar to para-amino benzoic acid (pABA), an intermediate metabolite in the folate synthesis pathway. The study has identified mutations in genes in folate pathway and their intergenic regions for their possibilities in responsible for PAS resistance. Genomic DNA from 120 PAS-resistant and 49 PAS-sensitive M. tuberculosis isolated from tuberculosis patients in Thailand were studied by whole genome sequencing. Twelve genes in the folate synthesis pathway were investigated for variants associated with PAS resistance. Fifty-one SNVs were found in nine genes and their intergenic regions (pabC, pabB, folC, ribD, thyX, dfrA, thyA, folK, folP). Functional correlation test confirmed mutations in RibD, ThyX, and ThyA are responsible for PAS resistance. Detection of mutation in thyA, folC, intergenic regions of thyX, ribD, and double deletion of thyA dfrA are proposed for determination of PAS resistant M. tuberculosis.

Genetic Modification of a Hox Locus Drives Mimetic Color Pattern Variation in a Highly Polymorphic Bumble Bee.

Müllerian mimicry provides natural replicates ideal for exploring mechanisms underlying adaptive phenotypic divergence and convergence, yet the genetic mechanisms underlying mimetic variation remain largely unknown. The current study investigates the genetic basis of mimetic color pattern variation in a highly polymorphic bumble bee, Bombus breviceps (Hymenoptera, Apidae). In South Asia, this species and multiple comimetic species converge onto local Müllerian mimicry patterns by shifting the abdominal setal color from orange to black. Genetic crossing between the orange and black phenotypes suggested the color dimorphism being controlled by a single Mendelian locus, with the orange allele being dominant over black. Genome-wide association suggests that a locus at the intergenic region between 2 abdominal fate-determining Hox genes, abd-A and Abd-B, is associated with the color change. This locus is therefore in the same intergenic region but not the same exact locus as found to drive red black midabdominal variation in a distantly related bumble bee species, Bombus melanopygus. Gene expression analysis and RNA interferences suggest that differential expression of an intergenic long noncoding RNA between abd-A and Abd-B at the onset setal color differentiation may drive the orange black color variation by causing a homeotic shift late in development. Analysis of this same color locus in comimetic species reveals no sequence association with the same color shift, suggesting that mimetic convergence is achieved through distinct genetic routes. Our study establishes Hox regions as genomic hotspots for color pattern evolution in bumble bees and demonstrates how pleiotropic developmental loci can drive adaptive radiations in nature.