Field and Laboratory Evidence That Chlorpyrifos Exposure Reduced the Population Density of a Freshwater Snail by Increasing Juvenile Mortality.

Pesticides have been frequently detected in global freshwater ecosystems, but attempts to document changes in population dynamics of organisms upon exposure to pesticides, establish a causal relationship between exposure and population effects, and identify the key toxic events within individuals under natural field conditions remain rare. Here, we used a field survey, a reciprocal cross-transplant experiment, and a laboratory toxicity experiment to build a compelling case that exposure to the insecticide chlorpyrifos was responsible for differences in snail (Bellamya aeruginosa) densities in eastern (ELL) and western basins of Liangzi Lake in China. Our field survey and reciprocal cross-transplant experiment revealed significant differences in snail densities, juvenile percentage, survival, and relative telomere length (RTL) in the two basins. The insecticide chlorpyrifos detected in snail tissues was negatively correlated with snail densities, the percentage of juvenile snails, and RTL and had an extremely high risk quotient in ELL. In the laboratory experiment, tissue concentrations of chlorpyrifos detected in ELL were associated with reduced RTL and increased juvenile mortality in B. aeruginosa. These results support the hypothesis that chlorpyrifos exposure in ELL reduced the density of snails by reducing juvenile survival and, consequently, recruitment to the adult population.

Research note: differential heterosis of spent laying hens' carcass characteristics and meat quality in reciprocal crosses between White Leghorn and Beijing-You chickens.

Hybridization is used extensively in commercial layer production. However, heterosis for carcass performance and meat quality of spent laying hens remains unclear, especially under the trend of extended laying cycles. In this study, indigenous Beijing-You chickens (Y) and elite White Leghorn layers (W) were selected to generate purebreds (WW and YY) and reciprocal crosses (WY and YW). Data on traits including carcass compositions, meat quality, and main nutrients for breast muscle were collected when chickens were fed to 100 wk of age. Results showed that body weight (BW) and dressed weight for WY and YW with positive heterosis were significantly higher than WW (P < 0.05). YW had the heaviest breast and thigh of 232.28 g and 278.48 g, respectively. The abdominal fat weight for WY and YW were greatly higher than that for WW (P > 0.05). The yields of carcass compositions, including the dressed yield, half eviscerated yield, eviscerated yield, breast yield and thigh yield, did not differ among the four genetic groups (P > 0.05), except for the yield of abdominal fat. The largest heterosis differences appeared in breast weight (12.26% in YW vs. -0.46% in WY) and abdominal fat yield (15.26% in YW vs. 24.55% in WY). Although BW for crossbreds were similar, the specific parts of the carcass between them were different. For meat quality, WY had negative heterosis (P < 0.05) with the lowest lightness and yellowness, whereas YW had the completely opposite trend. Neither pH1h nor pH24h values had differences among purebreds and reciprocal crossbreds (P > 0.05). The drip loss and cooking loss were 4.01%-4.77% and 15.59%-21.31% respectively among the four genetic groups. The main nutrients of breast, including moisture, crude protein, intramuscular fat and unsaturated fatty acid, did not differ for purebreds and crossbreds (P > 0.05), except for saturated fatty acid. In general, the crossbreds even at the later laying period still showed divergent heterosis on carcass performance and meat characteristics. In view of the heterosis, Beijing-You chickens can be used as the sire line in the crossbreeding to improve carcass compositions of spent hens.

A multivariate gametic model for the analysis of purebred and crossbred data. An example between two populations of Iberian pigs.

