The transcriptome of early chicken embryos reveals signaling pathways governing rapid asymmetric cellularization and lineage segregation

A maternal dorsoventral prepattern revealed by an asymmetric distribution of ventralizing molecules before fertilization in Xenopus laevis

The establishment of the embryonic dorsoventral axis in Xenopus occurs when the radial symmetry around the egg's animal-vegetal axis is broken to give rise to the typical symmetry of Bilaterians. We have previously shown that the Notch1 protein is ventrally enriched during early embryogenesis in Xenopus laevis and zebrafish and exerts ventralizing activity through β-Catenin destabilization and the positive regulation of ventral center genes in X. laevis. These findings led us to further investigate when these asymmetries arise. In this work, we show that the asymmetrical distribution of Notch1 protein and mRNA precedes cortical rotation and even fertilization in X. laevis. Moreover, we found that in unfertilized eggs transcripts encoded by the ventralizing gene bmp4 are also asymmetrically distributed in the animal hemisphere and notch1 transcripts accumulate consistently on the same side of the eccentric maturation point. Strikingly, a Notch1 asymmetry orthogonal to the animal-vegetal axis appears during X. laevis oogenesis. Thus, we show for the first time a maternal bias in the distribution of molecules that are later involved in ventral patterning during embryonic axialization, strongly supporting the hypothesis of a dorsoventral prepattern or intrinsic bilaterality of Xenopus eggs before fertilization.

Comprehensive genome‑wide analysis of the chicken heat shock protein family: identification, genomic organization, and expression profiles in indigenous chicken with highly pathogenic avian influenza infection

BACKGROUND

Heat shock proteins (HSPs) function as molecular chaperones with critical roles in chicken embryogenesis, immune response to infectious diseases, and response to various environmental stresses. However, little is known on HSP genes in chicken. In this study, to understand the roles of chicken HSPs, we performed genome-wide identification, expression, and functional analyses of the HSP family genes in chicken.

RESULTS

A total of 76 HSP genes were identified in the chicken genome, which were further classified into eight distinct groups (I-VIII) based on phylogenetic tree analysis. The gene-structure analysis revealed that the members of each clade had the same or similar exon-intron structures. Chromosome mapping suggested that HSP genes were widely dispersed across the chicken genome, except in chromosomes 16, 18, 22, 25, 26, and 28-32, which lacked chicken HSP genes. On the other hand, the interactions among chicken HSPs were limited, indicating that the remaining functions of HSPs could be investigated in chicken. Moreover, KEGG pathway analysis showed that the HSP gene family was involved in the regulation of heat stress, apoptotic, intracellular signaling, and immune response pathways. Finally, RNA sequencing data revealed that, of the 76 chicken HSP genes, 46 were differentially expressed at 21 different growth stages in chicken embryos, and 72 were differentially expressed on post-infection day 3 in two indigenous Ri chicken lines infected with highly pathogenic avian influenza.

CONCLUSIONS

This study provides significant insights into the potential functions of HSPs in chicken, including the regulation of apoptosis, heat stress, chaperone activity, intracellular signaling, and immune response to infectious diseases.

Genome‑wide identification, organization, and expression profiles of the chicken fibroblast growth factor genes in public databases and Vietnamese indigenous Ri chickens against highly pathogenic avian influenza H5N1 virus infection

OBJECTIVE

Fibroblast growth factors (FGFs) play critical roles in embryo development, and immune responses to infectious diseases. In this study, to investigate the roles of FGFs, we performed genome-wide identification, expression, and functional analyses of FGF family members in chickens.

METHODS

Chicken FGFs genes were identified and analyzed by using bioinformatics approach. Expression profiles and Hierarchical cluster analysis of the FGFs genes in different chicken tissues were obtained from the genome-wide RNA-seq.

RESULTS

