Mitochondrial DNA and Y chromosome reveal the genetic structure of the native Polish Konik horse population.

Polish Konik remains one of the most important horse breeds in Poland. The primitive, native horses with a stocky body and mouse-like coat color are protected by a conservation program, while their Polish population consists of about 3,480 individuals, representing 16 dam and six sire lines. To define the population's genetic structure, mitochondrial DNA and Y chromosome sequence variables were identified. The mtDNA whole hypervariable region analysis was carried out using the Sanger sequencing method on 233 Polish Koniks belonging to all dam lines, while the Y chromosome analysis was performed with the competitive allele-specific PCR genotyping method on 36 horses belonging to all sire lines. The analysis of the mtDNA hypervariable region detected 47 SNPs, which assigned all tested horses to 43 haplotypes. Most dam lines presented more than one haplotype; however, five dam lines were represented by only one haplotype. The haplotypes were classified into six (A, B, E, J, G, R) recognized mtDNA haplogroups, with most horses belonging to haplogroup A, common among Asian horse populations. Y chromosome analysis allocated Polish Koniks in the Crown group, condensing all modern horse breeds, and divided them into three haplotypes clustering with coldblood breeds (28 horses), warmblood breeds (two horses), and Duelmener Pony (six horses). The clustering of all Wicek sire line stallions with Duelmener horses may suggest a historical relationship between the breeds. Additionally, both mtDNA and Y chromosome sequence variability results indicate crossbreeding before the studbooks closure or irregularities in the pedigrees occurred before the DNA testing introduction.

Y-Chromosome Haplotype Report among Eight Italian Horse Breeds.

Horse domestication and breed selection processes have profoundly influenced the development and transformation of human society and civilization over time. Therefore, their origin and history have always attracted much attention. In Italy, several local breeds have won prestigious awards thanks to their unique traits and socio-cultural peculiarities. Here, for the first time, we report the genetic variation of three loci of the male-specific region of the Y chromosome (MSY) of four local breeds and another one (Lipizzan, UNESCO) well-represented in the Italian Peninsula. The analysis also includes data from three Sardinian breeds and another forty-eight Eurasian and Mediterranean horse breeds retrieved from GenBank for comparison. Three haplotypes (HT1, HT2, and HT3) were found in Italian stallions, with different spatial distributions between breeds. HT1 (the ancestral haplotype) was frequent, especially in Bardigiano and Monterufolino, HT2 (Neapolitan/Oriental wave) was found in almost all local breeds, and HT3 (Thoroughbred wave) was detected in Maremmano and two Sardinian breeds (Sardinian Anglo-Arab and Sarcidano). This differential distribution is due to three paternal introgressions of imported stallions from foreign countries to improve local herds; however, further genetic analyses are essential to reconstruct the genetic history of native horse breeds, evaluate the impact of selection events, and enable conservation strategies.

GISH painting of the Y chromosomes suggests advanced phases of sex chromosome evolution in three dioecious Cannabaceae species (Humulus lupulus, H. japonicus, and Cannabis sativa).

In plants, dioecy is relatively rare, and it involves sex chromosome systems that often developed independently over time. These characteristics make dioecious plants an attractive model to study sex chromosome evolution. To clarify the patterns of plant sex chromosome evolution, studies should be performed on a wide range of dioecious species. It is interesting to study the sex chromosomes in related species that evolved during a long period of independent sex chromosome evolution. The Cannabaceae family includes three dioecious species with heteromorphic sex chromosomes. Cannabis sativa and Humulus lupulus use the XX/XY chromosome system, whereas Humulus japonicus contains multiple sex chromosomes (XX/XY1Y2). To better understand sex chromosome evolution and the level of genomic divergence of these three related species, we undertook self-GISH and comparative GISH analyses. The self-GISH allowed visualization of the Y chromosomes of C. sativa, H. lupulus, and H. japonicus. The self-GISH signal was distributed along the entire Y chromosome, excluding the pseudo-autosomal region (PAR). Our results indicate that the male-specific region of the Y chromosome (MSY) spans the overwhelming majority of the Y chromosomes of all three species studied. The self-GISH results reveal the accumulation of repetitive DNA sequences in the Y chromosomes of all three species studied. This sequences presented in autosomes and/or chromosome X at a lower copy number than in Y. In comparative GISH experiments where the probe DNA of one species was applied to another species, a weak signal was exclusively detected on 45S rDNA sites, indicating a high level of genomic differentiation of the species used in this study. We demonstrate small PAR size and opposing large MSY and its positions on Y chromosomes. We also found that these genomes are highly differentiated. Furthermore, the data obtained in this study indicate a long period of independent and advanced sex chromosome evolution. Our study provides a valuable basis for future genomic studies of sex and suggests that the Cannabaceae family offers a promising model to study sex chromosome evolution.