Crossbreeding plays a pivotal role within pig breeding programmes, aiming to maximize heterosis and improve reproductive traits in crossbred maternal lines. Nevertheless, there is evidence indicating that the performance of reciprocal crosses between two genetic lines might exhibit variability. These variations in performance can be attributed to differences in the correlations between gametic effects, acting as either sire or dam, within purebred and crossbred populations. To address this issue, we propose a multivariate gametic model that incorporates up to four correlated gametic effects for each parental population. The model is employed on a data set comprising litter size data (total number of piglets born-TNB- and number of piglets born alive-NBA-) derived from a reciprocal cross involving two Iberian pig populations: Entrepelado and Retinto. The data set comprises 6933 records from 1564 purebred Entrepelado (EE) sows, 4995 records from 1015 Entrepelado × Retinto (ER) crosses, 2977 records from 756 Retinto × Entrepelado (RE) crosses and 7497 records from 1577 purebred Retinto (RR) sows. The data set is further supplemented by a pedigree encompassing 6007 individual-sire-dam entries. The statistical model also included the order of parity (with six levels), the breed of the service sire (five levels) and the herd-year-season effects (141 levels). Additionally, the model integrates random dominant and permanent environmental sow effects. The analysis employed a Bayesian approach, and the results revealed all the posterior estimates of the gametic correlations to be positive. The range of the posterior mean estimates of the correlations varied across different gametic effects and traits, with a range between 0.04 (gametic correlation between the paternal effects for purebred and the maternal for crossbred in Retinto) and 0.53 (gametic correlation between the paternal effects for purebred and the paternal for crossbred in Entrepelado). Furthermore, the posterior mean variance estimates of the maternal gametic effects were consistently surpassed those for paternal effects within all four populations. The results suggest the possible influence of imprinting effects on the genetic control of litter size, and underscore the importance of incorporating crossbred data into the breeding value predictions for purebred individuals.

Co-expression network analysis of environmental canalization in the ascidian Ciona.

BACKGROUND: Canalization, or buffering, is defined as developmental stability in the face of genetic and/or environmental perturbations. Understanding how canalization works is important in predicting how species survive environmental change, as well as deciphering how development can be altered in the evolutionary process. However, how developmental gene expression is linked to buffering remains unclear. We addressed this by co-expression network analysis, comparing gene expression changes caused by heat stress during development at a whole-embryonic scale in reciprocal hybrid crosses of sibling species of the ascidian Ciona that are adapted to different thermal environments. RESULTS: Since our previous work showed that developmental buffering in this group is maternally inherited, we first identified maternal developmental buffering genes (MDBGs) in which the expression level in embryos is both correlated to the level of environmental canalization and also differentially expressed depending on the species' gender roles in hybrid crosses. We found only 15 MDBGs, all of which showed high correlation coefficient values for expression with a large number of other genes, and 14 of these belonged to a single co-expression module. We then calculated correlation coefficients of expression between MDBGs and transcription factors in the central nervous system (CNS) developmental gene network that had previously been identified experimentally. We found that, compared to the correlation coefficients between MDBGs, which had an average of 0.96, the MDBGs are loosely linked to the CNS developmental genes (average correlation coefficient 0.45). Further, we investigated the correlation of each developmental to MDBGs, showing that only four out of 62 CNS developmental genes showed correlation coefficient > 0.9, comparable to the values between MDBGs, and three of these four genes were signaling molecules: BMP2/4, Wnt7, and Delta-like. CONCLUSIONS: We show that the developmental pathway is not centrally located within the buffering network. We found that out of 62 genes in the developmental gene network, only four genes showed correlation coefficients as high as between MDBGs. We propose that loose links to MDBGs stabilize spatiotemporally dynamic development.

Wolbachia Strain wGri From the Tea Geometrid Moth Ectropis grisescens Contributes to Its Host's Fecundity.