A total of 20 FGF genes were identified in the chicken genome, which were classified into seven distinct groups (A-F) in the phylogenetic tree. Gene structure analysis revealed that members of the same clade had the same or similar exon-intron structure. Chromosome mapping suggested that FGF genes were widely dispersed across the chicken genome and were located on chromosomes 1, 4-6, 9-10, 13, 15, 28, and Z. In addition, the interactions among FGF proteins and between FGFs and mitogen‑activated protein kinase (MAPK) proteins are limited, indicating that the remaining functions of FGF proteins should be further investigated in chickens. Kyoto encyclopedia of genes and genomes pathway analysis showed that FGF gene interacts with MAPK genes and are involved in stimulating signaling pathway and regulating immune responses. Furthermore, this study identified 15 differentially expressed genes (DEG) in 21 different growth stages during early chicken embryo development. RNA-sequencing data identified the DEG of FGFs on 1- and 3-days post infection in two indigenous Ri chicken lines infected with the highly pathogenic avian influenza virus H5N1 (HPAIV). Finally, all the genes examined through quantitative real-time polymerase chain reaction and RNA-Seq analyses showed similar responses to HPAIV infection in indigenous Ri chicken lines (R2 = 0.92- 0.95, p<0.01).

CONCLUSION

This study provides significant insights into the potential functions of FGFs in chickens, including the regulation of MAPK signaling pathways and the immune response of chickens to HPAIV infections.

Transcriptomic analysis reveals the dynamic changes of transcription factors during early development of chicken embryo

BACKGROUND

The transition from fertilized egg to embryo in chicken requires activation of hundreds of genes that were mostly inactivated before fertilization, which is accompanied with various biological processes. Undoubtedly, transcription factors (TFs) play important roles in regulating the changes in gene expression pattern observed at early development. However, the contribution of TFs during early embryo development of chicken still remains largely unknown that need to be investigated. Therefore, an understanding of the development of vertebrates would be greatly facilitated by study of the dynamic changes in transcription factors during early chicken embryo.

RESULTS

In the current study, we selected five early developmental stages in White Leghorn chicken, gallus gallus, for transcriptome analysis, cover 17,478 genes with about 807 million clean reads of RNA-sequencing. We have compared global gene expression patterns of consecutive stages and noted the differences. Comparative analysis of differentially expressed TFs (FDR < 0.05) profiles between neighboring developmental timepoints revealed significantly enriched biological categories associated with differentiation, development and morphogenesis. We also found that Zf-C2H2, Homeobox and bHLH were three dominant transcription factor families that appeared in early embryogenesis. More importantly, a TFs co-expression network was constructed and 16 critical TFs were identified.

CONCLUSION

Our findings provide a comprehensive regulatory framework of TFs in chicken early embryo, revealing new insights into alterations of chicken embryonic TF expression and broadening better understanding of TF function in chicken embryogenesis.

Identification of chicken LOC420478 as Bucky ball equivalent and potential germ plasm organizer in birds

Bucky ball was identified as germ plasm organizer in zebrafish and has proven crucial for Balbiani body condensation. A synteny comparison identified an uncharacterized gene locus in the chicken genome as predicted avian counterpart. Here, we present experimental evidence that this gene locus indeed encodes a ‘Bucky ball’ equivalent in matured oocytes and early embryos of chicken. Heterologous expression of Bucky ball fusion proteins both from zebrafish and chicken with a fluorescent reporter revealed unique patterns indicative for liquid–liquid phase separation of intrinsically disordered proteins. Immuno-labeling detected Bucky ball from oocytes to blastoderms with diffuse distribution in matured oocytes, aggregation in first cleavage furrows, and co-localization to the chicken vasa homolog (CVH). Later, Bucky ball translocated to the cytoplasm of first established cells, and showed nuclear translocation during the major zygotic activation together with CVH. Remarkably, during the phase of area pellucida formation, Bucky ball translocated back into the cytoplasm at stage EGK VI, whereas CVH remained within the nuclei. The condensation of Bucky ball and co-localization with CVH in cleavage furrows and nuclei of the centrally located cells strongly suggests chicken Bucky ball as a germ plasm organizer in birds, and indicate a special importance of the major zygotic activation for germline specification.

Female Germ Cell Development in Chickens and Humans: The Chicken Oocyte Enriched Genes Convergent and Divergent with the Human Oocyte

