Duplicated Gene Families Research Articles

A genome assembly and transcriptome atlas of the inbred Babraham pig to illuminate porcine immunogenetic variation.

The inbred Babraham pig serves as a valuable biomedical model for research due to its high level of homozygosity, including in the major histocompatibility complex (MHC) loci and likely other important immune-related gene complexes, which are generally highly diverse in outbred populations. As the ability to control for this diversity using inbred organisms is of great utility, we sought to improve this resource by generating a long-read whole genome assembly and transcriptome atlas of a Babraham pig. The genome was de novo assembled using PacBio longreads and error-corrected using Illumina shortreads. Assembled contigs were then mapped to the porcine reference assembly, Sscrofa11.1, to generate chromosome-level scaffolds. The resulting TPI_Babraham_pig_v1 assembly is nearly as contiguous as Sscrofa11.1 with a contig N50 of 34.95 Mb and contig L50 of 23. The remaining sequence gaps are generally the result of poor assembly across large and highly repetitive regions such as the centromeres and tandemly duplicated gene families, including immune-related gene complexes, that often vary in gene content between haplotypes. We also further confirm homozygosity across the Babraham MHC and characterize the allele content and tissue expression of several other immune-related gene complexes, including the antibody and T cell receptor loci, the natural killer complex, and the leukocyte receptor complex. The Babraham pig genome assembly provides an alternate highly contiguous porcine genome assembly as a resource for the livestock genomics community. The assembly will also aid biomedical and veterinary research that utilizes this animal model such as when controlling for genetic variation is critical.

Structural polymorphism and diversity of human segmental duplications.

Segmental duplications (SDs) contribute significantly to human disease, evolution, and diversity yet have been difficult to resolve at the sequence level. We present a population genetics survey of SDs by analyzing 170 human genome assemblies where the majority of SDs are fully resolved using long-read sequence assembly. Excluding the acrocentric short arms, we identify 173.2 Mbp of duplicated sequence (47.4 Mbp not present in the telomere-to-telomere reference) distinguishing fixed from structurally polymorphic events. We find that intrachromosomal SDs are among the most variable with rare events mapping near their progenitor sequences. African genomes harbor significantly more intrachromosomal SDs and are more likely to have recently duplicated gene families with higher copy number when compared to non-African samples. A comparison to a resource of 563 million full-length Iso-Seq reads identifies 201 novel, potentially protein-coding genes corresponding to these copy number polymorphic SDs.

Halophytes and heavy metals: A multi-omics approach to understand the role of gene and genome duplication in the abiotic stress tolerance of Cakile maritima.

The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole-genome triplication with closely related salt-sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy. Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium. Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance. These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response.

Chromosomal inversion polymorphisms shape human brain morphology

The impact of chromosomal inversions on human brain morphology remains underexplored. We studied 35 common inversions classified from genotypes of 33,018 adults with European ancestry. The inversions at 2p22.3, 16p11.2, and 17q21.31 reach genome-wide significance, followed by 8p23.1 and 6p21.33, in their association with cortical and subcortical morphology. The 17q21.31, 8p23.1, and 16p11.2 regions comprise the LRRC37, OR7E, and NPIP duplicated gene families. We find the 17q21.31 MAPT inversion region, known for harboring neurological risk, to be the most salient locus among common variants for shaping and patterning the cortex. Overall, we observe the inverted orientations decreasing brain size, with the exception that the 2p22.3 inversion is associated with increased subcortical volume and the 8p23.1 inversion is associated with increased motor cortex. These significant inversions are in the genomic hotspots of neuropsychiatric loci. Our findings are generalizable to 3,472 children and demonstrate inversions as essential genetic variation to understand human brain phenotypes.

Phylogenetic placement of whole genome duplications in yeasts through quantitative analysis of hierarchical orthologous groups