Members of the Wolbachia genus manipulate insect-host reproduction and are the most abundant bacterial endosymbionts of insects. The tea Geometrid moth Ectropis grisescens (Warren) (Lepidoptera: Geometridae) is the most devastating insect pest of tea plants [Camellia sinensis (L.) O. Kuntze] in China. However, limited data on the diversity, typing, or phenotypes of Wolbachia in E. grisescens are available. Here, we used a culture-independent method to compare the gut bacteria of E. grisescens and other tea Geometridae moths. The results showed that the composition of core gut bacteria in larvae of the three Geometridae moth species was similar, except for the presence of Wolbachia. Moreover, Wolbachia was also present in adult female E. grisescens samples. A Wolbachia strain was isolated from E. grisescens and designated as wGri. Comparative analyses showed that this strain shared multilocus sequence types and Wolbachia surface protein hypervariable region profiles with cytoplasmic incompatibility (CI)-inducing strains in supergroup B; however, the wGri-associated phenotypes were undetermined. A reciprocal cross analysis showed that Wolbachia-uninfected females mated with infected males resulted in 100% embryo mortality (0% eggs hatched per female). Eggs produced by mating between uninfected males and infected females hatched normally. These findings indicated that wGri induces strong unidirectional CI in E. grisescens. Additionally, compared with uninfected females, Wolbachia-infected females produced approximately 30-40% more eggs. Together, these results show that this Wolbachia strain induces reproductive CI in E. grisescens and enhances the fecundity of its female host. We also demonstrated that wGri potential influences reproductive communication between E. grisescens and Ectropis obliqua through CI.

Parent of origin gene expression in the bumblebee, Bombus terrestris, supports Haig's kinship theory for the evolution of genomic imprinting.

Genomic imprinting is the differential expression alleles in diploid individuals, with the expression being dependent on the sex of the parent from which it was inherited. Haig's kinship theory hypothesizes that genomic imprinting is due to an evolutionary conflict of interest between alleles from the mother and father. In social insects, it has been suggested that genomic imprinting should be widespread. One recent study identified parent-of-origin expression in honey bees and found evidence supporting the kinship theory. However, little is known about genomic imprinting in insects and multiple theoretical predictions must be tested to avoid single-study confirmation bias. We, therefore, tested for parent-of-origin expression in a primitively eusocial bee. We found equal numbers of maternally and paternally biased expressed genes. The most highly biased genes were maternally expressed, offering support for the kinship theory. We also found low conservation of potentially imprinted genes with the honey bee, suggesting rapid evolution of genomic imprinting in Hymenoptera.

Maternal effects shape the alternative splicing of parental alleles in reciprocal cross hybrids of Megalobrama amblycephala × Culter alburnus.

BACKGROUND: Maternal effects contribute to adaptive significance for shaping various phenotypes of many traits. Potential implications of maternal effects are the cause of expression diversity, but these effects on mRNA expression and alternative splicing (AS) have not been fully elucidated in hybrid animals. RESULTS: Two reciprocal cross hybrids following hybridization of Megalobrama amblycephala (blunt snout bream, BSB) and Culter alburnus (topmouth culter, TC) were used as a model to investigate maternal effects. By comparing the expression of BSB- and TC- homoeologous genes between the two reciprocal cross hybrids, we identified 49-348 differentially expressed BSB-homoeologous genes and 54-354 differentially expressed TC-homoeologous genes. 2402, 2959, and 3418 AS events between the two reciprocal cross hybrids were detected in Illumina data of muscle, liver, and gonad, respectively. Moreover, 21,577 (TC-homoeologs) and 30,007 (BSB-homoeologs) AS events were found in the 20,131 homoeologous gene pairs of TBF3 based on PacBio data, while 30,561 (TC-homoeologs) and 30,305 (BSB-homoeologs) AS events were found in BTF3. These results further improve AS prediction at the homoeolog level. The various AS patterns in bmpr2a belonging to the bone morphogenetic protein family were selected as AS models to investigate the expression diversity and its potential effects to body shape traits. CONCLUSIONS: The distribution of differentially expressed genes and AS in BSB- and TC-subgenomes exhibited various changes between the two reciprocal cross hybrids, suggesting that maternal effects were the cause of expression diversity. These findings provide a novel insight into mRNA expression changes and AS under maternal effects in lower vertebrates.

Effect of the genetic line and oleic acid-enriched mixed diets on the subcutaneous fatty acid composition and sensory characteristics of dry-cured shoulders from Iberian pig.