The development of germ cells and other physiological events in the differentiated ovary of humans are highly conserved with several mammalian species, except for the differences in timing. However, comparative knowledge on this topic is very scarce with respect to humans and lower vertebrates, such as chickens. In chickens, female germ cells enter into meiosis around embryonic day (E) 15.5 and are arrested in meiotic prophase I as primary oocytes. The oocytes arrested in meiosis I are accumulated in germ-cell cysts; shortly after hatching, they are enclosed by flattened granulosa cells in order to form primordial follicles. In humans, the process of meiotic recombination in female germ cells begins in the 10–11th week of gestation, and primordial follicles are formed at around week 20. In this review, we comprehensively elucidate both the conservation and the species-specific differences between chickens and humans with respect to germ cell, oocyte, and follicle development. Importantly, we provide functional insights into a set of chicken oocyte enriched genes (from E16 to 1 week post-hatch) that show convergent and divergent expression patterns with respect to the human oocyte (from week 11 to 26).

Chicken blastoderms and primordial germ cells possess a higher expression of DNA repair genes and lower expression of apoptosis genes to preserve their genome stability

DNA is susceptible to damage by various sources. When the DNA is damaged, the cell repairs the damage through an appropriate DNA repair pathway. When the cell fails to repair DNA damage, apoptosis is initiated. Although several genes are involved in five major DNA repair pathways and two major apoptosis pathways, a comprehensive understanding of those gene expression is not well-understood in chicken tissues. We performed whole-transcriptome sequencing (WTS) analysis in the chicken embryonic fibroblasts (CEFs), stage X blastoderms, and primordial germ cells (PGCs) to uncover this deficiency. Stage X blastoderms mostly consist of undifferentiated progenitor (pluripotent) cells that have the potency to differentiate into all cell types. PGCs are also undifferentiated progenitor cells that later differentiate into male and female germ cells. CEFs are differentiated and abundant somatic cells. Through WTS analysis, we identified that the DNA repair pathway genes were expressed more highly in blastoderms and high in PGCs than CEFs. Besides, the apoptosis pathway genes were expressed low in blastoderms and PGCs than CEFs. We have also examined the WTS-based expression profiling of candidate pluripotency regulating genes due to the conserved properties of blastoderms and PGCs. In the results, a limited number of pluripotency genes, especially the core transcriptional network, were detected higher in both blastoderms and PGCs than CEFs. Next, we treated the CEFs, blastoderm cells, and PGCs with hydrogen peroxide (H₂O₂) for 1 h to induce DNA damage. Then, the H₂O₂ treated cells were incubated in fresh media for 3–12 h to observe DNA repair. Subsequent analyses in treated cells found that blastoderm cells and PGCs were more likely to undergo apoptosis along with the loss of pluripotency and less likely to undergo DNA repair, contrasting with CEFs. These properties of blastoderms and PGCs should be necessary to preserve genome stability during the development of early embryos and germ cells, respectively.

Transcriptomic Analysis of Laying Hens Revealed the Role of Aging-Related Genes during Forced Molting

Molting in birds provides us with an ideal genetic model for understanding aging and rejuvenation since birds present younger characteristics for reproduction and appearance after molting. Forced molting (FM) by fasting in chickens causes aging of their reproductive system and then promotes cell redevelopment by providing water and feed again. To reveal the genetic mechanism of rejuvenation, we detected blood hormone indexes and gene expression levels in the hypothalamus and ovary of hens from five different periods during FM. Three hormones were identified as participating in FM. Furthermore, the variation trends of gene expression levels in the hypothalamus and ovary at five different stages were found to be basically similar using transcriptome analysis. Among them, 45 genes were found to regulate cell aging during fasting stress and 12 genes were found to promote cell development during the recovery period in the hypothalamus. In addition, five hub genes (INO80D, HELZ, AGO4, ROCK2, and RFX7) were identified by WGCNA. FM can restart the reproductive function of aged hens by regulating expression levels of genes associated with aging and development. Our study not only enriches the theoretical basis of FM but also provides insights for the study of antiaging in humans and the conception mechanism in elderly women.

Cellular evolution of language

The evolutionary origin of language remains unknown despite many efforts to determine the origin of this signature human trait. Based on epigenetic inheritance, the current article hypothesizes that language evolved from cell-cell communication as the basis for generating structure and function embryologically and phylogenetically, as did all physiologic traits. Beginning with lipids forming the first micelle, a vertical integration of the evolved properties of the cell, from multicellular organisms to the introduction of cholesterol into the cell membrane, to the evolution of the peroxisome, the water-land transition and duplication of the βAdrenergic Receptor, the evolution of endothermy, leading to bipedalism, freeing the forelimbs for toolmaking and language, selection pressure for myelinization of the central nervous system to facilitate calcium flux, bespeaks human expression, culminating in the evolution of civilization. This process is epitomized by the Area of Broca as the structural-functional site for both motor control and language formation. The mechanistic interrelationship between motor control and language formation is underscored by the role of FoxP2 gene expression in both bipedalism and language. The effect of endothermy on bipedalism, freeing the forelimbs for toolmaking and language as the vertical integration from Cosmology to Physiology as the basis for language bespeaks human expression.