Unlocking Horse Y Chromosome Diversity.

The present equine genetic variation mirrors the deep influence of intensive breeding programs during the last 200 years. Here, we provide a comprehensive current state of knowledge on the trends and prospects on the variation in the equine male-specific region of the Y chromosome (MSY), which was assembled for the first time in 2018. In comparison with the other 12 mammalian species, horses are now the most represented, with 56 documented MSY genes. However, in contrast to the high variability in mitochondrial DNA observed in many horse breeds from different geographic areas, modern horse populations demonstrate extremely low genetic Y-chromosome diversity. The selective pressures employed by breeders using pedigree data (which are not always error-free) as a predictive tool represent the main cause of this lack of variation in the Y-chromosome. Nevertheless, the detailed phylogenies obtained by recent fine-scaled Y-chromosomal genotyping in many horse breeds worldwide have contributed to addressing the genealogical, forensic, and population questions leading to the reappraisal of the Y-chromosome as a powerful genetic marker to avoid the loss of biodiversity as a result of selective breeding practices, and to better understand the historical development of horse breeds.

Model of two competing populations in two habitats with migration: Application to optimal marine protected area size.

The standard model of a single population fragmented into two patches connected by migration, was first introduced in the 1970s by Freedman and Waltman, since generating long-term research interest, though its full analysis for arbitrary values of migration rate has only been completed relatively recently. Here, we present a model of two competing species in a two-patch habitat with migrations between patches. We derive equilibrium solutions of this model for three cases of migration rate resulting in isolated, well-mixed and semi-isolated habitats. We evaluate the full range of effects of habitat, life-history and migration parameters on population sizes. Finally, we add harvesting mortality and define conditions under which introduction of a no-harvesting (protected) area may lead to increased maximum sustainable yield. The results have applications in mixed fishery management and the design of wildlife protection zones, including marine protected areas (MPAs).

Sex-specific gene expression in the blood of four primates.

Blood is an important non-reproductive tissue, but little is known about the sex-specific gene expressions in the blood. Therefore, we investigated sex-specific gene expression differences in the blood tissues of four primates, rhesus macaques (Macaca mulatta), Tibetan macaques (M. thibetana), yellow baboons (Papio cynocephalus), and humans. We identified seven sex-specific differentially expressed genes (SDEGs) in each non-human primate and 31 SDEGs in humans. The four primates had only one common SDEG, MAP7D2. In humans, immune-related SDEGs were identified as up-regulated, but also down-regulated in females. We also found that most of the X-Y gene pairs had similar expression levels between species, except pair EIF1AY/EIF1AX. The expression level of X-Y gene pairs of rhesus and Tibetan macaques showed no significant differential expression levels, while humans had six significant XY-biased and three XX-biased X-Y gene pairs. Our observed sex differences in blood should increase understanding of sex differences in primate blood tissue.

MFN2 interacts with nuage-associated proteins and is essential for male germ cell development by controlling mRNA fate during spermatogenesis.

Mitochondria play a crucial role in spermatogenesis and are regulated by several mitochondrial fusion proteins. However, their functional importance associated with their structure formation and mRNA fate regulation during spermatogenesis remains unclear. Here, we show that mitofusin 2 (MFN2), a mitochondrial fusion protein, interacts with nuage-associated proteins (including MIWI, DDX4, TDRKH and GASZ) in mice. Conditional mutation of Mfn2 in postnatal germ cells results in male sterility due to germ cell developmental defects. Moreover, MFN2 interacts with MFN1, another mitochondrial fusion protein with a high-sequence similarity to MFN2, in testes to facilitate spermatogenesis. Simultaneous mutation of Mfn1 and Mfn2 in testes causes very severe infertile phenotypes. Importantly, we show that MFN2 is enriched in polysome fractions of testes and interacts with MSY2, a germ cell-specific DNA/RNA-binding protein, to control gamete-specific mRNA (such as Spata19) translational activity during spermatogenesis. Collectively, our findings demonstrate that MFN2 interacts with nuage-associated proteins and MSY2 to regulate male germ cell development by controlling several gamete-specific mRNA fates.