Background: Whole genome duplications (WGD) are genomic events leading to formation of polyploid organisms. Resulting duplicated genes play important roles in driving species evolution and diversification. After such events, the initial ploidy is usually restored, complicating their detection across evolution. With the advance of bioinformatics and the rising number of new well-assembled genomes, new detection methods are ongoingly being developed to overcome the weaknesses of different approaches. Results: Here we propose a novel method for detecting WGD in yeast lineages based on the quantitative and comparative analysis of hierarchical orthologous groups (HOGs) of duplicated genes for a given set of organisms. We reconstruct ancestral genomes to obtain evolutionary information for each phylogenetic branch. This reconstruction relies on the inference of HOGs from the selected species’ proteomes. To estimate WGD events, the number of HOGs of duplicated genes across all taxonomic ranges are adjusted according to the molecular clock hypothesis and by the average genome size. Branches with a significant increase in the adjusted number of duplicated gene families are kept as candidates for WGD placement. The developed method was tested on two real datasets and showed promising results in phylogenetic WGD placements on the yeast lineage.

Expanding duplication of the testis PHD Finger Protein 7 (PHF7) gene in the chicken genome

Gene duplications increase genetic and phenotypic diversity and occur in complex genomic regions that are still difficult to sequence and assemble. PHD Finger Protein 7 (PHF7) acts during spermiogenesis for histone-to-histone protamine exchange and is a determinant of male fertility in Drosophila and the mouse. We aimed to explore and characterise in the chicken genome the expanding family of the numerous orthologues of the unique mouse Phf7 gene (highly expressed in the testis), observing the fact that this information is unclear and/or variable according to the versions of databases. We validated nine primer pairs by in silico PCR for their use in screening the chicken bacterial artificial chromosome (BAC) library to produce BAC-derived probes to detect and localise PHF7-like loci by fluorescence in situ hybridisation (FISH). We selected nine BAC that highlighted nine chromosomal regions for a total of 10 distinct PHF7-like loci on five Gallus gallus chromosomes: Chr1 (three loci), Chr2 (two loci), Chr12 (one locus), Chr19 (one locus) and ChrZ (three loci). We sequenced the corresponding BAC by using high-performance PacBio technology. After assembly, we performed annotation with the FGENESH program: there were a total of 116 peptides, including 39 PHF7-like proteins identified by BLASTP. These proteins share a common exon-intron core structure of 8–11 exons. Phylogeny revealed that the duplications occurred first between chromosomal regions and then inside each region. There are other duplicated genes in the identified BAC sequences, suggesting that these genomic regions exhibit a high rate of tandem duplication. We showed that the PHF7 gene, which is highly expressed in the rooster testis, is a highly duplicated gene family in the chicken genome, and this phenomenon probably concerns other bird species.

Sequence Diversity and Differences at the Highly Duplicated MHC-I Gene Reflect Viral Susceptibility in Sympatric Pinniped Species.

Differences in disease susceptibility among species can result from rapid host-pathogen coevolution and differences in host species ecology that affect the strength and direction of natural selection. Among 2 sympatric pinniped species that differ in sociality and putative disease exposure, we investigate observed differences in susceptibility through an analysis of a highly variable, duplicated gene family involved in the vertebrate immune response. Using high-throughput amplicon sequencing, we characterize diversity at the 2 exons that encode the peptide binding region of the major histocompatibility complex class I (MHC-I) gene in harbor (N = 60) and gray (N = 90) seal populations from the Northwest Atlantic. Across species, we identified 106 full-length exon 2 and 103 exon 3 sequence variants and a minimum of 11 duplicated MHC-I loci. The sequence variants clustered in 15 supertypes defined by the physiochemical properties of the peptide binding region, including a putatively novel Northwest Atlantic MHC-I diversity sublineage. Trans-species polymorphisms, dN/dS ratios, and evidence of gene conversion among supertypes are consistent with balancing selection acting on this gene. High functional redundancy suggests particularly strong selection among gray seals at the novel Northwest Atlantic MHC-I diversity sublineage. At exon 2, harbor seals had a significantly greater number of variants per individual than gray seals, but fewer supertypes. Supertype richness and private supertypes are hypothesized to contribute to observed differences in disease resistance between species, as consistently, across the North Atlantic and many disease outbreaks, gray seals appear to be more resistant to respiratory viruses than harbor seals.