Ultrafast functional profiling of RNA-seq data for nonmodel organisms

Computational time and cost remain a major bottleneck for RNA-seq data analysis of nonmodel organisms without reference genomes. To address this challenge, we have developed Seq2Fun, a novel, all-in-one, ultrafast tool to directly perform functional quantification of RNA-seq reads without transcriptome de novo assembly. The pipeline starts with raw read quality control: sequencing error correction, removing poly(A) tails, and joining overlapped paired-end reads. It then conducts a DNA-to-protein search by translating each read into all possible amino acid fragments and subsequently identifies possible homologous sequences in a well-curated protein database. Finally, the pipeline generates several informative outputs including gene abundance tables, pathway and species hit tables, an HTML report to visualize the results, and an output of clean reads annotated with mapped genes ready for downstream analysis. Seq2Fun does not have any intermediate steps of file writing and loading, making I/O very efficient. Seq2Fun is written in C++ and can run on a personal computer with a limited number of CPUs and memory. It can process >2,000,000 reads/min and is >120 times faster than conventional workflows based on de novo assembly, while maintaining high accuracy in our various test data sets.

Dynamic changes of transcriptome of fifth-instar spodoptera litura larvae in response to insecticide

Spodoptera litura is a major insect with a cosmopolitan distribution and strong resistance to multiple insecticides. Determining the molecular basis and key candidate genes of the insecticide resistance of S. litura may help in managing this insect. In this study, fifth-instar S. litura larvae were subjected to transcriptome analysis at 6, 12, 24, 48, and 72 h after feeding on an LC20 dose of avermectin. The result showed that genes responding to avermectin changed dynamically with different gene counts and resistance mechanisms at the fifth instar based on a metabolic pathway map. These responses included degrading the insecticide by a series of P450 and glutathione-S-transferase enzymes starting at the 12 h time point, with subsequent increases in the number of genes involved and shifts to TOLL and immune deficiency (IMD) pathways at 48 h after feeding the insecticide. Weighted correlation network analysis (WGCNA) determined a co-expression module related to the avermectin response at 12 and 24 h (r = 0.403, p = 0.0371; r = 0.436, p = 0.023), in which a hub gene (LOC111358940) related to metalloproteinase activity was identified. In addition, Analysis of the genes in the co-expression module further revealed that eight genes encoding UDP-glucuronosyltransferases were directly associated with insecticide response in S. litura. These results provide better understanding of the avermectin response mechanism of S. litura and may be useful in developing improved control strategies for this species.

SUPPLEMENTARY INFORMATION

The online version of this article (10.1007/s13205-021-02651-9) contains supplementary material, which is available to authorized users.

Whole-Transcriptome Sequencing-Based Analysis of DAZL and Its Interacting Genes during Germ Cells Specification and Zygotic Genome Activation in Chickens

The deleted in azoospermia like (DAZL) is required for germ cells development and maintenance. In chickens, the mRNA and protein of DAZL, a representative maternally inherited germ plasm factor, are detected in the germ plasm of oocyte, zygote, and all stages of the intrauterine embryos. However, it is still insufficient to explain the origin and specification process of chicken germ cells, because the stage at which the zygotic transcription of DAZL occurs and the stage at which the maternal DAZL RNA/protein clears have not yet been fully identified. Moreover, a comprehensive understanding of the expression of DAZL interacting genes during the germ cells specification and development and zygotic genome activation (ZGA) is lacking in chickens. In this study, we identified a set of DAZL interacting genes in chickens using in silico prediction method. Then, we analyzed the whole-transcriptome sequencing (WTS)-based expression of DAZL and its interacting genes in the chicken oocyte, zygote, and Eyal-Giladi and Kochav (EGK) stage embryos (EGK.I to EGK.X). In the results, DAZL transcripts are increased in the zygote (onset of transcription), maintained the increased level until EGK.VI, and decreased from EGK.VIII (possible clearance of maternal RNAs). Among the DAZL interacting genes, most of them are increased either at 1st ZGA or 2nd ZGA, indicating their involvement in germ cells specification and development.

