A novel PCR-based genotyping method for Proteus mirabilis - Intergenic region polymorphism analysis.

Proteus mirabilis is a predominant species in cases of food poisoning associated with meat products and is also an opportunistic pathogen causing numerous infections in humans. This study aimed to differentiate P. mirabilis isolates using intergenic region polymorphism analysis (IRPA). The IRPA typing scheme was developed to amplify polymorphic fragments in intergenic regions (IGRs). The presence, absence, or size change of amplified products were identified and utilized as genetic markers for rapid differentiation of strains. A total of 75 P. mirabilis isolates were isolated from 63 fresh poultry and pork samples were subtyped using the IRPA and ERIC-PCR methods, and their antibiotic resistance profiles were tested. The majority of P. mirabilis isolates showed resistance to tetracycline (85.3%), doxycycline (93.3%), chloramphenicol (82.7%), streptomycin (92.0%), spectinomycin (80.0%), trimethoprim (97.3%); trimethoprim-sulfalleth (82.7%), and erythromycin (100.0%). In contrast, resistance rates to ceftriaxon, cefoxitin, cefepime, and cefotaxim were lower at only 17.3%, 5.3%, 6.7%, and 13.3%, respectively, among P. mirabilis isolates. Eleven loci were selected for analysis of the genetic diversity of 75 P. mirabilis isolates. A combination of 4 loci was determined as the optimal combination. The results compared to those obtained using ERIC-PCR for the same isolates. The Simpson's index of diversity was 0.999 for IRPA and 0.923 for ERIC-PCR, indicating that IRPA has a higher discriminatory power than ERIC-PCR. The concordance between IRPA and ERIC-PCR methods was low, primarily because IRPA classified isolates from the same ERIC cluster into separate clusters due to its high resolution. The IRPA method presented in this study offers a rapid, simple, reproducible, and economical approach for genotyping P. mirabilis.

Meta-analysis of genome-wide association studies for cancer therapy-related cardiovascular dysfunction and functional mapping highlight an intergenic region close to TP63.

Cancer therapy-related cardiac dysfunction (CTRCD), which commonly includes left ventricular dysfunction and heart failure, is the main adverse effect of anticancer therapy. In recent years several candidate genes studies and genome-wide association studies have identified common genetic variants associated with CTRCD, but evidence remains limited and few genetic variants are robust. A genome-wide meta-analysis of CTRCD was performed with 852 oncology patients receiving cancer therapy. DNA samples were genotyped and imputed to perform a GWAS meta-analysis for case-control (N = 852 (380 cases and 472 controls) and extreme phenotypes (N = 618 (78 cases and 472 controls) looking for genetic variants that predispose to CTRCD. The results were validated in a replicate cohort of 1,191 oncology patients (245 cases and 946 controls). Functional mapping of the replicated loci was then performed. The meta-analysis showed 9 and 17 loci suggestively associated (P-value < 1 × 10-5) with CTRCD in case-control and extreme phenotypes analyses, respectively. The 3q28 locus (rs rs7652759, P = 5.64 × 10-6) in the case-control analysis was the strongest signal, with up to 64 SNPs above the suggestive significance threshold. The rs7652759, an intergenic variant between TPRG1 and TP63 genes, was the only variant validated in the replication cohort (P-value = 0.01). Functional mapping of this significant locus revealed up to 5 new genes potentially involved in the CTRCD. We identified the intergenic region near TP63 as a novel CTRCD susceptibility locus. In the future, the genotyping of these markers could be considered in new CTRCD risk scores to improve preventive strategies in cardio-oncology.

Engineering spacer specificity of the Cre/loxP system.

Translational research on the Cre/loxP recombination system focuses on enhancing its specificity by modifying Cre/DNA interactions. Despite extensive efforts, the exact mechanisms governing Cre discrimination between substrates remains elusive. Cre recognizes 13 bp inverted repeats, initiating recombination in the 8 bp spacer region. While literature suggests that efficient recombination proceeds between lox sites with non-loxP spacer sequences when both lox sites have matching spacers, experimental validation for this assumption is lacking. To fill this gap, we investigated target site variations of identical pairs of the loxP 8 bp spacer region, screening 6000 unique loxP-like sequences. Approximately 84% of these sites exhibited efficient recombination, affirming the plasticity of spacer sequences for catalysis. However, certain spacers negatively impacted recombination, emphasizing sequence dependence. Directed evolution of Cre on inefficiently recombined spacers not only yielded recombinases with enhanced activity but also mutants with reprogrammed selective activity. Mutations altering spacer specificity were identified, and molecular modelling and dynamics simulations were used to investigate the possible mechanisms behind the specificity switch. Our findings highlight the potential to fine-tune site-specific recombinases for spacer sequence specificity, offering a novel concept to enhance the applied properties of designer-recombinases for genome engineering applications.

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.

Development of DNA Insertion-specific Markers Based on the Intergenic Region of Oplegnathus punctatus Cdkn1/srsf3 for Sex Identification.

Spotted knifejaw (Oplegnathus punctatus) is a marine economic fish with high food and ecological value, and its growth process has obvious male and female sexual dimorphism, with males growing significantly faster than females. However, the current sex identification technology is not yet mature, which will limit the growth rate of O. punctatus aquaculture and the efficiency of separate sex breeding, so the development of efficient sex molecular markers is imperative. This study identified a 926 bp DNA insertion fragment in the cdkn1/srsf3 intergenic region of O. punctatus males through whole-genome scanning, comparative genomics, and structural variant analysis. A pair of primers was designed based on the insertion information of the Y chromosome intergenic region in male individuals. Agarose gel electrophoresis revealed the amplification of two DNA fragments, 1118 bp and 192 bp, in male O. punctatus individuals. The 926 bp fragment was identified as the insertion in the intergenic region of cdkn1/srsf3 in males, while only a single 192 bp DNA fragment was amplified in females. The biological sex of the individuals identified in this manner was consistent with their known phenotypic sex. In this study, we developed a method to detect DNA insertion variants in the intergenic region of O. punctatus. Additionally, we introduced a new DNA marker for the rapid identification of the sex of O. punctatus, which enhances detection efficiency. The text has important reference significance and application value in sex identification, all-male breeding, and lineage selection. It provides new insights into the regulation of variation in the intergenic region of cdkn1/srsf3 genes and the study of RNA shearing.

Unraveling noncoding DNA variants and epimutations: a paradigm shift in hereditary cancer research.

Exhaustive efforts have been dedicated to uncovering genomic aberrations linked to cancer susceptibility. Noncoding sequence variants and epigenetic alterations significantly influence gene regulation and could contribute to cancer development. However, exploring noncoding regions in hereditary cancer susceptibility demands cutting-edge methodologies for functionally characterizing genomic discoveries. Additionally, comprehending the impact on cancer development of variants in noncoding DNA and the epigenome necessitates integrating diverse data through bioinformatic analyses. As novel technologies and analytical methods continue to advance, this realm of research is rapidly gaining traction. Within this mini-review, we delve into future research domains concerning aberrations in noncoding DNA regions, such as pseudoexons, promoter variants and cis-epimutations.

[Box: see text].

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.