High quality genome assembly of theanhydrobiotic midgeprovides insights on a single chromosome-based emergenceof extreme desiccation tolerance.

Non-biting midges (Chironomidae) are known to inhabit a wide range of environments, and certain species can tolerate extreme conditions, where the rest of insects cannot survive. In particular, the sleeping chironomid Polypedilum vanderplanki is known for the remarkable ability of its larvae to withstand almost complete desiccation by entering a state called anhydrobiosis. Chromosome numbers in chironomids are higher than in other dipterans and this extra genomic resource might facilitate rapid adaptation to novel environments. We used improved sequencing strategies to assemble a chromosome-level genome sequence for P. vanderplanki for deep comparative analysis of genomic location of genes associated with desiccation tolerance. Using whole genome-based cross-species and intra-species analysis, we provide evidence for the unique functional specialization of Chromosome 4 through extensive acquisition of novel genes. In contrast to other insect genomes, in the sleeping chironomid a uniquely high degree of subfunctionalization in paralogous anhydrobiosis genes occurs in this chromosome, as well as pseudogenization in a highly duplicated gene family. Our findings suggest that the Chromosome 4 in Polypedilum is a site of high genetic turnover, allowing it to act as a ‘sandbox’ for evolutionary experiments, thus facilitating the rapid adaptation of midges to harsh environments.

Machine learning discovery of missing links that mediate alternative branches to plant alkaloids

Engineering the microbial production of secondary metabolites is limited by the known reactions of correctly annotated enzymes. Therefore, the machine learning discovery of specialized enzymes offers great potential to expand the range of biosynthesis pathways. Benzylisoquinoline alkaloid production is a model example of metabolic engineering with potential to revolutionize the paradigm of sustainable biomanufacturing. Existing bacterial studies utilize a norlaudanosoline pathway, whereas plants contain a more stable norcoclaurine pathway, which is exploited in yeast. However, committed aromatic precursors are still produced using microbial enzymes that remain elusive in plants, and additional downstream missing links remain hidden within highly duplicated plant gene families. In the current study, machine learning is applied to predict and select plant missing link enzymes from homologous candidate sequences. Metabolomics-based characterization of the selected sequences reveals potential aromatic acetaldehyde synthases and phenylpyruvate decarboxylases in reconstructed plant gene-only benzylisoquinoline alkaloid pathways from tyrosine. Synergistic application of the aryl acetaldehyde producing enzymes results in enhanced benzylisoquinoline alkaloid production through hybrid norcoclaurine and norlaudanosoline pathways.

Functional buffering via cell-specific gene expression promotes tissue homeostasis and cancer robustness

Functional buffering that ensures biological robustness is critical for maintaining tissue homeostasis, organismal survival, and evolution of novelty. However, the mechanism underlying functional buffering, particularly in multicellular organisms, remains largely elusive. Here, we proposed that functional buffering can be mediated via expression of buffering genes in specific cells and tissues, by which we named Cell-specific Expression-BUffering (CEBU). We developed an inference index (C-score) for CEBU by computing C-scores across 684 human cell lines using genome-wide CRISPR screens and transcriptomic RNA-seq. We report that C-score-identified putative buffering gene pairs are enriched for members of the same duplicated gene family, pathway, and protein complex. Furthermore, CEBU is especially prevalent in tissues of low regenerative capacity (e.g., bone and neuronal tissues) and is weakest in highly regenerative blood cells, linking functional buffering to tissue regeneration. Clinically, the buffering capacity enabled by CEBU can help predict patient survival for multiple cancers. Our results suggest CEBU as a potential buffering mechanism contributing to tissue homeostasis and cancer robustness in humans.

New telomere to telomere assembly of human chromosome 8 reveals a previous underestimation of G-quadruplex forming sequences and inverted repeats.

