Absence of Figla-like Gene Is Concordant with Femaleness in Cichlids Harboring the LG1 Sex-Determination System.

Oreochromis niloticus has been used as a reference genome for studies of tilapia sex determination (SD) revealing segregating genetic loci on linkage groups (LGs) 1, 3, and 23. The master key regulator genes (MKR) underlying the SD regions on LGs 3 and 23 have been already found. To identify the MKR in fish that segregate for the LG1 XX/XY SD-system, we applied short variant discovery within the sequence reads of the genomic libraries of the Amherst hybrid stock, Coptodon zillii and Sarotherodon galilaeus, which were aligned to a 3-Mbp-region of the O. aureus genome. We obtained 66,372 variants of which six were concordant with the XX/XY model of SD and were conserved across these species, disclosing the male specific figla-like gene. We further validated this observation in O. mossambicus and in the Chitralada hybrid stock. Genome alignment of the 1252-bp transcript showed that the figla-like gene's size was 2664 bp, and that its three exons were capable of encoding 99 amino acids including a 45-amino-acid basic helix-loop-helix domain that is typical of the ovary development regulator-factor-in-the-germline-alpha (FIGLA). In Amherst gonads, the figla-like gene was exclusively expressed in testes. Thus, the figla-like genomic presence determines male fate by interrupting the female developmental program. This indicates that the figla-like gene is the long-sought SD MKR on LG1.

Gene Variant of Barrier to Autointegration Factor 2 (Banf2w) Is Concordant with Female Determination in Cichlids.

Oreochromis fishes exhibit variability of sex-determination (SD) genes whose characterization contributes to understanding of the sex differentiation network, and to effective tilapia farming, which requires all-male culture. However, O. niloticus (On) amh is the only master-key regulator (MKR) of SD that has been mapped (XY/XX SD-system on LG23). In O. aureus (Oa), LG3 controls a WZ/ZZ SD-system that has recently been delimited to 9.2 Mbp, with an embedded interval rich with female-specific variation, harboring two paics genes and banf2. Developing genetic markers within this interval and using a hybrid Oa stock that demonstrates no recombination repression in LG3, we mapped the critical SD region to 235 Kbp on the orthologous On physical map (p < 1.5 × 10-26). DNA-seq assembly and peak-proportion analysis of variation based on Sanger chromatograms allowed the characterization of copy-number variation (CNV) of banf2. Oa males had three exons capable of encoding 90-amino-acid polypeptides, yet in Oa females, we found an extra copy with an 89-amino-acid polypeptide and three non-conservative amino acid substitutions, designated as banf2w. CNV analysis suggested the existence of two to five copies of banf2 in diploidic Cichlidae. Disrupting the Hardy-Weinberg equilibrium (p < 4.2 × 10-3), banf2w was concordant with female determination in Oa and in three cichlids with LG3 WZ/ZZ SD-systems (O. tanganicae, O. hornorum and Pelmatolapia mariae). Furthermore, exclusive RNA-seq expression in Oa females strengthened the candidacy of banf2w as the long-sought LG3 SD MKR. As banf genes mediate nuclear assembly, chromatin organization, gene expression and gonad development, banf2w may play a fundamental role inducing female nucleus formation that is essential for WZ/ZZ SD.

Seasonal Dynamics Are the Major Driver of Microbial Diversity and Composition in Intensive Freshwater Aquaculture.

Aquaculture facilities such as fishponds are one of the most anthropogenically impacted freshwater ecosystems. The high fish biomass reared in aquaculture is associated with an intensive input into the water of fish-feed and fish excrements. This nutrients load may affect the microbial community in the water, which in turn can impact the fish health. To determine to what extent aquaculture practices and natural seasonal cycles affect the microbial populations, we characterized the microbiome of an inter-connected aquaculture system at monthly resolution, over 3 years. The system comprised two fishponds, where fish are grown, and an operational water reservoir in which fish are not actively stocked. Clear natural seasonal cycles of temperature and inorganic nutrients concentration, as well as recurring cyanobacterial blooms during summer, were observed in both the fishponds and the reservoir. The structure of the aquatic bacterial communities in the system, characterized using 16S rRNA sequencing, was explained primarily by the natural seasonality, whereas aquaculture-related parameters had only a minor explanatory power. However, the cyanobacterial blooms were characterized by different cyanobacterial clades dominating at each fishpond, possibly in response to distinct nitrogen and phosphate ratios. In turn, nutrient ratios may have been affected by the magnitude of fish feed input. Taken together, our results show that, even in strongly anthropogenically impacted aquatic ecosystems, the structure of bacterial communities is mainly driven by the natural seasonality, with more subtle effects of aquaculture-related factors.