Endodermal Maternal Transcription Factors Establish Super-Enhancers during Zygotic Genome Activation

Elucidation of the sequence of events underlying the dynamic interaction between transcription factors and chromatin states is essential. Maternal transcription factors function at the top of the regulatory hierarchy to specify the primary germ layers at the onset of zygotic genome activation (ZGA). We focus on the formation of endoderm progenitor cells and examine the interactions between maternal transcription factors and chromatin state changes underlying the cell specification process. Endoderm-specific factors Otx1 and Vegt together with Foxh1 orchestrate endoderm formation by coordinated binding to select regulatory regions. These interactions occur before the deposition of enhancer histone marks around the regulatory regions, and these TFs recruit RNA polymerase II, regulate enhancer activity, and establish super-enhancers associated with important endodermal genes. Therefore, maternal transcription factors Otx1, Vegt, and Foxh1 combinatorially regulate the activity of super-enhancers, which in turn activate key lineage-specifying genes during ZGA.

How do maternal transcription factors interact with chromatin regions to coordinate the endodermal gene regulatory program? Paraiso et al. demonstrate that combinatorial binding of maternal Otx1, Vegt, and Foxh1 to select enhancers and super-enhancers in the genome controls endodermal cell fate specification during zygotic gene activation.

Current Approaches and Applications in Avian Genome Editing

Advances in genome-editing technologies and sequencing of animal genomes enable researchers to generate genome-edited (GE) livestock as valuable animal models that benefit biological researches and biomedical and agricultural industries. As birds are an important species in biology and agriculture, their genome editing has gained significant interest and is mainly performed by using a primordial germ cell (PGC)-mediated method because pronuclear injection is not practical in the avian species. In this method, PGCs can be isolated, cultured, genetically edited in vitro, and injected into a recipient embryo to produce GE offspring. Recently, a couple of GE quail have been generated by using the newly developed adenovirus-mediated method. Without technically required in vitro procedures of the PGC-mediated method, direct injection of adenovirus into the avian blastoderm in the freshly laid eggs resulted in the production of germ-line chimera and GE offspring. As more approaches are available in avian genome editing, avian research in various fields will progress rapidly. In this review, we describe the development of avian genome editing and scientific and industrial applications of GE avian species.

Zygotic genome activation in the chicken: a comparative review

Maternal RNAs and proteins in the oocyte contribute to early embryonic development. After fertilization, these maternal factors are cleared and embryonic development is determined by an individual's own RNAs and proteins, in a process called the maternal-to-zygotic transition. Zygotic transcription is initially inactive, but is eventually activated by maternal transcription factors. The timing and molecular mechanisms involved in zygotic genome activation (ZGA) have been well-described in many species. Among birds, a transcriptome-based understanding of ZGA has only been explored in chickens by RNA sequencing of intrauterine embryos. RNA sequencing of chicken intrauterine embryos, including oocytes, zygotes, and Eyal-Giladi and Kochav (EGK) stages I-X has enabled the identification of differentially expressed genes between consecutive stages. These studies have revealed that there are two waves of ZGA: a minor wave at the one-cell stage (shortly after fertilization) and a major wave between EGK.III and EGK.VI (during cellularization). In the chicken, the maternal genome is activated during minor ZGA and the paternal genome is quiescent until major ZGA to avoid transcription from supernumerary sperm nuclei. In this review, we provide a detailed overview of events in intrauterine embryonic development in birds (and particularly in chickens), as well as a transcriptome-based analysis of ZGA.

Evolutionary History of GLIS Genes Illuminates Their Roles in Cell Reprograming and Ciliogenesis