Taking advantage of evolving and improving sequencing methods, human chromosome 8 is now available as a gapless, end-to-end assembly. Thanks to advances in long-read sequencing technologies, its centromere, telomeres, duplicated gene families and repeat-rich regions are now fully sequenced. We were interested to assess if the new assembly altered our understanding of the potential impact of non-B DNA structures within this completed chromosome sequence. It has been shown that non-B secondary structures, such as G-quadruplexes, hairpins and cruciforms, have important regulatory functions and potential as targeted therapeutics. Therefore, we analysed the presence of putative G-quadruplex forming sequences and inverted repeats in the current human reference genome (GRCh38) and in the new end-to-end assembly of chromosome 8. The comparison revealed that the new assembly contains significantly more inverted repeats and G-quadruplex forming sequences compared to the current reference sequence. This observation can be explained by improved accuracy of the new sequencing methods, particularly in regions that contain extensive repeats of bases, as is preferred by many non-B DNA structures. These results show a significant underestimation of the prevalence of non-B DNA secondary structure in previous assembly versions of the human genome and point to their importance being not fully appreciated. We anticipate that similar observations will occur as the improved sequencing technologies fill in gaps across the genomes of humans and other organisms.

Adaptive evolution driving the young duplications in six Rosaceae species

BackgroundIn plant genomes, high proportions of duplicate copies reveals that gene duplications play an important role in the evolutionary processes of plant species. A series of gene families under positive selection after recent duplication events in plant genomes indicated the evolution of duplicates driven by adaptive evolution. However, the genome-wide evolutionary features of young duplicate genes among closely related species are rarely reported.ResultsIn this study, we conducted a systematic survey of young duplicate genes at genome-wide levels among six Rosaceae species, whose whole-genome sequencing data were successively released in recent years. A total of 35,936 gene families were detected among the six species, in which 60.25% were generated by young duplications. The 21,650 young duplicate gene families could be divided into two expansion types based on their duplication patterns, species-specific and lineage-specific expansions. Our results showed the species-specific expansions advantaging over the lineage-specific expansions. In the two types of expansions, high-frequency duplicate domains exhibited functional preference in response to environmental stresses.ConclusionsThe functional preference of the young duplicate genes in both the expansion types showed that they were inclined to respond to abiotic or biotic stimuli. Moreover, young duplicate genes under positive selection in both species-specific and lineage-specific expansions suggested that they were generated to adapt to the environmental factors in Rosaceae species.

Global Evolutionary Analysis of 11 Gene Families Part of Reactive Oxygen Species (ROS) Gene Network in Four Eucalyptus Species.

Eucalyptus is a worldwide hard-wood species which increasingly focused on. To adapt to various biotic and abiotic stresses, Eucalyptus have evolved complex mechanisms, increasing the cellular concentration of reactive oxygen species (ROS) by numerous ROS controlling enzymes. To better analyse the ROS gene network and discuss the differences between four Eucalyptus species, ROS gene network including 11 proteins families (1CysPrx, 2CysPrx, APx, APx-R, CIII Prx, Diox, GPx, Kat, PrxII, PrxQ and Rboh) were annotated and compared in an expert and exhaustive manner from the genomic data available from E. camaldulensis, E. globulus, E. grandis, and E. gunnii. In addition, a specific sequencing strategy was performed in order to determine if the missed sequences in at least one organism are the results of gain/loss events or only sequencing gaps. We observed that the automatic annotation applied to multigenic families is the source of miss-annotation. Base on the family size, the 11 families can be categorized into duplicated gene families (CIII Prx, Kat, 1CysPrx, and GPx), which contain a lot of gene duplication events and non-duplicated families (APx, APx-R, Rboh, DiOx, 2CysPrx, PrxII, and PrxQ). The gene family sizes are much larger in Eucalyptus than most of other angiosperms due to recent gene duplications, which could give higher adaptability to environmental changes and stresses. The cross-species comparative analysis shows gene gain and loss events during the evolutionary process. The 11 families possess different expression patterns, while in the Eucalyptus genus, the ROS families present similar expression patterns. Overall, the comparative analysis might be a good criterion to evaluate the adaptation of different species with different characters, but only if data mining is as exhaustive as possible. It is also a good indicator to explore the evolutionary process.