Association of MSY haplotype background with nonobstructive azoospermia is AZF-dependent: A case-control study.

Identifying causal genes of spermatogenic failure on the male-specific region of Y chromosome (MSY) has been a challenging process. Due to the nonrecombining nature of MSY, haplotype-based approaches have recently been shown to be promising in identifying associated MSY haplogroups. We conducted an MSY analysis of nonobstructive azoospermia (NOA) patients in a case-control setting (N = 278 and 105 respectively) to identify modal haplogroups strongly associated with NOA. Patients with AZF deletions (AZF+) and no AZF deletions (AZF-) were compared with the control group. Given the larger sample set of AZF- NOA patients, we further investigated the association based on histopathological severity, namely Sertoli cell-only syndrome and maturation arrest subtypes. We observed no significant enrichment of MSY haplogroups in AZF- azoospermic patients (or its subtypes). However, we observed a strongly significant association between haplogroup J2a* and AZF+ patients (FDR-corrected p = .0056; OR = 7.02, 95%CI 1.89 to 39.20), a haplogroup which also showed significant enrichment for AZFa/b deletions (p = 4x10-4 ). We conclude that unlike AZF+ patients, AZF- NOA are less likely to have an MSY causative factor with large effect size, thus indicating that the aetiology of AZF- NOA, and to some extent AZFc NOA, is more likely to be based on non-MSY factors.

Population Genetic Analysis of the Estonian Native Horse Suggests Diverse and Distinct Genetics, Ancient Origin and Contribution from Unique Patrilines.

The Estonian Native Horse (ENH) is a medium-size pony found mainly in the western islands of Estonia and is well-adapted to the harsh northern climate and poor pastures. The ancestry of the ENH is debated, including alleged claims about direct descendance from the extinct Tarpan. Here we conducted a detailed analysis of the genetic makeup and relationships of the ENH based on the genotypes of 15 autosomal short tandem repeats (STRs), 18 Y chromosomal single nucleotide polymorphisms (SNPs), mitochondrial D-loop sequence and lateral gait allele in DMRT3. The study encompassed 2890 horses of 61 breeds, including 33 ENHs. We show that the expected and observed genetic diversities of the ENH are among the highest within 52 global breeds, and the highest among 8 related Northern European ponies. The genetically closest breeds to the ENH are the Finn Horse, and the geographically more distant primitive Hucul and Konik. ENH matrilines are diverse and relate to draught and Pontic-Caspian breeds. ENH patrilines relate to draught breeds, and to a unique haplogroup not described before. None of the 33 ENHs carried the "gait" mutation, but the mutation was found in 2 Huculs. The study demonstrates that the ENH is a genetically distinct and diverse breed of ancient origin with no notable pressure of selective breeding.

Cytogenetic and genomic organization analyses of chloroplast DNA invasions in the nuclear genome of Asparagus officinalis L. provides signatures of evolutionary complexity and informativity in sex chromosome evolution.

BACKGROUND: The transfer of chloroplast DNA into nuclear genome is a common process in plants. These transfers form nuclear integrants of plastid DNAs (NUPTs), which are thought to be driving forces in genome evolution, including sex chromosome evolution. In this study, NUPTs in the genome of a dioecious plant Asparagus officinalis L. were systematically analyzed, in order to investigate the characteristics of NUPTs in the nuclear genome and the relationship between NUPTs and sex chromosome evolution in this species. RESULTS: A total of 3155 NUPT insertions were detected, and they represented approximated 0.06% of the nuclear genome. About 45% of the NUPTs were organized in clusters. These clusters were derived from various evolutionary events. The Y chromosome contained the highest number and largest proportion of NUPTs, suggesting more accumulation of NUPTs on sex chromosomes. NUPTs were distributed widely in all of the chromosomes, and some regions preferred these insertions. The highest density of NUPTs was found in a 47 kb region in the Y chromosome; more than 75% of this region was occupied by NUPTs. Further cytogenetic and sequence alignment analysis revealed that this region was likely the centromeric region of the sex chromosomes. On the other hand, the male-specific region of the Y chromosome (MSY) and the adjacent regions did not have NUPT insertions. CONCLUSIONS: These results indicated that NUPTs were involved in shaping the genome of A. officinalis through complicated process. NUPTs may play important roles in the centromere shaping of the sex chromosomes of A. officinalis, but were not implicated in MSY formation.

The Red Fox Y-Chromosome in Comparative Context.