Temporal Resistome and Microbial Community Dynamics in an Intensive Aquaculture Facility with Prophylactic Antimicrobial Treatment.

Excessive use of antimicrobials in aquaculture is concerning, given possible environmental ramifications and the potential contribution to the spread of antimicrobial resistance (AR). In this study, we explored seasonal abundance of antimicrobial resistance genes and bacterial community composition in the water column of an intensive aquaculture pond stocked with Silver Carp (Hypophthalmichthys molitrix) prophylactically treated with sulfamethoprim (25% sulfadiazine; 5% trimethoprim), relative to an adjacent unstocked reservoir. Bacterial community composition was monitored using high-throughput sequencing of 16S rRNA gene amplicons in eight sampling profiles to determine seasonal dynamics, representing principal stages in the fish fattening cycle. In tandem, qPCR was applied to assess relative abundance of selected antimicrobial resistance genes (sul1, sul2, dfrA1, tetA and blaTEM) and class-1 integrons (int1). Concomitantly, resistomes were extrapolated from shotgun metagenomes in representative profiles. Analyses revealed increased relative abundance of sulfonamide and tetracycline resistance genes in fishpond-03, relative to pre-stocking and reservoir levels, whereas no significant differences were observed for genes encoding resistance to antimicrobials that were not used in the fishpond-03. Seasons strongly dictated bacterial community composition, with high abundance of cyanobacteria in summer and increased relative abundance of Flavobacterium in the winter. Our results indicate that prophylactic use of sulfonamides in intensive aquaculture ponds facilitates resistance suggesting that prophylactic use of these antimicrobials in aquaculture should be restricted.

Evidence of Increased Antibiotic Resistance in Phylogenetically-Diverse Aeromonas Isolates from Semi-Intensive Fish Ponds Treated with Antibiotics.

The genus Aeromonas is ubiquitous in aquatic environments encompassing a broad range of fish and human pathogens. Aeromonas strains are known for their enhanced capacity to acquire and exchange antibiotic resistance genes and therefore, are frequently targeted as indicator bacteria for monitoring antimicrobial resistance in aquatic environments. This study evaluated temporal trends in Aeromonas diversity and antibiotic resistance in two adjacent semi-intensive aquaculture facilities to ascertain the effects of antibiotic treatment on antimicrobial resistance. In the first facility, sulfadiazine-trimethoprim was added prophylactically to fingerling stocks and water column-associated Aeromonas were monitored periodically over an 11-month fish fattening cycle to assess temporal dynamics in taxonomy and antibiotic resistance. In the second facility, Aeromonas were isolated from fish skin ulcers sampled over a 3-year period and from pond water samples to assess associations between pathogenic strains to those in the water column. A total of 1200 Aeromonas isolates were initially screened for sulfadiazine resistance and further screened against five additional antimicrobials. In both facilities, strong correlations were observed between sulfadiazine resistance and trimethoprim and tetracycline resistances, whereas correlations between sulfadiazine resistance and ceftriaxone, gentamicin, and chloramphenicol resistances were low. Multidrug resistant strains as well as sul1, tetA, and intI1 gene-harboring strains were significantly higher in profiles sampled during the fish cycle than those isolated prior to stocking and these genes were extremely abundant in the pathogenic strains. Five phylogenetically distinct Aeromonas clusters were identified using partial rpoD gene sequence analysis. Interestingly, prior to fingerling stocking the diversity of water column strains was high, and representatives from all five clusters were identified, including an A. salmonicida cluster that harbored all characterized fish skin ulcer samples. Subsequent to stocking, diversity was much lower and most water column isolates in both facilities segregated into an A. veronii-associated cluster. This study demonstrated a strong correlation between aquaculture, Aeromonas diversity and antibiotic resistance. It provides strong evidence for linkage between prophylactic and systemic use of antibiotics in aquaculture and the propagation of antibiotic resistance.