The Reproducibility of an Inferred Tree and the Diploidization of Gene Segregation after Genome Duplication.

We previously introduced a numerical quantity called the stability (Ps) of an inferred tree and showed that for the tree to be reliable this stability as well as the reliability of the tree, which is usually computed as the bootstrap probability (Pb), must be high. However, if genome duplication occurs in a species, a gene family of the genome also duplicates, and for this reason alone some Ps values can be high in a tree of the duplicated gene families. In addition, the topology of the duplicated gene family can be similar to that of the original gene family if such gene families are identifiable. After genome duplication, however, the gene families are often partially deleted or partially duplicated, and the duplicated gene family may not show the same topology as that of the original family. It is therefore necessary to compute the similarity of the topologies of the duplicated and the original gene families. In this paper, we introduce another quantity called the reproducibility (Pr) for measuring the similarity of the two gene families. To show how to compute the Pr values, we first compute the Pb and Ps values for each of the MHC class II α and β chain gene families, which were apparently generated by genome duplication. We then compute the Pr values for the MHC class II α and β chain gene families. The Pr values for the α and β chain gene families are now low, and this suggests that the diploidization of gene segregation has occurred after the genome duplication. Currently higher animals, defined as animals with complex phenotypic characters, generally have a higher genome size, and this increase in genome size appears to have been caused by genome duplication and diploidization of gene segregation after genome duplication.

Accurate Prediction of RH Genotypes Using Whole Genome Sequencing Data

IntroductionRBC alloimmunization is common in patients with sickle cell disease (SCD). Despite serological matching RBCs for major Rh antigens, Rh alloimmunization remains problematic. The Rh blood group is encoded by two genes RHD and RHCE, which exhibit extensive nucleotide polymorphism and chromosome structural changes, resulting in the formation of Rh variant antigens. Rh variants can result in loss of protein epitopes or expression of neo-epitopes, and are common in SCD patients. Hence SCD patients harboring Rh variants can be predisposed to Rh alloimmunization. Given the limitation of traditional serologic antigen typing for detection of Rh variants, molecular genotyping has become required. A DNA microarray-based platform, BioArray RHCE and RHD BeadChip (Immuncor) is available for RH genotyping. However, it detects the most common, but not all, variants. Whole exome sequence data have been used for prediction of Rh variants (Chou, et. al, Blood Adv., 2017), offer some advantages, including detection of rare variants, structural rearrangements and copy number variation. However, whole genome sequence (WGS) analysis of RHD/RHCE is challenging due to difficulties in mapping next generation sequencing (NGS) reads to this duplicated gene family. We developed a computational algorithm to identify RH variants using WGS data.MethodsThe pipeline included three major components, RH allele database construction, RH variant calling, and classification of Rh blood group according the identified variants. The RH allele database was built based on NCBI Blood Group Antigen Gene Mutation (BGMUT) and International Society of Blood Transfusion (ISBT) database. Since the alleles in the BGMUT and ISBT databases were specified according to conventional RH genes (RHD, L08429; RHCE, DQ322275) that are different from those on reference human genome, we first called the variations based on the reference human genome. The positions of the identified variations were subsequently corrected to match with the BGMUT and ISBT annotation system. Next, the NGS reads with low base quality and/or mapping quality were discarded during the variation calling step. Synonymous and non-synonymous amino acid changes were characterized for each polymorphism. Haplotypes were constructed for the segments with NGS read support. Gene sequencing coverage was calculated to determine gene deletions or amplifications. Lastly, we implemented an algorithm to predict RH genotypes based on a selection of candidate alleles by read-mapping profile which considers both sequence variations and sequence consistency followed by a likelihood-based ranking of all pairwise combinations of the selected alleles. The allele combination with the highest likelihood is considered the most likely pair of alleles at a given locus. Patient specimens used in this study were from participants of the Sickle Cell Clinical Research and Intervention Program (SCCRIP, Hankins et al. Pediatr Blood Cancer. 2018).ResultsWe validated our method in a cohort of 58 SCD patients whose RH genotypes had been determined by BioArray RhCE and RhD BeadChip and supplementary molecular tests that identify the most common variants among individuals of African descent. In this validation cohort including a total of 11 RHD and 13 RHCE alleles, our approach achieved a concordance rate of 85.85% (91 of 106 alleles) for RHD and 83.02% (88 of 106 alleles) for RHCE genotyping. WGS was highly sensitive in distinguishing homozygosity from heterozygosity of genes. By comparing the numbers of NGS reads on RH regions and whole genome average coverage, heterozygous deletion can be determined. Since WGS provides comprehensive genotyping, our analysis identified single nucleotide polymorphisms that were not identified by the BeadChip and supplemental molecular testing. The final source of discordance was likely due to the short read length of NGS such that haplotype phases cannot be correctly predicted if the variations are separated by thousands of base pairs, for which long read DNA sequencing or RNA/cDNA sequencing are required. Evaluation of the identified discrepancies is ongoing.ConclusionsWe developed and validated a diagnostic method for RH genotyping that leveraged the accuracy and flexibility of RH genotyping based on WGS data. With further optimization of our method, this may be useful for RBC genotype matching sickle cell patients to blood donors in the future. DisclosuresHankins:Novartis: Research Funding; Global Blood Therapeutics: Research Funding; NCQA: Consultancy; bluebird bio: Consultancy.