The aim of this study was to research the effect of the genetic background (Retinto, Torbiscal, and their reciprocal crosses) on the subcutaneous fatty acids and the sensory characteristics of dry-cured shoulders from Iberian pig, and also to investigate whether there is some interaction between genotype and diet composition when pigs are reared indoors, to obtain information to improve the selection strategies for purebred Iberian pig. The genetic background affected both the fatty acid composition (C17:0, C17:1 n-7, C18:3 n-3 and C20:0 were significantly different) and the sensory characteristics (marbling, lean fibrousness, and flavour intensity and persistence were significantly influenced), which indicates that they should be taken into account in the selection strategies for purebred Iberian pig. In a similar way, the genotype × diet composition interaction also should be taken into account when selecting a genetic line or cross to be fed indoors on a particular diet because of its repercussion on the sensory characteristics.

AIJ: joint test for simultaneous detection of imprinting and non-imprinting allelic expression imbalance.

Epigenetics is the study of heritable changes in gene expression or cellular phenotype caused by mechanisms other than changes in the underlying DNA sequence. Genomic imprinting is an epigenetically regulated process by which imprinted genes are expressed in a parent-of-origin-specific manner. It can be confounded with a phenomenon, allelic expression imbalance (AEI), which, in this paper, refers to asymmetric expression of the two alleles of a heterozygous subject at a single nucleotide polymorphism not caused by imprinting (non-imprinting AEI). Since existing methods in the literature are not amenable to distinguishing imprinting from non-imprinting AEI for data without replicates, we propose AIJ, a joint test for simultaneous detection of imprinting and non-imprinting AEI that accounts for potential confounding using RNA-seq data based on a reciprocal cross design. Through a simulation study, we show that AIJ is more powerful compared to two frequently used methods that do not account for confounding. To illustrate the practical utility of AIJ, we applied the method to a mouse dataset and identified genes with the imprinting effect and/or non-imprinting AEI phenomenon, with some already confirmed in an existing database. The results are also largely consistent with a study on human data for a set of orthologous genes, affirming earlier conclusion in the literature that non-imprinting AEI events are evolutionarily conserved.

Analysis of allele-specific expression using RNA-seq of the Korean native pig and Landrace reciprocal cross.

OBJECTIVE: We tried to analyze allele-specific expression in the pig neocortex using bioinformatic analysis of high-throughput sequencing results from the parental genomes and offspring transcriptomes from reciprocal crosses between Korean Native and Landrace pigs. METHODS: We carried out sequencing of parental genomes and offspring transcriptomes using next generation sequencing. We subsequently carried out genome scale identification of SNPs in two different ways using either individual genome mapping or joint genome mapping of the same breed parents that were used for the reciprocal crosses. Using parent-specific SNPs, allele-specifically expressed genes were analyzed. RESULTS: Because of the low genome coverage (~4x) of the sequencing results, most SNPs were non-informative for parental lineage determination of the expressed alleles in the offspring and were thus excluded from our analysis. Consequently, 436 SNPs covering 336 genes were applicable to measure the imbalanced expression of paternal alleles in the offspring. By calculating the read ratios of parental alleles in the offspring, we identified seven genes showing allele-biased expression (P &lt; 0.05) including three previously reported and four newly identified genes in this study. CONCLUSION: The newly identified allele-specifically expressing genes in the neocortex of pigs should contribute to improving our knowledge on genomic imprinting in pigs. To our knowledge, this is the first study of allelic imbalance using high throughput analysis of both parental genomes and offspring transcriptomes of the reciprocal cross in outbred animals. Our study also showed the effect of the number of informative animals on the genome level investigation of allele-specific expression using RNA-seq analysis in livestock species.

Glucose Tolerance and Plasma Non-Esterified Fatty Acid Levels in Chickens Selected for Low Body Weight, Red Junglefowl, and their Reciprocal Cross.