While the number of mammalian genome assemblies has proliferated, Y-chromosome assemblies have lagged behind. This discrepancy is caused by biological features of the Y-chromosome, such as its high repeat content, that present challenges to assembly with short-read, next-generation sequencing technologies. Partial Y-chromosome assemblies have been developed for the cat (Feliscatus), dog (Canislupusfamiliaris), and grey wolf (Canislupuslupus), providing the opportunity to examine the red fox (Vulpesvulpes) Y-chromosome in the context of closely related species. Here we present a data-driven approach to identifying Y-chromosome sequence among the scaffolds that comprise the short-read assembled red fox genome. First, scaffolds containing genes found on the Y-chromosomes of cats, dogs, and wolves were identified. Next, analysis of the resequenced genomes of 15 male and 15 female foxes revealed scaffolds containing male-specific k-mers and patterns of inter-sex copy number variation consistent with the heterogametic chromosome. Analyzing variation across these two metrics revealed 171 scaffolds containing 3.37 Mbp of putative Y-chromosome sequence. The gene content of these scaffolds is consistent overall with that of the Y-chromosome in other carnivore species, though the red fox Y-chromosome carries more copies of BCORY2 and UBE1Y than has been reported in related species and fewer copies of SRY than in other canids. The assignment of these scaffolds to the Y-chromosome serves to further characterize the content of the red fox draft genome while providing resources for future analyses of canid Y-chromosome evolution.

Managing yield and resilience in a harvested tri-trophic food chain model.

In this article, we compare the two ecological services known as yield and resilience, for a tri-trophic food chain model consisting of a prey, an intermediate predator and a top predator. For this comparison process, we use both analytical and numerical techniques. It is shown that a variety of patterns are possible based on the intensity of efforts distributed among different trophic levels. Thus we may suggest that fishing down the food chain, as suggested by Pauly et al. (1998) is not bound to happen. Our analysis also shows that balancing the harvest between prey, intermediate predator and top predator could give more yield and stabilizing the ecosystem, than the selective harvesting of any one species. This balanced harvesting may not be a win-win situation for yield and resilience, but it could be a most favourable strategy to balance them. This research would help to correlate resilience with yield and determines the desirable selection of two policies, resilience maximizing yield or maximum sustainable yield to safeguard ecological communities.

The peopling of the last Green Sahara revealed by high-coverage resequencing of trans-Saharan patrilineages.

BACKGROUND: Little is known about the peopling of the Sahara during the Holocene climatic optimum, when the desert was replaced by a fertile environment. RESULTS: In order to investigate the role of the last Green Sahara in the peopling of Africa, we deep-sequence the whole non-repetitive portion of the Y chromosome in 104 males selected as representative of haplogroups which are currently found to the north and to the south of the Sahara. We identify 5,966 mutations, from which we extract 142 informative markers then genotyped in about 8,000 subjects from 145 African, Eurasian and African American populations. We find that the coalescence age of the trans-Saharan haplogroups dates back to the last Green Sahara, while most northern African or sub-Saharan clades expanded locally in the subsequent arid phase. CONCLUSIONS: Our findings suggest that the Green Sahara promoted human movements and demographic expansions, possibly linked to the adoption of pastoralism. Comparing our results with previously reported genome-wide data, we also find evidence for a sex-biased sub-Saharan contribution to northern Africans, suggesting that historical events such as the trans-Saharan slave trade mainly contributed to the mtDNA and autosomal gene pool, whereas the northern African paternal gene pool was mainly shaped by more ancient events.

Role for PADI6 in securing the mRNA-MSY2 complex to the oocyte cytoplasmic lattices.

The oocyte cytoplasmic lattices (CPLs) have long been predicted to function as a storage form for the maternal contribution of ribosomes to the early embryo. Our previous studies have demonstrated that ribosomal component S6 is stored in the oocyte CPLs and peptidylarginine deiminase 6 (PADI6) is critical for CPLs formation. Additionally, we found that depletion of PADI6 reduced de novo protein synthesis prior to the maternal-to-embryonic transition, therefore causing embryos to arrest at the 2-cell stage. Here, we present evidence further supporting the association of ribosomes with the CPLs by demonstrating that rRNAs are dramatically decreased in Padi6 KO oocytes. We also show that the abundance and localization of mRNAs is affected upon PADI6 depletion, suggesting that mRNAs are very possibly associated with CPLs. Consistent with this observation, the amount of the major RNA binding protein, MSY2, that is associated with the insoluble fraction of the oocytes after Triton X-100 extraction is also markedly decreased in the Padi6 KO oocytes. Furthermore, treatment of the oocytes with RNase A followed by Triton X-100 extraction severely impairs the localization of PADI6 and MSY2 in oocytes. These results indicate that mRNAs, possibly in a complex with MSY2 and PADI6, are bound in the CPLs and may play a role in securing the mRNA-MSY2 complex to the CPLs.