Frequent nonallelic gene conversion on the human lineage and its effect on the divergence of gene duplicates

Gene conversion is the copying of a genetic sequence from a "donor" region to an "acceptor." In nonallelic gene conversion (NAGC), the donor and the acceptor are at distinct genetic loci. Despite the role NAGC plays in various genetic diseases and the concerted evolution of gene families, the parameters that govern NAGC are not well characterized. Here, we survey duplicate gene families and identify converted tracts in 46% of them. These conversions reflect a large GC bias of NAGC. We develop a sequence evolution model that leverages substantially more information in duplicate sequences than used by previous methods and use it to estimate the parameters that govern NAGC in humans: a mean converted tract length of 250 bp and a probability of [Formula: see text] per generation for a nucleotide to be converted (an order of magnitude higher than the point mutation rate). Despite this high baseline rate, we show that NAGC slows down as duplicate sequences diverge-until an eventual "escape" of the sequences from its influence. As a result, NAGC has a small average effect on the sequence divergence of duplicates. This work improves our understanding of the NAGC mechanism and the role that it plays in the evolution of gene duplicates.

A Fast and Exact Algorithm for the Exemplar Breakpoint Distance.

A fundamental problem in comparative genomics is to compute the distance between two genomes. For two genomes without duplicate genes, we can easily compute a variety of distance measures in linear time, but the problem is NP-hard under most models when genomes contain duplicate genes. Sankoff proposed the use of exemplars to tackle the problem of duplicate genes and gene families: each gene family is represented by a single gene (the exemplar for that family), chosen so as to optimize some metric. Unfortunately, choosing exemplars is itself an NP-hard problem. In this article, we propose a very fast and exact algorithm to compute the exemplar breakpoint distance, based on new insights in the underlying structure of genome rearrangements and exemplars. We evaluate the performance of our algorithm on simulation data and compare its performance to the best effort to date (a divide-and-conquer approach), showing that our algorithm runs much faster and scales much better. We also devise a new algorithm for the intermediate breakpoint distance problem, which can then be applied to assign orthologs. We compare our algorithm with the state-of-the-art method MSOAR by assigning orthologs among five well annotated mammalian genomes, showing that our algorithm runs much faster and is slightly more accurate than MSOAR.

Epigenetic silencing in transgenic plants.