Responses of an individual to food deprivation, such as a 16-h fast, are complex, and are influenced by environmental and genetic factors. Domestication is an ongoing process during which adaptations to changing environments occur over generations. Food deprivation by their caretakers is less for domestic chickens than for their junglefowl ancestors. Unlike domestic chicken, the junglefowl adapted over generations to periods of food deprivation, which may be reflected in differences in metabolic responses to brief periods without food. Here, we compared the blood glucose and plasma levels of non-esterified fatty acids (NEFA) among four populations when deprived of feed for 16 h. The four populations included a domestic White Rock experimental line (LWS) maintained for generations under ad libitum feeding, adult red junglefowl (RJF), and a reciprocal cross of the lines. Although there were significant differences in adult (31-week) body weight between the RJF (683 g) and LWS (1282 g), with the weight of F1 crosses being intermediate, the amount of abdominal fat relative to body weight was similar for all populations. Patterns for blood glucose responses to a glucose bolus after a 16-h fast were similar for the initial and final points in the parental and cross populations. However, RJF reached their peak faster than LWS, with the reciprocal cross intermediate to the parental populations. Plasma NEFA concentrations were higher after the 16-h fast than in fed states, with no population differences for the fasting state. However, in the fed state, NEFA levels were lesser for LWS than for others, which was reflected further in percentage change from fed to fasted. This larger change in LWS suggests differences in mobilization of energy substrates and implies that during domestication or development of the LWS line, thresholds for responses to acute stressors may have increased.

CBCT image based segmentation method for tooth pulp cavity region extraction.

OBJECTIVES:: A method was proposed to segment the tooth pulp cavity region in cone beam CT) images, which aimed to make the extraction process more efficient and generate more reliable results for further research. METHODS:: Cone beam CT images of 50 teeth from 10 patients were randomly collected with the help of Peking University Hospital of Stomatology. All slice images have a ground truth tooth pulp cavity region delineated by two doctors manually. After necessary gamma transform in pre-processing stage, three kinds of information in an image such as greyscale, neighbour average greyscale and gradient were fused to search an optimal segmentation threshold by using plane intercept histogram of reciprocal cross entropy algorithm. With the optimal threshold, binarization was conducted and the tooth pulp cavity regions in slice images can be extracted. Qualitative and quantitative analyses compared to ground truth are involved with the evaluation criterion of average non-coincidence rate ( RANOA ). Independent repeated experiments were carried out to test the stability of this segmentation method. RESULTS:: Accurate and complete segmentation results are obtained. The proposed method reaches the lowest RANOA values in most cases and owns more competitive robustness under various interferences compared with the other popular segmentation methods like reciprocal cross entropy method, active contour-based method, region growing method and level set method. Quantitative analysis verified the effectiveness of this method. CONCLUSIONS:: The proposed method can extract tooth pulp cavity regions from teeth efficiently. The segmentation results of this method are more accurate compared to other popular methods under different circumstances and can be used for subsequent applications.

Genetic analysis assessed by microsatellites for a diallel mating design of two geographical stocks of the blood clam Tegillarca granosa.

To determine the potential for productive efficiency and genetic improvement in the blood clam Tegillarca granosa, four offspring populations (ZZ, ZK, KZ and KK) were produced from a diallel mating of two different geographical stocks (Z and K). The levels of genetic diversity and population structures of four populations were analyzed using 14 polymorphic microsatellites. The results showed that the mean observed heterozygosities (Ho) of reciprocal cross populations (ZK and KZ) was higher than those of pure populations(ZZ and KK). The largest values of genetic differentiation coefficient (Fst = 0.067) and Nei's unbiased genetic distance (Dc = 0.263) were between ZK and KZ, and the smallest (Fst = 0.020, Dc = 0.116) were between ZZ and KK, which revealed that the largest genetic divergence was between the two reciprocal cross populations and the smallest was between two pure populations. This study demonstrated that the reciprocal cross populations of T. granosa had an extensive genetic difference and improvement, which may be advantageous for future breeding studies.

Quantifying pearl millet response to high temperature stress: thresholds, sensitive stages, genetic variability and relative sensitivity of pollen and pistil.