The GLIS family transcription factors, GLIS1 and GLIS3, potentiate generation of induced pluripotent stem cells (iPSCs). In contrast, another GLIS family member, GLIS2, suppresses cell reprograming. To understand how these disparate roles arose, we examined evolutionary origins and genomic organization of GLIS genes. Comprehensive phylogenetic analysis shows that GLIS1 and GLIS3 originated during vertebrate whole genome duplication, whereas GLIS2 is a sister group to the GLIS1/3 and GLI families. This result is consistent with their opposing functions in cell reprograming. Glis1 evolved faster than Glis3, losing many protein-interacting motifs. This suggests that Glis1 acquired new functions under weakened evolutionary constraints. In fact, GLIS1 induces induced pluripotent stem cells more strongly. Transcriptomic data from various animal embryos demonstrate that glis1 is maternally expressed in some tetrapods, whereas vertebrate glis3 and invertebrate glis1/3 genes are rarely expressed in oocytes, suggesting that vertebrate (or tetrapod) Glis1 acquired a new expression domain and function as a maternal factor. Furthermore, comparative genomic analysis reveals that glis1/3 is part of a bilaterian-specific gene cluster, together with rfx3, ndc1, hspb11, and lrrc42. Because known functions of these genes are related to cilia formation and function, the last common ancestor of bilaterians may have acquired this cluster by shuffling gene order to establish more sophisticated epithelial tissues involving cilia. This evolutionary study highlights the significance of GLIS1/3 for cell reprograming, development, and diseases in ciliated organs such as lung, kidney, and pancreas.

Bioinformatic identification and expression analysis of the chicken B cell lymphoma (BCL) gene

BACKGROUND

B cell lymphoma (BCL) families play an important role in apoptosis as a growth factor, cell death programming, cytokine expression and immune-related genes expression.

OBJECTIVES

In this study, to investigate the roles of BCLs, we performed genome-wide identification, expression and functional analyses of the BCL family in chicken.

METHODS

Chicken BCLs genes were identified and analyzed by using bioinformatics approach. Expression profiles and Hierarchical cluster analysis of the BCLs genes in different chicken tissues were obtained from the genome-wide RNA-seq in the GEO, and Cluster and Java Treeview, respectively.

RESULTS

A total of 16 BCLs genes were identified from the chicken genome, which could be further classified into five distinct groups in the phylogenetic tree. On the other hand, the interaction among BCLs proteins and between BCLs proteins with NF-κB subunits are limited, indicating that the remaining the functions of BCLs protein could be investigated in chicken. Moreover, KEGG pathway analysis indicated that BCL gene family was involved in regulation of apoptotic and immune response. Finally, BCL gene family was differentially expressed in chicken tissues, pathogen infection and growth stages of early chicken early embryo.

CONCLUSION

This study provides significant insights into the potential functions of BCLs in chicken, including the regulation of apoptosis, cell death and expression of immune-related genes.

Zygotic gene activation in the chicken occurs in two waves, the first involving only maternally derived genes

The first wave of transcriptional activation occurs after fertilisation in a species-specific pattern. Despite its importance to initial embryonic development, the characteristics of transcription following fertilisation are poorly understood in Aves. Here, we report detailed insights into the onset of genome activation in chickens. We established that two waves of transcriptional activation occurred, one shortly after fertilisation and another at Eyal-Giladi and Kochav Stage V. We found 1544 single nucleotide polymorphisms across 424 transcripts derived from parents that were expressed in offspring during the early embryonic stages. Surprisingly, only the maternal genome was activated in the zygote, and the paternal genome remained silent until the second-wave, regardless of the presence of a paternal pronucleus or supernumerary sperm in the egg. The identified maternal genes involved in cleavage that were replaced by bi-allelic expression. The results demonstrate that only maternal alleles are activated in the chicken zygote upon fertilisation, which could be essential for early embryogenesis and evolutionary outcomes in birds.

The early stages of animal development involve a handover of genetic control. Initially, the egg cell is maintained by genetic information inherited from the mother, but soon after fertilization it starts to depend on its own genes instead. Activating genes inside the fertilized egg cell (zygote) so that they can take control of development is known as zygotic genome activation.

Despite the fact that birds are often used to study how embryos develop, zygotic genome activation in birds is not well understood. Fertilization in birds, including chickens, is different to mammals in that it requires multiple sperm to fertilize an egg cell. As such, zygotic genome activation in birds is likely to differ from that in mammals.

By examining gene expression in embryos from mixed-breed chickens, Hwang, Seo et al. showed that there are two stages of zygotic genome activation in chickens. The genes derived from the mother become active in the first stage, while genes from the father become active in the second stage. Genome activation in birds is therefore very different to the same process in mammals, which involves genome activation of both parents from the first stage. This extra level of control may help to prevent genetic complications resulting from the presence of multiple sperm, each of which carries a different set of genes from the father.