An explicit solution for calculating optimum spawning stock size from Ricker's stock recruitment model.

Stock-recruitment models have been used for decades in fisheries management as a means of formalizing the expected number of offspring that recruit to a fishery based on the number of parents. In particular, Ricker's stock recruitment model is widely used due to its flexibility and ease with which the parameters can be estimated. After model fitting, the spawning stock size that produces the maximum sustainable yield (S MSY) to a fishery, and the harvest corresponding to it (U MSY), are two of the most common biological reference points of interest to fisheries managers. However, to date there has been no explicit solution for either reference point because of the transcendental nature of the equation needed to solve for them. Therefore, numerical or statistical approximations have been used for more than 30 years. Here I provide explicit formulae for calculating both S MSY and U MSY in terms of the productivity and density-dependent parameters of Ricker's model.

The phenotypic impact of the male-specific region of chromosome-Y in inbred mating: the role of genetic variants and gene duplications in multiple inbred rat strains.

BACKGOUND: The male-specific region of chromosome-Y (MSY) contributes to phenotypes outside of testis development and has a high rate of evolution between mammalian species. With a lack of genomic crossover, MSY is one of the few genomic areas under similar variation and evolutionary selection in inbred and outbred animal populations, allowing for an assessment of evolutionary mechanisms to translate between the populations. METHODS: Using next-generation sequencing, MSY consomic strains, molecular characterization, and large-scale phenotyping, we present here regions of MSY that contribute to inbred strain phenotypes. RESULTS: We have shown that (1) MSY of rat has nine autosomal gene transposition events with strain-specific selection; (2) sequence variants in MSY occur with a 1.98-fold higher number of variants than other chromosomes in seven sequenced rat strains; (3) Sry, the most studied MSY gene, has undergone extensive gene duplications, driving ubiquitous expression not seen in human or mouse; (4) the expression profile of Sry in the rat is driven by the insertion of the Sry2 copy into an intron of the ubiquitously expressed Kdm5d gene in antisense orientation, but due to several loss of function mutations in the Sry2 protein, nuclear localization and transcriptional control are decreased; (5) expression of Sry copies other than Sry2 in the rat overlaps with the expression profile for human SRY; (6) gene duplications and sequence variants (P76T) of Sry can be selected for phenotypes such as high blood pressure and androgen receptor signaling within inbred mating; and most importantly, (7) per chromosome size, MSY contributes to higher strain-specific phenotypic variation relative to all other chromosomes, with 53 phenotypes showing both a male to female and consomic cross significance. CONCLUSION: The data presented supports a high probability of MSY genetic variation altering a broad range of inbred rat phenotypes.

The Biology and Evolution of Mammalian Y Chromosomes.

Mammals have the oldest sex chromosome system known: the mammalian X and Y chromosomes evolved from ordinary autosomes beginning at least 180 million years ago. Despite their shared ancestry, mammalian Y chromosomes display enormous variation among species in size, gene content, and structural complexity. Several unique features of the Y chromosome--its lack of a homologous partner for crossing over, its functional specialization for spermatogenesis, and its high degree of sequence amplification--contribute to this extreme variation. However, amid this evolutionary turmoil many commonalities have been revealed that have contributed to our understanding of the selective pressures driving the evolution and biology of the Y chromosome. Two biological themes have defined Y-chromosome research over the past six decades: testis determination and spermatogenesis. A third biological theme begins to emerge from recent insights into the Y chromosome's roles beyond the reproductive tract--a theme that promises to broaden the reach of Y-chromosome research by shedding light on fundamental sex differences in human health and disease.

Akt-mediated phosphorylation controls the activity of the Y-box protein MSY3 in skeletal muscle.