The objectives were to (1) quantify high temperature (HT) stress impacts at different growth stages (season long, booting to seed-set and booting to maturity) on various yield components; (2) identify the most sensitive stage(s) to short episodes of HT stress during reproductive development; (3) understand the genetic variations for HT stress tolerance based on cardinal temperatures for pollen germination; and (4) determine relative sensitivity of pollen and pistil to HT stress and associated tolerance or susceptible mechanisms in pearl millet. High temperature stress (≥36/26°C) imposed at different stages and durations caused decrease in number of seeds, individual seed weight and seed yield. Two periods (10-12 days and 2-0 days before anthesis) were identified as most sensitive to short episodes of stress, causing maximum decreases in pollen germination percentage and seeds numbers. HT stresses of ≥36/26°C results in floret sterility. Pistils were relatively more sensitive than pollen grains, causing decreased number of seeds and seed yield. HT stress increased the reactive oxygen species contents and decreased the activity of the antioxidant enzymes in both pollen and pistils. Under HT stress, pistils had relatively higher reactive oxygen species and lower antioxidant enzymes activity compared with pollen grains, which explains greater susceptibility of pistils.

Is Longevity a Heritable Trait? Evidence for Non-genomic Influence from an Extended Longevity Phenotype of Drosophila melanogaster.

BACKGROUND: Although genetic variations are heritable, some quantitative traits like longevity may have non-genomic influence on heritability. Laboratory-selected inbred strains of extended longevity phenotype of Drosophila offer an opportunity to study the inheritance of longevity. OBJECTIVE: The aim of the study was to examine the heritability of longevity in an extended longevity phenotype of Drosophila melanogaster using reciprocal cross effects in F1 and F2 generations. METHODS: Lifespan variations of virgin and mated flies in parent, F1 and F2 generations were investigated using reciprocal crosses between normal and long lifespan lines of inbred population of D. melanogaster. Heterosis, narrow-sense heritability, recombination loss, maternal effect and overdominance with respect to survivorship in virgin and mated flies were analyzed. RESULTS: Virgin flies lived longer than mated flies. There was no significant effect of mid-parent heterosis, recombination loss and overdominance on variations in longevity, whereas, significant maternal effect and narrow-sense heritability were observed in mated and virgin flies, respectively. CONCLUSION: Absence of heterosis in our study population of Drosophila phenotypes could be due to the lack of genetic heterogeneity. The heritability of the longevity trait in an inbred extended longevity phenotype depends on the variations in genetic and environmental factors.

Reciprocal cross-regulation of VND and SND multigene TF families for wood formation in Populus trichocarpa.

Secondary cell wall (SCW) biosynthesis is the biological process that generates wood, an important renewable feedstock for materials and energy. NAC domain transcription factors, particularly Vascular-Related NAC-Domain (VND) and Secondary Wall-Associated NAC Domain (SND) proteins, are known to regulate SCW differentiation. The regulation of VND and SND is important to maintain homeostasis for plants to avoid abnormal growth and development. We previously identified a splice variant, PtrSND1-A2IR , derived from PtrSND1-A2 as a dominant-negative regulator, which suppresses the transactivation of all PtrSND1 family members. PtrSND1-A2IR also suppresses the self-activation of the PtrSND1 family members except for its cognate transcription factor, PtrSND1-A2, suggesting the existence of an unknown factor needed to regulate PtrSND1-A2 Here, a splice variant, PtrVND6-C1IR , derived from PtrVND6-C1 was discovered that suppresses the protein functions of all PtrVND6 family members. PtrVND6-C1IR also suppresses the expression of all PtrSND1 members, including PtrSND1-A2, demonstrating that PtrVND6-C1IR is the previously unidentified regulator of PtrSND1-A2 We also found that PtrVND6-C1IR cannot suppress the expression of its cognate transcription factor, PtrVND6-C1PtrVND6-C1 is suppressed by PtrSND1-A2IR Both PtrVND6-C1IR and PtrSND1-A2IR cannot suppress their cognate transcription factors but can suppress all members of the other family. The results indicate that the splice variants from the PtrVND6 and PtrSND1 family may exert reciprocal cross-regulation for complete transcriptional regulation of these two families in wood formation. This reciprocal cross-regulation between families suggests a general mechanism among NAC domain proteins and likely other transcription factors, where intron-retained splice variants provide an additional level of regulation.