Epigenetic silencing is a natural phenomenon in which the expression of genes is regulated through modifications of DNA, RNA, or histone proteins. It is a mechanism for defending host genomes against the effects of transposable elements and viral infection, and acts as a modulator of expression of duplicated gene family members and as a silencer of transgenes. A major breakthrough in understanding the mechanism of epigenetic silencing was the discovery of silencing in transgenic tobacco plants due to the interaction between two homologous promoters. The molecular mechanism of epigenetic mechanism is highly complicated and it is not completely understood yet. Two different molecular routes have been proposed for this, that is, transcriptional gene silencing, which is associated with heavy methylation of promoter regions and blocks the transcription of transgenes, and post-transcriptional gene silencing (PTGS), the basic mechanism is degradation of the cytosolic mRNA of transgenes or endogenous genes. Undesired transgene silencing is of major concern in the transgenic technologies used in crop improvement. A complete understanding of this phenomenon will be very useful for transgenic applications, where silencing of specific genes is required. The current status of epigenetic silencing in transgenic technology is discussed and summarized in this mini-review.

Pleiotropy of the Drosophila JAK pathway cytokine Unpaired 3 in development and aging.

The Janus kinase (JAK) pathway is an essential, highly re-utilized developmental signaling cascade found in most metazoans. In vertebrates, the JAK intracellular cascade mediates signaling by dozens of cytokines and growth factors. In Drosophila, the Unpaired (Upd) family, encoded by three tandemly duplicated genes, is the only class of ligands associated with JAK stimulation. Unpaired has a central role in activation of JAK for most pathway functions, while Unpaired 2 regulates body size through insulin signaling. We show here that the third member of the family, unpaired 3 (upd3), overlaps upd in expression in some tissues and is essential for a subset of JAK-mediated developmental functions. First, consistent with the known requirements of JAK signaling in gametogenesis, we find that mutants of upd3 show an age-dependent impairment of fertility in both sexes. In oogenesis, graded JAK activity stimulated by Upd specifies the fates of the somatic follicle cells. As upd3 mutant females age, defects arise that can be attributed to perturbations of the terminal follicle cells, which require the highest levels of JAK activation. Therefore, in oogenesis, the activities of Upd and Upd3 both appear to quantitatively contribute to specification of those follicle cell fates. Furthermore, the sensitization of upd3 mutants to age-related decline in fertility can be used to investigate reproductive senescence. Second, loss of Upd3 during imaginal development results in defects of adult structures, including reduced eye size and abnormal wing and haltere posture. The outstretched wing and small eye phenotypes resemble classical alleles referred to as outstretched (os) mutations that have been previously ascribed to upd. However, we show that os alleles affect expression of both upd and upd3 and map to untranscribed regions, suggesting that they disrupt regulatory elements shared by both genes. Thus the upd region serves as a genetically tractable model for coordinate regulation of tandemly duplicated gene families that are commonly found in higher eukaryotes.

Global Dissection of Alternative Splicing in Paleopolyploid Soybean

Alternative splicing (AS) is common in higher eukaryotes and plays an important role in gene posttranscriptional regulation. It has been suggested that AS varies dramatically among species, tissues, and duplicated gene families of different sizes. However, the genomic forces that govern AS variation remain poorly understood. Here, through genome-wide identification of AS events in the soybean (Glycine max) genome using high-throughput RNA sequencing of 28 samples from different developmental stages, we found that more than 63% of multiexonic genes underwent AS. More AS events occurred in the younger developmental stages than in the older developmental stages for the same type of tissue, and the four main AS types, exon skipping, intron retention, alternative donor sites, and alternative acceptor sites, exhibited different characteristics. Global computational analysis demonstrated that the variations of AS frequency and AS types were significantly correlated with the changes of gene features and gene transcriptional level. Further investigation suggested that the decrease of AS within the genome-wide duplicated genes were due to the diminution of intron length, exon number, and transcriptional level. Altogether, our study revealed that a large number of genes were alternatively spliced in the soybean genome and that variations in gene structure and transcriptional level may play important roles in regulating AS.