BACKGROUND: The Y-box protein MSY3/Csda represses myogenin transcription in skeletal muscle by binding a highly conserved cis-acting DNA element located just upstream of the myogenin minimal promoter (myogHCE). It is not known how this MSY3 activity is controlled in skeletal muscle. In this study, we provide multiple lines of evidence showing that the post-translational phosphorylation of MSY3 by Akt kinase modulates the MSY3 repression of myogenin. METHODS: Skeletal muscle and myogenic C2C12 cells were used to study the effects of MSY3 phosphorylation in vivo and in vitro on its sub-cellular localization and activity, by blocking the IGF1/PI3K/Akt pathway, by Akt depletion and over-expression, and by mutating potential MSY3 phosphorylation sites. RESULTS: We observed that, as skeletal muscle progressed from perinatal to postnatal and adult developmental stages, MSY3 protein became gradually dephosphorylated and accumulated in the nucleus. This correlated well with the reduction of phosphorylated active Akt. In C2C12 myogenic cells, blocking the IGF1/PI3K/Akt pathway using LY294002 inhibitor reduced MSY3 phosphorylation levels resulting in its accumulation in the nuclei. Knocking down Akt expression increased the amount of dephosphorylated MSY3 and reduced myogenin expression and muscle differentiation. MSY3 phosphorylation by Akt in vitro impaired its binding at the MyogHCE element, while blocking Akt increased MSY3 binding activity. While Akt over-expression rescued myogenin expression in MSY3 overexpressing myogenic cells, ablation of the Akt substrate, (Ser126 located in the MSY3 cold shock domain) promoted MSY3 accumulation in the nucleus and abolished this rescue. Furthermore, forced expression of Akt in adult skeletal muscle induced MSY3 phosphorylation and myogenin derepression. CONCLUSIONS: These results support the hypothesis that MSY3 phosphorylation by Akt interferes with MSY3 repression of myogenin circuit activity during muscle development. This study highlights a previously undescribed Akt-mediated signaling pathway involved in the repression of myogenin expression in myogenic cells and in mature muscle. Given the significance of myogenin regulation in adult muscle, the Akt/MSY3/myogenin regulatory circuit is a potential therapeutic target to counteract muscle degenerative disease.

Phylogeographic Refinement and Large Scale Genotyping of Human Y Chromosome Haplogroup E Provide New Insights into the Dispersal of Early Pastoralists in the African Continent.

Haplogroup E, defined by mutation M40, is the most common human Y chromosome clade within Africa. To increase the level of resolution of haplogroup E, we disclosed the phylogenetic relationships among 729 mutations found in 33 haplogroup DE Y-chromosomes sequenced at high coverage in previous studies. Additionally, we dissected the E-M35 subclade by genotyping 62 informative markers in 5,222 samples from 118 worldwide populations. The phylogeny of haplogroup E showed novel features compared with the previous topology, including a new basal dichotomy. Within haplogroup E-M35, we resolved all the previously known polytomies and assigned all the E-M35* chromosomes to five new different clades, all belonging to a newly identified subhaplogroup (E-V1515), which accounts for almost half of the E-M35 chromosomes from the Horn of Africa. Moreover, using a Bayesian phylogeographic analysis and a single nucleotide polymorphism-based approach we localized and dated the origin of this new lineage in the northern part of the Horn, about 12 ka. Time frames, phylogenetic structuring, and sociogeographic distribution of E-V1515 and its subclades are consistent with a multistep demic spread of pastoralism within north-eastern Africa and its subsequent diffusion to subequatorial areas. In addition, our results increase the discriminative power of the E-M35 haplogroup for use in forensic genetics through the identification of new ancestry-informative markers.

A medium-term, stochastic forecast model to support sustainable, mixed fisheries management in the Mediterranean Sea.

A medium-term (10 year) stochastic forecast model is developed and presented for mixed fisheries that can provide estimations of age-specific parameters for a maximum of 10 stocks and 10 fisheries. Designed to support fishery managers dealing with complex, multi-annual management plans, the model can be used to quantitatively test the consequences of various stock-specific and fishery-specific decisions, using non-equilibrium stock dynamics. Such decisions include fishing restrictions and other strategies aimed at achieving sustainable mixed fisheries consistent with the concept of maximum sustainable yield (MSY). In order to test the model, recently gathered data on seven stocks and four fisheries operating in the Ligurian and North Tyrrhenian Seas are used to generate quantitative, 10 year predictions of biomass and catch trends under four different management scenarios. The results show that using the fishing mortality at MSY as the biological reference point for the management of all stocks would be a strong incentive to reduce the technical interactions among concurrent fishing strategies. This would optimize the stock-specific exploitation and be consistent with sustainability criteria.