[Pre-programmed genes]. / Eloprogramozott gének.

Cells feel good and carry on perfect functions when they contain the right types of proteins in the right concentration, at the right time and sites. There are mechanisms that ensure the right level of gene expression in the different cell types: the formation of protein molecules based on the DNA-encoded genetic information. Gene expression can also be regulated through the compactness of chromatin, i.e. the accessibility of the genes. The chromosomes are repositories of the genetic information - the sequence of base pairs - and also of the so-called epigenetic mechanisms that control gene expression through the regulation of chromatin compactness. The epigenetic mechanisms operate through DNA methylation and/or the regulation of chromatin compactness. The present overview takes a look into the phenomenon of epigenesis. It summarizes how genetic crosses reveal the involvement of epigenesis, explains its meaning and impact on life of the organisms. An understanding of epigenesis provides guidance to improve our life. Orv Hetil. 2017; 158(34): 1323-1330.

Temporal small RNA transcriptome profiling unraveled partitioned miRNA expression in developing maize endosperms between reciprocal crosses.

In angiosperms, the endosperm nurtures the embryo and provides nutrients for seed germination. To identify the expression pattern of small interfering RNA in the developing maize endosperm, we have performed high-throughput small RNA transcriptome sequencing of kernels at 0, 3, and 5 days after pollination (DAP) and endosperms at 7, 10, and 15 DAP using B73 and Mo17 reciprocal crosses in previous study. Here, we further explored these small RNA-seq data to investigate the potential roles of miRNAs in regulating the gene expression process. In total, 57 conserved miRNAs and 18 novel miRNAs were observed highly expressed in maize endosperm. Temporal expression profiling indicated that these miRNAs exhibited dynamic and partitioned expression patterns at different developmental stages between maize reciprocal crosses, and quantitative RT-PCR results further confirmed our observation. In addition, we found a subset of distinct tandem miRNAs are generated from a single stem-loop structure in maize that might be conserved in monocots. Furthermore, a SNP variation of Zma-miR408-5p at 11th base position was characterized between B73 and Mo17 which might lead to completely different functions in repressing targets. More interestingly, Zma-miR408-5p exhibited B73-biased expression pattern in the B73 and Mo17 reciprocal hybrid endosperms at 7, 10, and 15 DAP according to the reads abundance with SNPs and CAPS experiment. Together, this study suggests that miRNA plays a crucial role in regulating endosperm development, and exhibited distinct expression patterns in developing endosperm between maize reciprocal crosses.

A reciprocal cross design to map the genetic architecture of complex traits in apomictic plants.

Many higher plants of economic and biological importance undergo apomixis in which the maternal tissue of the ovule forms a seed, without experiencing meiosis and fertilization. This feature of apomixis has made it difficult to perform linkage mapping which relies on meiotic recombination. Here, we describe a computational model for mapping quantitative trait loci (QTLs) that control complex traits in apomictic plants. The model is founded on the mixture model-based likelihood in which maternal genotypes are dissolved into two possible components generated by meiotic and apomictic processes, respectively. The EM algorithm was implemented to discern meiotic and apomictic genotypes and, therefore, allow the marker-QTL linkage relationship to be estimated. By capitalizing on reciprocal crosses, the model is renovated to estimate and test imprinting effects of QTLs, providing a better gateway to characterize the genetic architecture of complex traits. The model was validated through computer simulation and further demonstrated for its usefulness by analyzing a real data for an apomictic woody plant. The model has for the first time provided a unique tool for genetic mapping in apomictic plants.